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Preface 


The  proceedings  of  the  Resonance  Meeting  consist  of  two  volumes.  Volume  1  contains  a 
verbatim  transcript  of  each  presentation.  Volume  2  contains  accompanying  transparencies  for  each 
transcript.  Anyone  interested  in  obtaining  more  information  on  a  particular  presentation  may 
contact  the  author(s).  Volume  1  contains  a  complete  list  of  participants  and  addresses. 


Copies  of  these  proceedings  may  be  obtained  from  the 
National  Center  for  Physical  Acoustics 
University  of  Mississippi 
University,  MS  38677 
601-232-5889 
NCPA@OLEMISS.EDU 
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1.  Howto  do  it 

2.  What  it  does 

3.  Non-destructive  testing 
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Why  study  ultrasound? 

^<3  =  Gfto,  -  GeM 

AG  is  the  Gibbs  free  energy 
change  across  a  phase  boundary 

Al/  =  -^AG 

A\/  is  the  volume  change  across 
the  boundary,  and  is  zero  for  second 
order  phase  transitions 


_1 _ 1_ 

^hot  ^cold 


I  _  1  ^AG 

V  ^  ~  V  dP^ 


Ax:  is  the  discontinuity  in  bulk  modulus 


across  the  boundary,  which  ties  directly 


to  the  speed  of  sound,  and  which  is  not 


zero  for  a  second  order  transition 
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gijnai/nnifi«^.cQm p  prison  of  pulsed  and  resonant  measurements 


parameter 

Ihipulse 

Swept  Sine 

drive  power  per  unit 
bandwidth 

peak  power/full 
bandwidth 
10®/10®=.001 

peak  power/sweep 
rate 

1/100=0.01 

noise  bandwidth  for 
complete 

measurement  using 
optimum  receiver 

10® 

number  of  modes  x  . 
width  of  each  mode  x 
10=10^Hz 

drive  duty  cycle 
(typical) 

10'^ 

1 

detect  duty  cycle 

1 

1 

square  root  of  all 
factors,  which  is  a 
measure  of  S/N 

3x10-' 

10'^ 

Table  I.  Signal-to-noise  comparison  between  impulse  (pulse-echo)  and 
swept-sine  (RUS)  resonance  measurement  methods  for  a  measurement  of 
a  1  cm  sample  with  resonances  having  a  Q=10^  using*  10  modes  over  0.5 
MHz-1. 5MH2=10®  Hz  bandwidth  to  obtain  an  elastic  modulus.  Note  that  the 
pulse-echo  measurement  provides  about  0.1%  absolute  accuracy  at  best, 
compared  with  about  0.01%  for  the  best  RUS  measurements. 
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John  William  Strutt,  the  Baron  Rayleigh  tried  to  do  this 
computation.  Without  a  90MHz  Pentium,  he  found  that 


In  the 

cftse  of  A  short  rod  end  of  a  perticle  ntUAted  neer  the  ^lindricAl 
boundaiy,  this  lAteral  motion  would  be  eompeiAble  in  msgnitude 
with  the  longitudinAl  motion,  And  could  not  be  overlooked  without 
risk  of  considerable  error. 

226.  The  problem  of  a  rectAngular  pUte,  whose  edges  are 
free,  is  one  of  greet  di£5eulty,  and  has  for  the  most  part  resisted 
Attack*. 


Even  with  a  Pentium,  if  you  try  this  using  finite  element  methods, 
the  computation  time  goes  like  the  cube  of  the  numerical  accuracy 
and  you  can’t  compute  as  well  as  you  can  measure  in  a 
reasonable  time  on  a  reasonable  computer.  However,  if  you  are 
careful,  and  smart,  as  were  Orson  Anderson  and  his  postdoc 
Harold  Demarest  at  Bell  labs  30  years  ago,  then . 
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1 2  ELASTIC  MODULI  AND  SOLID-BODY  RESONANCES 
3.3  The  problem  with  3-D  resonances  in  soiids 

As  suggested  above,  the  resonances  of  3-dimensjonai  solids  with  small  aspect 
ratios  (all  the  dimensions  are  comparable)  are  very  difficult  to  compute.  To  see 
why,  let’s  set  up  a  problem  that  John  William  Strutt,  Baron  Rayleigh^  worked 
on  his  honeymoon.  We  will  not  solve  the  problem  here,  nore  we  will  follow 
through  with  Rayleigh’s  approximate  solution  because,  as  Rayleigh  notes'’A 
more  complete  solution.. .has  been  given  by  Pochhammer..."®.  Furthermore,  in 
Chapter  4,  we  will  set  up  all  the  machinery  to  solve  this  and  other  nasty 
problems  to  as  high  an  accuracy  as  anyone  might  care  about  (6  digits  is  easy 
on  a  PC). 

The  problem  here  is  to  redo  the  computation  that  resulted  in  2.2.7  but 
now  with  a  short  fat  rod.  When  such  a  rod  resonates,  it  alternately  bulges  and 
constricts  at  some  places  along  its  length  but  not  at  others. 


Fig.  3.3  A  short  fat  bar  undergoing  a  longitudinal  oscillation.  As  the  bar 
oscillates  between  being  longer  and  shorter,  it  also  alternately  bulges  and 
constricts  near  its  middle. 

Stresses  along  the  bar  associated  with  longitudinal  oscillations,  combined  with 
the  effect  of  the  Lam6  coefficient  A,  cause  lateral  strains  as  depicted  in  Fig. 
3.3.  The  problem  is  to  understand  completely  all  the  strains  and  velocities  of 
all  the  material  in  the  bar.  This  is  very  difficult  to  do  simply.  What  we  can  do 
is  use  o  to  estimate  the  lateral  motion  of  the  bar,  include  this  in  the  total  kinetic 
energy,  and  then  estimate  the  frequency  shift.  This  approach  relies  on  an 
oscillator  alternating  between  brief  moments  when  everything  is  stopped  (think 
of  the  top  of  a  pendulum’s  swing)  to  brief  moments  when  ail  that  moves  is 
moving  at  its  peak  speed  and  no  strians  are  present.  Because  we  are  not 
considering  any  dissipation  to  be  present,  the  total  kinetic  energy  at  peak 
speeds  equals  the  total  potential  energy  when  things  stop  briefly.  If  in  the  case 
of  a  short  fat  bar,  more  mass  moves  than  in  a  long  thin  bar  of  eqaul  mass,  then 
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Properties  and  Composition  of  Lunar  Materials#  Earth  Analofies 


Abstnet  The  sound  velocity  data  for  the  lunar  rocks  were  compared  $o  numer^ 
ous  terrestrial  rock  types  and  were  found  to  deviate  widely  from  them,  A  group  of 
terrestrial  materials  were  found  which  have  velocities  comparabUt  to  those  of  the 
lunar  rocks,  but  they  do  obey  velocity^densiry  relations  proposed  for  earth  rocks. 


Certain  dau  from  Apollo  11  and 
Apollo  12  missions  present  some  dilB* 
mlties  in  that  they  require  explanations 
for  the  signals  received  by  the  lunar 
seismograph  as  a  result  of  the  impact 
of  the  lunar  module  (LEM)  on  the  lunar 
surface  (/).  In  particular,  the  observed 
signal  does  not  resemble  one  due  to  an 
impulsive  source,  but  exhibits  a  geoer- 
ally  slow  build-up  of  energy  with  time, 
in  spite  of  the  appearance  of  the  re* 
turned  lunar  samples,  the  lunar  seismic 
signal  continued  to  ring  for  a  remark* 
ably  long  time — a  characteristic  of  very 
high  Q  material.  The  lunar  rocks,  when 
studied  in  the  laboratory,  exhibited  a 
low  Q  (2).  Perhaps  most  startling  of  all. 


however,  was  the  very  low  sound  veloc¬ 
ity  indicated  for  the  oiHer  lunar  layer 
deduced  from  the  LEM  impact  signal. 
The  dau  obtained  on  the  lunar  rocks 
and  fines  agree  well  with  the  muKs  of 
the  Apollo  12  seismic  experiment  (2, 
J).  These  rock  velocities  are  surtlingly 
low.  The  measured  velocities  on  a  vesic¬ 
ular  medium  grained,  igneous  rock 
(10017)  having  a  bulk  density  of  3.2 
gfesn^  were  v,  =  1.S4,  and  v,  =  1.05 
km/sec.  Hie  results  for  a  microbreeda 
(10046)  with  a  bulk  density  of  2.2  gf 
cm*  were  Vp  =  1.25  and  v.  =  0.74  km/ 
sec  for  the  compressional  (v,)  and 
shear  (v.)  velocities. 

It  was  of  some  interest  to  consider 


the  behavior  of  these  lunar  rocks  in 
terms  of  the  expected  behavior  based 
on  measurements  of  earth  materials. 
Birch  (4)  first  proposed  a  simple  linear 
relation  between  compressional  velocity 
and  density  for  rocks.  This  relation  was 
examined  further  by  Anderson  (5)  who 
showed  that  this  was  a  first  approxima¬ 
tion  to  a  more  general  relaiion,  deriva- 
Me  from  a  dependence  of  the  elastic 
moduli  with  the  density  through  a 
power  function.  Comparison  of  the  re¬ 
sults  obtained  from  the  returned  lunar 
rocks  with  the  predictions  of  these 
relationships  expresses  graphically  the 
manner  they  Aviate  from  the  be¬ 
havior  of  rocks  found  on  earth.  The 
velocities  are  remarkably  lower  than 
what  would  be  predicted  from  either 
the  Birch  or  Anderson  relationships. 

To  account  for  this  very  low  velocity, 
we  decided  to  consider  maurials  other 
than  those  listed  initially  by  Birch  (4) 
or  more  detailed  compilation  of  Ander¬ 
son  and  liebermann  (6).  The  search 
wu  aided  by  considerations  of  much 
earlier  speculations  concerning  the  na- 


Tsble  1.  Comparison  of  eorapmskmal  vclockics  of  lunar  rocks  and  various  earth  materials. 
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Computation  of  resonances 
from  Migliori  et.al.  Physica  B  183,1,1993 


The  procedure  for  solving  the  direct  problem 
for  an  arbitrarily  shaped  elastic  solid  with  vol¬ 
ume  V,  elastic  tensor  density  p,  and  with  a 
free  surface  S  begins  with  the  Lagrangian 

L  =  J  (KE  -  PE)  dV  (4) 

V 

where  the  kinetic  energy,  KE,  is  given  by 

KE=  ,  (5) 


and  the  potential  energy,  PE,  by 


Following  Hamilton,  we  allow  «,  to  vary  arbit¬ 
rarily  in  the  volume  V  and  on  the  surface  5 
(u,  — >  «,  +  8m,)  and  calculate  the  variation  8L  in 

L.  The  result  is 

8L  =  J  (left  side  of  eq.  (8)),-8m,  dV 

V 

+  f  (left  side  of  eq.  (9)), 8m,-  d5  (7) 

s 

The  immediate  results  are  two  equations,  the 
elastic  wave  equation  and  the  vanishing  of 

surface  traction 

+  ,*/“*.// =  0,  (8) 


,  =  0  (9) 

where  {nj  is  the  unit  outer  normal  to  S. 

Because  of  the  arbitrariness  of  8w.  in  V  and  on 
S,  the  u/s  which  correspond  to  stationary  points 
of  L  (i.e.  8L  =0)  must  satisfy  eq.  (8)  in  V  and 
eq.  (9)  on  5.  There  are  no  such  m/s,  of  course, 
unless  is  one  of  a  discrete  set  of  eigenvalues, 
the  normal  mode  frequencies  of  free  vibration  of 
the  system. 


10 
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Following  the  Rayleigh-Ritz  prescription,  we 
expand  the  displacement  vector  in  a  complete  set 
of  functions 

’  (10) 

and  choose  as  our  basis  functions  powers  of 
cartesian  coordinates: 

^.=x'y'^z\  (11) 

where  A  =  (/,  m,  n)  is  the  function  label,  a  set  of 

three  nonnegative  integers.  After  substituting 

eq.  (10)  into  eq.  (4),  we  obtain  (c  becomes  a 
column  vector) 

L  =  (12) 


11 


TRANSPARENCY  11 


Vibrations  of  a  Rectangular  Paraflefipiped 
*Not  a  plane  wave  among  them 
•Sufficiently  complex  to  provide  all  the  elastic  moduli 
•Completely  understood  mechanics  problem 
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The  “Calvin  and  Hobbes”  mode!  of  the  vibrations  of  a  rectangular 

parallelepiped 
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FIG.  8.  Frequency  spectra  of  a  number  of  objects  in  the  potato  family.  The 
seven  stations  correspond  to  shapes  as  labeled,  with  semiaxes  as  given  in 
Table  II.  The  sphere  frequencies  agree  well  with  those  in  the  literature  for 
these  material  parameters  (Poisson’s  ratio  =  1/4).^°  The  dimensional  pa¬ 
rameters  are  interpolated  linearly  between  the  seven  sta¬ 

tions  here.  Several  interesting  features  invite  comment.  First,  the  potato  has 
no  degenerate  lines,  because  of  its  low  symmetry,  and  the  sphere,  converee- 
ly,  has  few  lines  that  are  nondegenerate.  The  ellipsoid  has  no  degraeracies, 
and  the  spheroid,  the  egg,  and  the  hemisphere  (all  being  rotationally  sym¬ 
metric  )  do,  but  never  more  than  doubly  degenerate  lines.  Small  deviations 
from  the  sphere  in  the  egg  direction  do  not  change  any  of  the  frequencies  to 
first  order,  because  ^  increases  as  much  decreases,  compensating 

one  another  as  far  as  affecting  resonant  frequencies  is  concerned.  As  in  sev¬ 
eral  other  figures,  apparent  avoided  crossings  on  this  plot  should  be  viewed 
with  suspicion  because  the  plotting  program  does  not  interchange  line  iden¬ 
tities  when  physically  the  modes  do,  in  fact,  cross.  Spectra  are  computed  for 
24 1  abscissa  values  here  and  elsewhere,  which  sets  the  scale  on  which  avoid¬ 
ed  crossings  may  be  spurious. 
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Instrumentation  Block  Diagram 


Transducers 

Transmitting  Receiving 


Low  temperature  RUS  cell 
with  a  2mm  rectangular 
parallelipiped  sample  set 
between  2  diamond/LiNbOa 
transducers  using  no 
coupling  fluids.  The  force  on 
the  sample  is  equivalent  to  1 
gm  weight.  This  cell 
operates  between  1K  and 
400K. 
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Transducer  construction 


Lithium  Niobate  Disc 


Diamond 
Cylinder  ^ 


1.5  mm 


Silvered  Kapton  Film 


Kapton 
Electrical 
,  Lead 


Fig.  3.  Shown  is  a  schematic  of  the  diamond /polyimide/ 
LiNbO,  composite  transducer  used  for  all  the  measurements. 
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Resonances  of  a  5120  steel  RPR  at  38C  and  378C 
measure  using  metal  diffusion  bonded  LiNbOj/Alumina 
transducers 


38C 


25 

20 


1.5 


1.0 


0.5 


0.0 


5120ated378C 


f 


r  1 


JL 


730000 


740000 


750000 


760000 


Frequency  (ft) 


378C 
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Ultimate  accuracy  determined  by  geometry-for  a  Si3N4  ball 
bearing,  geometry  errors  are  less  that  1  part  In  10®.  So  are  the 
modulus  errors! 


Tabic  1 

Resonant  ultrasound  measurement  of  a  0.63500  cm  diameter  Si,N,  ceramic  sphere  with  a  density  of  3.2325  g/cm\  /«  are  measured 
frequencies.  /,  are  fitted,  n  is  the  mode  number,  k  is  our  designator  (to  be  discussed  below)  for  the  symmetry  of  the  mode  and  i  is 
in  essence  the  harmonic  number  of  each  symmetry  type.  Multiple  entries  indicate  the  mode  degeneracy.  The  fit  for 
^  =  1.2374x  10‘-dyne/cm"  and  o*  =  0.2703  ha.s  a  x'  (%)  =  0.0124.  This  is  sufficient  to  determine  m  to  about  0.01%  and  cr  to 
about  0.05% .  There  are  no  corrections  so  these  values  are  absolute. 


n 

f,  (MHz) 

L  (Mto) 

%  error 

(*.i) 

1 

0.775706 

0.775707 

-0.000138 

(6.1).(l.l).(4.  l).(4.2),(7,l) 

6 

0.819567 

0.819983 

-0.050778 

(5. 1).  (3.1).  (5.  2).  (8.1).  (2.  n 

11 

1.075664 

1.075399 

0.024614 

(1.2).  (7. 2).  (6.  2) 

14 

1.198616 

1.198505 

0.009239 

(5. 3).  (2. 2).  (3.  2).  (8. 2).  (3. 3).  (8. 3).  (2. 3) 

21 

1.217375 

1.217850 

-0.039042 

(1.3).  (6. 3).  (7.  3).  (1. 4).  (6. 4).  (7. 4).  (4. 3) 

28 

1.440760 

1.440750 

0.000712 

(5.4) 

29 

1.527080 

1.526474 

0.039695 

(5. 5).  (8. 4).  (3.  4).  (5. 6),  (2. 4) 

34 

1.558358 

1.558848 

-0.031448 

(5. 7).  (5. 8).  (5.  V).  (3. 5).  (8. 5).  (2. 5).  (3, 6).  (8. 6).  (2. 6) 

43 

1.580067 

1.579871 

0.012426 

(6. 5).  (7. 5).  (7.  6).  (1. 5).  (4. 4).  (1.6).  (6. 6).  (4. 5).  (4. 6) 
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Determination  of  high  temperature  moduii  to  support  the 
broad  DOE  initiative  in  heat  treatment  distortion 


5120 steel  -ground parallelepiped 

free  dimensions  are  d1,  d2,  d3  (initial  0.41205,  0.33973,  0.26380)  cm 
using  10  order  polynomials  mass=  .2875  gm;  denslty=  7.785  gm/cc 


n 

f-expt 

f-calc 

%err 

C11  C44 

df/d(moduli) 

1 

.340010 

.339607 

-.12 

.00 

1.00 

2 

.447200 

.447109 

-.02 

.13 

.87 

3 

.490240 

.490366 

.03 

.16 

.84 

4 

.541860 

.541482 

-.07 

.02 

.98 

5 

.557460 

.556663 

-.14 

.00 

1.00 

6 

.591740 

.592306 

.10 

.05 

.95 

7 

.608670 

.608417 

-.04 

.31 

.69 

8 

.615590 

.615290 

-.05 

.02 

.98 

9 

.623100 

.622662 

-.07 

.07 

.93 

10 

.642700 

.642897 

.03 

.12 

.88 

11 

.675510 

.676384 

.13 

.10 

.90 

12 

.683560 

.683316 

-.04 

.05 

.95 

13 

.690650 

.690436 

-.03 

.16 

.84 

14 

.746720 

.747335 

.08 

.10 

.90 

15 

.753860 

.754125 

.04 

.05 

.95 

16 

.827430 

.827569 

.02 

.06 

.94 

17 

.848300 

.849153 

.10 

.03 

.97 

18 

.855320 

.854953 

-.04 

.08 

.92 

19 

.870540 

.870865 

.04 

.12 

.88 

20 

.878390 

.878891 

.06 

.11 

.89 

21 

.882080 

.881734 

-.04 

.33 

.67 

22 

.884870 

.884342 

-.06 

.21 

.79 

23 

.887380 

.887639 

.03 

.20 

.80 

24 

.891190 

.891892 

.08 

.04 

.96 

25 

.916680 

.916689 

.00 

.31 

.69 

rms  errors  .0680  %, 

Fitted  values  for  elastic  constants  and  dimensions: 

Young’s  Modulus  =  30,39x10*  psi  (209.55  GPa) 

Shear  Modulus  s  11.82x10*  psi  (81.50  GPa)  Poisson’s  ratio- 0.285 
d1  =  0.4121  cm  d2  =  0.3395  cm  d3  *  0.2639  cm 

Perturbation  analysis  of  error  sensitivity  :  (goodness  of  fit) 

X*  increased  2%  by  the  following  %  changes  in  independent  parameters 
C11  C44  d1  d2  d3 

.48  .41  -.42  -.42  -.4  =>  errors  on  fitted  vaiues  -  0.5% 
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Elastic  constants  of  copper.  B  represents  bulk  modulus. 


Single 

crystal 

Polycrystal  (wire-drawn) 

Cij  (GPa) 

Literature 

average 

Measured® 

RP 

Cylinder 

Cll 

168.75 

170.88 

193.61 

194.25 

C33 

205.88 

203.98 

C12 

122.14 

124.63 

105.65 

106.84 

Cl3 

95.00 

95.93 

C44 

75.48 

74.01 

39.35 

39.46 

^66 

43.98 

43.71 

B 

137.68 

140.05 

131.65 

132.84 

“Crystal  rotated  45.4®  about  [100]. 
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Elastic  constants  of  tantalum  at  room  temperature. 


T(K) 

P  (g/cm^) 

Cii  (GPa) 

ci2  (GPa) 

C44  (GPa) 

300‘ 

16.678 

266.7 

160.8 

82.5 

0 

0 

CO 

16.678 

266.8 

161.4 

82.5 

3002 

16.633 

260.9 

157.4 

81.8 

298^ 

16.626 

260.2 

154.5 

82.6 

295^ 

16.641 

266.3 

160.5 

82.8 

iD.I.  Bolef,  J.  Appl.  Phys.  33,  2311  (1962). 

^F.H.  Featherstone  and  J.R.  Neighbours,  Phys.  Rev.  130,  1324  (1963). 

^N.  Soga,  J.  Appl.  Phys.  37,  3416  (1966). 

"^Euler  angles  (defined  in  Roe’s  convention)  were  determined  to  be  a  =  138.3®, 
/?  =:  29.7®,  and  7  =  155.1®  by  RUS  and  a  =  135®,  /?  =  33®,  and  7  =  158®  by 
X-ray. 
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Structural  phase  transition  in  La2.xSrxCu04 

RUS  simultaneously  sees  collapse  of  Cqq 
and  no  effect  on  C44  in  a  2mm  single  crystal 


220 


0.  7 

cTO.  61- 

£ 

^0.  5 

c 

^0-  4 

T3 
C\? 

•^0.3F 

Xo. 

^  0.  1 

0.0  _ 
220 


240  260  280 

Temperature(K) 


!•••  •  •  • 


300 


L32Cu  O4 


240  260  280 

Temperature  (K) 


300 
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P 

nn,  np,  it 
i  fax 

,  rho*  6, 
fr 

10,  1 

%err 

,  6.946 

wt  k 

i 

11 

33  23  12 

44 

00 

1 

0.550101 

0.548752 

-0.25 

0.0 

4 

1 

0.01 

0.01  0.00  0.00 

0.17 

0.33 

2 

0.055544 

0.053127 

rO.37 

0.0 

4 

2 

0.01 

0.00  0.00  0.00 

0.47 

0.02 

3 

0.707270 

0.707542 

0.03 

1.0 

3 

1 

0.02 

0.01  -0.01  0.00 

0.00 

0.48 

4 

0.840904 

0.840137 

-0.09 

1.0 

1 

1 

0.35 

0.17  -0.18  0.01 

0.01 

0.15 

5 

0.842859 

0.840724 

-0.25 

1.0 

8 

1 

0.02 

0.01  0.00  0.00 

0.48 

0.00 

0 

0.870588 

0.872288 

,  0.20 

1.0 

2 

0.02 

0.01  -0.01  0.00 

0.48 

0.00 

7 

0.882324 

0.882338 

0.00 

1.0 

0 

1 

0.41 

0.03  -0.05  -0.00* 

0.17 

0.00 

8 

0.914348 

0.910208 

0.21 

1.0 

7 

0.31 

0.27  -0.24  0.02 

0.04 

0.10 

9 

0.982734 

0.985644 

0.32 

1.0 

0 

2 

0.35 

0.00  -O.Il  0.01 

0.19 

0.00 

10 

1.017951 

1.018350 

0.04 

1.0 

5 

1 

0.40 

0.45  -0.41  0.00 

0.00 

0.00 

11 

^.031479 

1.028863 

-0.25 

1.0 

1 

2 

0.07 

0.30  -0.09  0.00 

0.14 

0.09 

12 

1.004057 

1.002290 

-0.17 

1.0 

2 

2 

0.30 

0.04  -0.09  -0.01 

0.19 

0.01 

13 

1.070904 

1.071351 

0.04 

1.0 

7 

2 

0.10 

0.19  -0.04  -0.01 

0.12 

0.14 

14 

1.071533 

1.072243 

0.07 

1.0 

5 

2 

0.47 

0.31  -0.27  -0.01 

0.00 

0.00 

15 

1.091915 

1.091434 

-0.04 

1.0 

3 

2 

0.09 

0.45  -0.21  0.02 

0.04 

0.12 

10 

1.101242 

1.101298 

0.01 

1.0 

2 

3 

0.05 

0.02  -0.02  0.00 

0.21 

0.24 

17 

1.119834 

1.110059 

-0.28 

1.0 

8 

2 

0.37 

0.05  -0.08  -0.02 

0.18 

0.01 

18 

0.000000 

1.135942 

0.00 

0.0 

5 

3 

0.03 

0.04  -0.02  -0.14 

0.00 

0.00 

19 

1.103975 

1.104571 

0.05 

1.0 

4 

3 

0.02 

0.01  0.00  0.00 

0.34 

0.13 

20 

1.186820 

1.188036 

0.10 

1.0 

8 

3 

0.08 

0.03  -0.04  0.00 

0.19 

0.23 

21 

1.237501 

1.235048 

-0.20 

1.0 

5 

4 

0.32 

0.38  -0.20  0.05 

0.00 

0.01 

22 

1.204551 

1.201810 

-0.22 

1.0 

3 

3 

0.07 

0.00  -0.04  -0.01 

0.40 

0.01 

23 

1.290307 

1.294200 

-0.10 

1.0 

1 

3 

0.20 

0.33  -0.20  0.01 

0.03 

0.08 

24 

1.317507 

1.310531 

-0.07 

1.0 

7 

3 

0.30 

0.31  -0.20  0.00 

0.03 

0.00 

25 

1.321005 

1.320288 

-0.10 

1.0 

5 

5 

0.40 

0.11  -0.02  -0.07 

0.02 

0.00 

20 

1.329940 

1.320090 

-0.29 

1.0 

0 

3 

0.45 

0.11  -0.15  -0.02 

0.11 

0.00 

27 

1.354530 

K 354440 

-0.01 

1.0 

1 

4 

0.20 

0.07  -0.08  -0.02 

0.11 

0.10 

28 

1.357242 

1.356113 

0.00 

1.0 

7 

4 

0.24 

0.10  -0.15  O.Ol 

0.15 

0.09 

29 

1.410719 

1.410488 

-0.02 

1.0 

0 

4 

0.27 

0.09  -0.09  -0.02 

0.24 

0.00 

30 

1.471757 

1.474150 

0.10 

1.0 

0 

5 

0.27 

0.15  -0.15  0.01 

0.21 

0.01 

31 

1.510419 

1.508057 

-0.10 

1.0 

3 

4 

0.10 

0.02  -0.03  -0.01 

0.03 

0.34 

32 

1.511018 

1.515919 

0.32 

1.0 

.7 

5 

0.35 

0.14  -0.04  -0.00 

0.04 

0.07 

33 

1.512835 

1.510014 

0.21 

1.0 

0 

0 

0.13 

0.04  -0.03  -0.01 

0.27 

0.11 

34 

1.510509 

1.518927 

0.10 

1.0 

1 

5 

0.23 

0.20  -0.09  -0.02 

0.14 

0.05 

35 

1.504802 

1.500320 

-0.29 

1.0 

1 

0 

0.11 

0.00  -0.03  -0.01 

0.27 

0.10 

30 

1.500090 

1.507781 

0.11 

1.0 

5 

0 

0.23 

0.21  0.01  0.02 

0.03 

0.00 

37 

1.580001 

1.585057 

-0.02 

1.0 

8 

4 

0.15 

0.09  -0.00  -0.01 

0.31 

0.01 

38 

1.592957 

1.590255 

-0.17 

1.0 

4 

4 

0.14 

0.00  -0.00  0.00 

0.15 

0.22 

39 

1.599704 

1.000803 

0.07 

1.0 

2 

4 

0.10 

0.10  -0.07  -0.01 

0.19 

0.13 

40 

1.007484 

1.010029 

0.10 

1.0 

3 

5 

0.17 

0.19  *0.10  0.03 

0.22 

0.00 

41 

1.008041 

1.011000 

0.18 

.1.0 

5 

7 

0.17 

0.24  -0.14  0.01 

0.07 

0.10 

42 

1.028330 

1.028521 

0.01 

1.0 

7 

0 

0.09 

0.05  -0.03  0.00 

0.27 

0.12 

43 

1.035094 

1.039071 

0.24 

1.0 

8 

5 

0.20 

0.11  -0.08  -0.01 

0.16 

0.11 

44 

1.072210 

1.074042 

0.15 

1.0 

2 

5 

0.15 

0.11  -0.00  -0.01 

0.30 

0.01 

45 

1.738091 

1.735809 

-0.13 

1.0 

4 

5 

0.13 

0.08  -0.00  0.00 

0.20 

0.07 

40 

1.740325 

1.743251 

0.17 

1.0 

5 

8 

0.28 

0.04  -0.05  -0.02 

0.25 

0.01 

47 

1.775304 

1.775280 

0.00 

1.0 

4 

0 

0.12 

0.07  -0.03  -0.01 

0.27 

0.08 

48 

1.782004 

K 779391 

-0.15 

1.0 

7 

7 

0.20 

0.10  -0.10  -0.01 

0.12 

0.13 

49 

0.000000 

1.818747 

0.00 

0.0 

3 

6 

0.07 

0.02  -0.02  0.00 

0.02 

0.41 

50 

0.000000 

1.822177 

0.00 

0.0 

5 

9 

0.28 

0.00  -0.07  -0.01 

0.07 

0.17 

51 

0.000000 

1.820381 

0.00 

0.0 

1 

7 

0.28 

0.07  -0.08  -0.02 

0.07 

0.16 

orthorhombic  olostic  eenttanis  (ell •c22tc33«c234cl3f el2«e44,c5S,c6£) : 


2.05932 

2.05932 

2.58209 

0.99080 

0.99080  0.03820 

0.07740 

0.07740 

.  0.58097 

dimensions 

0.11430 

0.12371 

0.10992 

mats  error  •  0*00000e^00  % 

rms*  .  0« 1^522  percent 

chisquart  increased  2%  by  %  change  in  x's  of 
‘  -0.25  -0.14  -0.06  -1.70  0.01  -0.01  0.00  0.00  0.00 

chisquare  increased  2%  by  %  change  in  s's  ©f 
0.06  -0.31  -0.28  0.53  0.00  0.02  0.00  0.00  0.00 

chisquare  increased  2^  by  %  change  in  k's  of 
-0.12  0.03  -0.13  0.39  0.05  0.00  0.00  0.00  0.00 

xohmrqOl  ctss  time  82.350  seconds 

cpu"  06.355  i/o*  7.158  mem*  8.837 


a11  done 
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Detection  of  microcrackii^  in  a 
single  crystal  sample  of  EuBe 

\)  Resonances  of  the  pristine  sample 
3)  Resonances  after  a  small  Ofees  has  induced 
mlerocracklng 
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We  have  not  succeeded  in  answering  all  of  our  questions. 
Indeed,  we  sometimes  feel  that  we  have  not  completely 
answered  any  of  them.  The  answers  we  have  found  only 
served  to  raise  a  whole  new  set  of  questions,  in  some 
ways  we  feel  that  we  are  as  confused  as  ever,  but  we  think 
we  are  now  confused  on  a  higher  level,  and  about  more 
important  things. 

-Author  unknown 
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ELASTIC  PROPERTIES  CALCULATED  FROM  THE 
HIGH-TEMPERATURE  DATABASE  FOUND  FROM  RUS 


Orson  L.  Anderson 


Supported  by  ONR 


31 


TRMSPMENCrr 


The  vibrational  part  of  the  free  energy  of  an  ensemble  of  inde¬ 
pendent  insulator  oscillators  obeying  Bose-Einstein  statistics 
is  [Landau  and  Lifshitz,  1958] 

(II)  fn  (1) 

where  there  are  3pN  modal  energies  given  by  uji,(jJ2,uJ3  . . .  SpN; 
N  is  Avogadro’s  number;  p  is  the  number  of  atoms  in  the  for¬ 
mula  unit;  h  is  Planck’s  constant;  and  k  is  the  Boltzmann 
constant. 

An  important  approximation  to  the  above  equation  is  pro¬ 
vided  by  the  assumption  that  all  frequencies  depend  on  F, 
but  not  explicitly  on  T.  This  is  called  the  quasiharmonic 
approximation,  defined  by 


and 


where  i  =  1, 2, 3 . . .  SpN. 
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Assume  that  the  acoustic  normal  modes,  /j,  measured  in 
RUS,  follow  the  quasiharmonic  approximation;  that  is 


dfi 

dT 


V 


(4) 


Expressing  fi  as  /j(P,T),  we  have  from  calculus 


By  means  of  a  calculus  equation  involving  P,  V,  and  T,  we 
have 

(6) 

where  a,  the  volume  coefficient  of  thermal  expansivity,  is 
{llV){dVldT)p^  and  Kt-,  the  isothermal  bulk  modulus,  is 
-V  {dP/dV)p. 


Therefore,  within  the  quasiharmonic  approximation,  the  pres¬ 
sure  derivatives  of  frequencies  are  easily  approximated  from 
the  temperature  derivatives,  being 

^ _ 

\dPjp  aKT\dTjp’ 
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Then,  assuming  each  fi  is  linear  in  P,  we  can  find 

MP)t  =  fiiP  =  0,  T)  +  (8) 

The  linear  pressure  dependencies  of  optic  phonon  frequency, 
including  Raman  and  infrared  vibration,  are  found  in  a-  and 
/5-olivine  up  to  at  least  14  GPa  [Chopelas,  1990, 1991;  Hofmeis- 
ter  et  a/.,  1989;  Cynn  and  Hofmeister,  1994]. 

Before  we  can  confidently  use  the  above  equations  to  find 
elastic  constants  versus  P,  we  must  be  sure  of  3  conditions: 
(1)  the  thermal  pressure  (or  aKx)  must  be  independent  of 
V;  (2)  {dfi/dT)y  =  0  must  be  fulfilled;  and  (3)  the  mode 
frequencies  must  be  linear  in  P.  In  principle,  aKx  could  vary 
with  T  and  P  (or  F),  implying  that  the  {dfi/dP)j^  obtained 
from  —  (1/aKT)  {dfi/dT)Y  is  appropriate  for  only  one  co¬ 
ordinate  of  P,  r  space.  Thus,  the  assumption  that  aKT  is 
independent  of  volume  is  necessary  for  condition  3  as  well. 
We  shall  examine  solids  for  which  conditions  1  and  2  are  ful¬ 
filled,  or  one  or  the  other  condition  is  fulfilled.  The  test  of 
the  first  condition  is  whether  aKT  is  independent  of  V. 

I 

I 
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Is  aKx  Independent  of  V? 

When  Pth  is  independent  of  V",  aKx  is  independent  of  V,  as 
is  now  shown.  The  general  form  of  the  equation  of  state  is 

P{V,T)  =  PiiV)  +  Pth{V,T),  (9) 

where  the  first  term  on  the  right,  Pi,  refers  to  the  isothermal 
equation  of  state  (EoS),  and  the  second  term  on  the  right, 
Pth^  refers  to  the  thermal  pressure.  Since  Pi  is  independent 
of  T,  (dP/dT)y  =  {dPTHldT)y^  and  therefore 

Pth  =  J  ocKt  dT  at  constant  V.  (10) 

Therefore,  aKr  must  vary  with  V  in  the  same  way  Pth  varies 
with  V.  The  thermodynamic  identity  relating  aKT  and  V  is 
given  by  Anderson  [1995] 

where  6t  is  the  Anderson-Griineisen  parameter,  a  dimension¬ 
less  thermoelastic  parameter  defined  as— (l/aiTT)  (dKT/dT)p. 

Is  There  Anharmonicity  in  Cy  cd  High  T  ? 

The  literature  has  several  examples,  e.g.,  for  NaCl  and  MgO, 
{dCvldT)p  =  0  at  high  T. 
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Table  1.  MgO:  Physical  Properties  at  High  Temperature 


T 

P 

a 

Ks 

G 

Cp 

7 

Cv 

Kt 

100 

3.602 

0.63 

165.7 

132.0 

0.194 

1.59 

0.194 

165.6 

200 

3.597 

2.24 

164.6 

130.3 

0.662 

1.55 

0.658 

163.5 

300 

3.585 

3.12 

163.9 

131.8 

0.928 

1.54 

0.915 

161.6 

400 

3.573 

3.57 

162.3 

129.4 

1.061 

1.53 

1.048 

158.9 

500 

3.559 

3.84 

160.7 

126.9 

1.130 

1.53 

1.098 

156.1 

600 

3.545 

4.02 

158.9 

124.4 

1.173 

1.54 

1.131 

153.2 

700 

3.531 

4.14 

157.1 

121.8 

1.204 

1.53 

1.153 

150.4 

800 

3.516 

4.26 

155.1 

119.2 

1.227 

1.53 

1.166 

147.4 

900 

3.501 

4.38 

153.1 

116.7 

1.246 

1.54 

1.175 

144.3 

1000 

3.486 

4.47 

151.1 

114.1 

1.262 

1.54 

1.181 

141.4 

1100 

3.470 

4.56 

148.9 

111.5 

1.276 

1.53 

1.185 

138.3 

1200 

3.454 

4.65 

146.7 

109.0 

1.289 

1.53 

1.188 

135.1 

1300 

3.438 

4.71 

144.4 

106.4 

1.301 

1.52 

1.190 

132.1 

1400 

3.422 

4.80 

142.0 

103.8 

1.312 

1.52 

1.191 

128.1 

1500 

3.405 

4.89 

139.7 

101.3 

1.323 

1.52 

1.191 

125.7 

1600 

3.388 

4.98 

137.3 

99.0 

1.334 

1.51 

1.191 

122.5 

1700 

3.371 

5.04 

134.9 

96.7 

1.346 

1.50 

1.193 

119.6 

1800 

3.354 

5.13 

132.7 

94.5 

1.358 

1.50 

1.193 

116.6 
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To  successfully  predict  Cij(P)  from  measured  {dCij/dT)p, 
we  need  three  conditions  satisfied: 


d  (aKr) 
dV 


2. 


=  0 


3. 


Cij{P)  =  Cij{P  =  0)-^ 


To  test  these  three  conditions: 


1. 


Is  Pth  independent  of  V  or  is  6t  =  K', 


2. 

3. 


(Dulong  and  Petit  limit)? 


Is 


fd^CjA 
\  dP^  ) 


(beyond  our  resolution)? 


(16) 

(17) 
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Quasiharmouic  theory  is  based  on  the  assumptions, 


Hierefoie,  we  assume  that  for  the  normal  modes  f,- 

#  «/  I  ) 


measured  in  RUS, 


(18) 
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Sodium  chloride  (NaCl), 
at  ambient  conditions. 


dimension  :  5.947x4.529x3.463  mm3 

density  :  2.160  g/cm3 

Cll  :  496.7  (3.7)  kbar 

Cl2  :  131.8  (4.4)  kbar 

C44  •  127.8  (0.2)  kbar 

aKr  =  (117x10-6  /K)*(242  kbar)  =  0.028557  kbar/K 


mode 

Au-1  (torsion) : 
Blu-1  (llexlure): 
Blg-1  (shear): 
Ag-l  (dilation): 


294  K  571  K 

(Hz)  (Hz) 

180600  176300 
255100  229100 
297300  286900 
370600  320200 


(Of/aT) 

(Or/ap) 

(Hz/K) 

(Hz/kbar) 

-15.5 

543.6 

-93.9 

3286.9 

-37.5 

1314.7 

-181.9 

6371.4 
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d  to.  q.K't 

d  In  V 


=  8t-K' 
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d  aKr  __ 
d  in  V  ^ 
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0-2  0.4  0.6  0.8  1.0 

Pressure  (GPa) 
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d  ccKt 

di-D.V 


Temperature  (K) 
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Temperature  (K) 


46 


TRANSPARENCY  17 


500  1000  1500 

Temperature  (K) 
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600  1000  1600  2000 

Temperature  (K) 
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-  Elastic  Moduli,  (6Pa) 


- a:  Kg  ,  Spetzler  (1970) 

- ’  calculated  result 

- -  c:  Cg  ,  calculated  result 

- d:  Cg ,  Spetzler  (1970) 


I  I  I  I* ■! '  r  I  I  I  I  I  I  I  ■  r"'i "  1  »  1  I  <  f' 

0.5  1.0  1.6  2.0 

Pressure  (GPa) 


Pressure  (GPa) 
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Table  5.  Calculated  RPR  Frequencies  for  Forsterite  at 
High  Pressure 

Mode  300  K  (BfildT)  (BfilBP)  10  kbar 
Hz  Hz/kbar  Hz/kbar  Hz 


An--1 

Ae-2 

B2g-1 

B3g-1 

B2u-1 

Blg-1 

Blu-1 

B3u-1 

B2u^2 

Blu-2 

Blg-2 

B2g-2 

B3g-2 

B3u-2 

B3g-3 

Ag^ 

Au-3 

Blg-3 

B3u-3 

Ag-5 

Blu-3 

B2g-3 

B2u-3 

B1u-4 

B2u-4 

Ag-6 

Blu~5 

B2u-5 

B3u-*5 

B2u“6 


647430 

771500 

959190 

965260 

973720 

1001740 

1004250 

1018320 

1061100 

1120120 

1194370 

1211820 

1219490 

1225260 

1302690 

1331730 

1368250 

1398060 

1418260 

1432110 

1459700 

1468400 

1494980 

1564840 

1575170 

1636520 

1668710 

1688220 

1701340 

1707460 


-55.3 
-76.2 
-78.7 
-93.9 
-63.9 
-95.2 
-61.7 
-80.8 
-78.3 
-88.4 
-91.2 
-93.0 
-86.0 
-85.9 
-105.9 
-96.2 
-123.1 
-121.8 
-104.7 
-  77.0 
-104.7 
-131.5 
-102.1 
-130.0 
-135.0 
-98.2 
-110.1 
-118.4 
-108.4 
-147.7 


1598.3 

2202.3 

2274.6 

2713.9 

1845.4 
2750.0 

1783.2 

2335.3 

2261.6 

2554.9 

2634.4 

2686.4 

2485.6 

2842.7 

3060.7 

2780.3 

3557.8 
3520.2 

3024.6 

2225.4 
3026.0 

3799.1 

2949.4 

3757.2 

3900.3 

2836.7 
3180.6 
3422.0 

3132.9 

4267.3 


663412 

793523 

981936 

992399 

992174 

1029240 

1022082 

1041672 

1083715 

1145669 

1220713 

1238684 

1244345 

1250086 

1333296 

1359533 

1403828 

1433262 

1448505 

1454364 

1489960 

1506391 

1524474 

1602412 

1614172 

1664887 

1700516 

1722439 

1732669 

1750133 
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Table  6.  Pressure  Dependence  of  the  Elastic  Constants  of  Forsterite 


f  ac,f  \ 

3P  )t 

PbIm* 

Saperpoaition 

PalM> 

SaperpontioB 

Bzilbnin’ 

Pnise^Echo* 

Overiap 

Bnllonin^ 

This  Study* 

This  Study 

C\\ 

8.47 

8.52 

9.83 

7X8 

5.45 

10X0 

3X2 

Cn 

6.55 

SJ3 

7J» 

5X9 

5.17 

8.11 

3X8 

Cu 

6.57 

6.21 

7.55 

5X9 

4.72 

8.68 

3.70 

2.12 

2.12 

L72 

1X0 

1.67 

4.09 

2X3 

C«i 

1.66 

1.65 

1.50 

1.40 

1.41 

4.15 

2X9 

a« 

ZZ7 

2.32 

2.20 

2X9 

L76 

4.84 

3X7 

4.11 

3J13 

... 

3.17 

X45 

1.73 

0.24 

C^\ 

4.84 

4.23 

... 

3.40 

3X4 

2X1 

1.46 

Cx2 

4.67 

4.30 

... 

3X6 

3.66 

3.19 

1.77 

^Kumaxawa  and  ^4 nderfon  [1969]. 
^  Graham  and  Barreh  [1969]. 
et  aL  [1982]. 

^  l^nedo  and  Aforroiba  [1992]. 
^Duffj;  et  oL  [2995]. 

^{dfJdT)^  «0. 

^(af/.3T)^  ^0. 
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Summary 

We  have  found  that  for  a  solid  (such  as  Nad)  in  which  aKx 
is  independent  of  V  and  Cy  is  quasiharmonic,  the  inversion  of 
the  measured  {dfi/dT)p  yields  reliable  values  of  \dfi/ dP)j> 
modes.  This  provides  an  accurate  calculation  of  Ks{P)  and 
shear  moduli  up  to  1.5  GPa. 


We  have  found  that  for  a  solid  (such  as  MgO),  where  Cy  is 
quasiharmonic  up  to  very  high  T  but  aKx  is  not  independent 
of  V  at  ambient  conditions,  the  inversion  leading  to  Ks{P) 
yields  good  results,  but  G(P)  departs  from  experiment. 


We  have  found  that  for  a  solid  (such  as  Mg2Si04),  where 
aKx  is  independent  of  V  but  Cy  is  not  quasiharmonic,  the 
inversion  leading  to  Ks{P)  and  ^(P)  yields  good  results,  but 
G{P)  departs  from  experiment. 
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RUS  Measurements  of  Nickel 
Alloys  at  High  Temperature 

J.D.  Games 
Don  Isaak 
Orson  Anderson 


Institute  for  Geophysics  and  Planetary  Physics 

UCLA 
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Overview  of  Topics 

•  Discussion  of  sources  of  funding  and 
interest  for  these  RUS  projects 

•  Materials  studied 

•  Resonant  mode  identification  issues 

•  RUS  spectra  and  resultant  moduli 
determined  as  a  function  of  temperature 

•  Results/conclusions 
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RUS  Projects  involving  nickel 
alloys  at  UCLA 

•  Project  in  collaboration  with  Concurrent 
Technologies  Corporation  (CTC),  an 
independent  nonprofit  organization,  serves 
as  a  national  resource  to  the  industrial  base 

•  Project  in  collaboration  with  Prof.  Alan 
Ardell  of  the  materials  science  dept., 

UCLA 
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(Ni-13.36%A) 
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CMSX4  Mode  ID  Using 
Cy  from  Nickel 

'  I  '  '  '  '  ' - • - 1 - 1 - 1 - 1 - 1—1 — I — I — I — r 

Cii:  250.8  GPa  (Alers  et  al.) 

Cij:  160.0 

C44:  123.5 

I  II  I  I  II  I  II  I  II  I  I  II  III  II  I  III  II  III! 
Top:  Spectrum  Calculated  from  Above  Cy 

I  III  I  II  III!  mil  III  1 1  III  nil  I II  1 1 1 

Bottom:  Spectrum  Observed  from  CMSX4 


' — • — ' — ' — ' — « — I — I — « — I — I — L—J _ I _ ■  ■  I  «  • 

300  400  500  600 

Frequency  (kHz) 


700 
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CMSX4  (Result  of  First  Mode  Assignments) 

“V  '“I  j  I  I  I  I  I  I  I  I  I  [  I  1—1.111  ,  I  II  I  I 

Cii:  223.07  GPa 
C12;  128.58 

C44:  126.68 

i  I  III  II  III!  II  II 1 1  I II 1 1  nil  I  III!  I  III  II I 

Top:  Spectrum  Calculated  from  Above  Cy 

1 1  III  I  II  III!  nil  u  III  1 1  iiiiiiiiii  1 1 II II 

Bottom:  Spectrum  Observed  from  CMSX4 

-< - « - 1 - 1 - « - 1 - 1 - 1 - 1 - 1 - 1—4 _ I _ I _ I _ ■  ■  I  <  ■  ■  ■ 

300  400  500  600  700 

Frequency  (kHz) 
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CMSX4  Mode  ID  (Final  Calculated  Spectrum) 

'  I  '  '  '  •  I — I — I 1 1 — I 1 — r— I 1 — I — I — I 1 1 

Cii.-  256.5  GPa 

Cij.'  166.0 

C44:  131.0 

II  III  II  111  II  mill  I  III  1 1  iiiiiiiiii  II II  III 

Top:  Spectrum  Calculated  from  Above  Cy 

1 1  in  I  II  I  III  INI  II  III  I  I  IIIIIIIIII  I  I  II  II 

Bottom:  Spectrum  Observed  from  CMSX4 

- ' - • - ' - > - « - 1 - 1 - 1 - 1 - 1 - 1 - 1 _ I _ ■  ■  ■  t  .  ■  ■  ■ 

300  400  500  600  700 

Frequency  (kHz) 
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CMSX6  (Comparison  with  Nickel) 

i  I  I - 1 - - — 

Cii:  250.8  GPa  (Alers  et  al.) 

C12:  150.0 
C44:  123.5 

mm  II  II  II  nil  I  I  I  II  i  I  II  I  III  II 

Top:  Spectrum  Calculated  from  Above  Cjj 
M  II  I  II  I II  nil  II  I  II  III  I  III 

Bottom:  Spectrum  Observed  from  CMSX6 


300  400  500  600  700 

Frequency  (kHz) 
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CMSX6  Mode  ID  (Final  Calculated  Spectrum) 


300  400  500  600  700 

Frequency  (kHz) 
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Temperature  (Before  Heating) 


(sijun  /liBjjiqjB)  BpnjjidLuv 
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Frequency  (Hz) 


Frequency  (Hz) 
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Amplitude  (arbitrary  units) 
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Frequency  Shift  with  Temp,  for  CMSX4 
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(sijun  /Uejiiqje)  epnmdujv 
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Frequency  (Hz) 


Frequency  (Hz) 
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Amplitude  (arbitrary  units) 


Frequency  Shift  with  Temp,  for  CMSX6 


Deterioration  of  Acoustic  Response,  CMSX4 
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(s)|un  Ajej)!qje)  epn)!|diuv 
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Frequency  (Hz) 


C44  (GPa) 
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C44  for  CTC  Nickel  Alloys 


Temperature  (C) 
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3  versus  Temperature 
Spec2:  Ni-13.36%AI 
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(BdO)  ®>l 
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T(K) 
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Derivatives  of  Elastic  Moduli  of 
Nickel  Alloys  Studied 


CWISX-4 

ClVISX-6 

SPEC2 

Q. 

o 

-0.043 

-0.045 

-0.038 

0.984 

0.973 

0,999 

dCi2/dT 

-0.010 

-0.013 

-0.010 

0.863 

0.882 

0.974 

dC44/dT 

-0.037 

-0.037 

-0.032 

0.995 

0.997 

0.997 

dGs/dT 

-0.016 

-0.016 

-0.014 

r2 

0.986 

0.986 

0;998 

dKs/dT 

-0.021 

-0.024 

-0.019 

R' 

0.961 

0  948 

0.948 
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Result/conslusions 

•  Temperature  derivatives  determined 

•  RUS  has  proven  useful  in  the 
characterization  of  materials  whose  initial 
Cjj's  were  not  well  known. 

•  RUS  studies  are  valued  in  the 
characterization  of  new  metal  alloys 
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Los  Alamos 

Los  Alamos  National  Laboratory 
Los  Alamos.  New  Mexico  87545 


Response  of  the  Elastic  Moduli  of  Lao.83Sro.i7Mn03 
to  Magnetic  fields  near  the  CMR  Transition 

Temperature 


T.W.  Darling.  A.  Migliori,  J  J.  Neumeier,  E.G.  Moshopoulou, 

J,D.  Thompson,  J.Sarrao’ 

Los  Alamos  National  Laboratory,  Los  Alamos,  NM  87545 
*NHMFL,  Florida  State  University,  Tallahassee,  FL) 


•  Clj  couple  to  free  energy  by  second  strain  derivative  - 
sensitive  to  symmetry  and  coupling 

•  RUS  measures  the  “complete”  elastic  tensor  (as 
complete  as  the  model) 

•  Strong  electron-spin-lattice  effects  in  CMR-  Jahn-Teller 
etc 


1  RUS  primer 

2  Sample  -  preparation  and  analysis 

3  RUS  data 

4  Discussion 
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9,^  -  ^  ''&c"  =  ^0-^5^ 

Asa^/hu  0^  oJ^  ZZoK 
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Lcu^/  pha>^C  Or-Unorho^b/c 
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Orthorhombic 


Figure  1  Shown  are  the  parent  FCC-Iike  unit  cell  and  the  rhombohedral  and  orthorhombic  derivative 
structures  and  their  orientation  with  respect  to  cubic  axes.  The  arrows  on  the  left  cube  indicate  the  direction 
of  the  distonion  that  develops  to  take  the  cubic  structure  to  rhombohedral.  The  two  sets  of  arrows  on  the 
right  cube  show  the  distortion  direaions  that  can  take  a  rhombohedral  structure  and  deform  it  to  monoclinic 
or  to  onhorhombic.  Mn  ions  are  at  the  comers,  the  center  of  the  faces,  the  center  of  the  edges  and  the  body 
center  of  the  large  cube. 
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Figure  4  Shown  are  the  Mn  ions,  a  few  of  the  oxygen 
octahedra  and  the  simple  cubic  parem  unit  ceil. 
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Sample  preparation  and  characterization 

•  Optical-float-zone  method  used 

•  Starting  materials  La203,  SrCOs  and  MnOa  mixed,  ground 
and  reacted  at  high  temperature  then  Isostaticaily  pressed 
into  rods 

•  The  crystai  was  grown  in  air  at  a  rate  of  6  mm/hour 

•  Powder  x-ray  diffraction  detected  no  second  phase 
Refined  cell  parameters  (Hex)  were:  a  =  5.545(1)  A  and  c 
=  13.382(2)  A.  The  x-ray  density  was  6.51427  g  cm‘^ 

•  A  single  grain  was  cut  and  Laue  diffraction  was  used  to 
orient  the  single  crystal  to  the  almost  cubic  lattice  planes. 
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.Aleasurtdi  •'^IHZ) 

.A:n«.(MHZ) 

%error 

k- 

i 

clnA'clnc.. 

2cln/7clnc 

■,j  2cliiAcInCj 

0,835100 

0.832982 

-0.25 

4 

1 

0.06 

-0.04 

0.98 

0.849000 

0.849199 

0.02 

6 

1 

2.77 

-1.88 

0.11 

0.912100 

0.911084 

-0.11 

7 

1 

2.75 

-1.86 

0.11 

1.025500 

1  028140 

0.26 

5 

1 

3.18 

-2.18 

0.00 

1.178600 

1.177794 

-0.07 

5 

2 

3.15 

-2.15 

0.01 

1.184900 

1.182407 

-0.21 

1 

1 

2.91 

-1.99 

0.07 

1.242000 

1  247010 

0.40 

4 

2 

0.30 

-0.21 

0.91 

0.000000 

1.265354 

0.00 

6 

2 

2.28 

-1.47 

0.19 

1.268200 

1.268338 

O.Ol 

3 

1 

0.62 

-0.42 

0.80 

1.278100 

1.280845 

0.21 

5 

3 

3.18 

-2.18 

0.00 

1.299200 

1.297507 

-0.13 

2 

1 

1.69 

-1.15 

0.47 

1.323600 

1.317075 

-0.49 

8 

1 

1.41 

-0.96 

0.55 

1.491900 

1.493765 

0.13 

2 

2 

0.77 

-0.49 

0.72 

1.541500 

1.542163 

0.04 

1 

2 

2.71 

-1.86 

0.15 

1.561300 

1.562077 

0.05 

7 

2 

2.47 

-1.63 

0.16 

1.586500 

1.588603 

0.13 

8 

2 

0.17 

-0.09 

0.93 

RMS  error  =0.2159% 

c„(  lO'Mynes/cm-)=  2.2006±2.28% 
t’,,(10*^dynes/cni')=  1.5142±3.35% 
c«(IO‘’d>Ties/cm‘)=  0.6683x0. 12  % 
B(10'-dynes/'cm')=  1.743  (computed  from  the  moduli) 


T=290K 

mass  =  0.0179  gm  (p=  6.364  g/cm’),  dimensions:  0. 17201  cm  x  0. 13944  cm  x  0. 1 1727  cm 

Table  i .  Shown  are  the  measured  and  fitted  resonance  frequencies,  the  mode  type  (k),  the  mode  harmonic 
(i),  and  the  logarithmic  derivatives  of  frequency  with  modulus.  These  derivatives  describe  accurately  the 
fraaional  dependence  of  each  mode  frequency  on  each  elastic  modulus.  Also  shown  are  the  fitting  errors, 
the  elastic  moduli  (B  is  the  computed  bulk  modulus)  and  their  (larger  than  the  RMS)  errors  determined 
from  the  tit.  and  the  sample  size  and  mass. 
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Modulus  temperature  variation  at  different  B 


INTRO  TRANSPARENCY 


Magnetic  Effects  and  RUS 


J.L.  Sarrao 
LANL 


Collaborators: 

A.  Migliori 
T.W.  Darling 
D.  Mandrus 
Z.  Fisk 


1.  Magneto-elastic  effects 

n.  Field-dependent  effects 

ni.  Details  and  problems  particular  to  RUS 

IV.  Prospects  for  the  future 
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Resonant  frequency  of  single-crystal  fragment 

Undetermined  mix  of  C44  and 

Undetermined  orientation  of  field  with  crystal  axes 
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Variation  of  Neel  temperature  and  softening  with  field  and  with 
C44/Cg6  ratio  (undetermined)  on  several  samples  with  different  mixes 

Tn  essentially  Independent  of  field  direction,  suggesting  free 
energy  couples  as 
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Effects  above  the  Neel  transition  in  applied  field 
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The  above  is  indicative  of  a  weak  phase  transition  of  unknown 
origin 


This  effect  disappears  above  about  70K 
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FIG.  3.  Phase  diagrams  obtained  from  velocity  measure¬ 
ments  for  three  different  orientations  of  the  field:  HJ.c  Md 
HDe.  The  data  points  for  the  phase  boundaries  are  denoted  as 
i,  •.  fftii  and  ^  T*,*-  ^  our  resets  for 

Te*(H)  and  from  velocity  measurements  (A)  in  am- 

rarison  with  the  heat-capadty  data  of  Refs.  9  and  16  (shown 
u  4.),  ji^ote  that  we  have  normalized  the  temperature  scales  in 

the  insets  for  the  different  samples. 
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Flux  lok-Wriee 

Sc,»  Cm  Oy 


Figure  1:  Change  of  sound  velocity  as  a  function 
of  tenq)erature  in  BsTTesla  for  diree  doping  levels 
(B|c,  k|c  and  ki.u). 


Figure  2:  Phase  diagram  of  the  effective  field 
B^^B'COs^  vs.  T*  with  B  |c  (□),  by  rotating  the 
field  BarTTesla  in  the  (IlO)-plane  (0,0  and  a) 
or  in  the  (110)-plane  (■).  Open  symbols:  28MHz, 
full  symbols:  8.9MHz, 

?.  iernmens  e-VoA.,  3  I^M  Iifc3 
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ELECTRON-PHONON  INTERACTION 


295 


Fig.  40.  Acoustic  quantum  oscillations  in  LaB*  for  the  (c,|  -  c,2)/2  and  modes  for  different 
geometries  (Suzuki  et  al.  1985a). 
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U,  -  I  Cj„— — 


i.iJ(J 


dV, 


The  “standard”  problem 
Computational  time  =1  (1) 


le«so<-  T»»o«UcV^  Wo<- 
J  I 

2  2  XT’  Iwi/ 


ijjkj 


mjn 


kjjn 


dx„  dx^ 


Electric  field  effects 

Computational  time  =  16  (7.6) 

solved  by  Ohno  [Phys.  Chem.  Minerals  17  (1990)  371] 
0.5%  effects  on  frequencies 


'•=/ 


^Xj  C/Xi  „jf 


UJfJ 


kjjn 


dV, 


Magnetic  field  effects 
Computational  time  =  64  (1 9) 

<i|j  #  Oji  due  to  torques 

(so  81  elastic  constants  —  generally  ignored  in  pulse-echo) 
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Experimental  Field  Alignment 
Pulse-echo 


RUS 

Rotatable  probe  head 


RST 


INTRO  TRANSPARENCY 


High  Temperature  RUS  and  the 
Elastic  Constants  of  Steel 


Tim  Darling 
Philip  Armstrong 
Albert  Migliori 

Raj  Vaidya 
Carolynn  Scherer 

Los  Alamos  National  Laboratory 


♦  RUS  Application 

♦  “Real  Time”  following  of  elastic  constants  during 
heat  treatment 

♦  Improved  elastic  constants  for  51XX  grade  steels 
as  a  function  of  temperature 


0  Comparison  of  two  resonance  methods  -  RUS  and 
rod  resonance 

0  Comparison  of  single  mode  vs  full  fit 
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HMnw  Risufck  UbantoriMt  Itthtohan  tUal  Cwp. 


Coolini  Timi,  See 


Diffusive  transformations  (slow,  high  temperature 
activated)  compete  with  non-diffusive  (Martensitic, 
fast,  low  temperature  activated)  transformations  to 
produce  a  complex  combination  of  phases, 
microstructures  and  strains. 
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Samples:  5180,  5140,  5120  steel  (assume  isotropic) 
RUS:  RPR  4.8  X  5.0  X  5.2  mm 
Rod:  Va  diameter  x  4”  long 

Procedure: 

1.  heat  treat  and  quench  to  get  “100%”  phases 

2.  measure  resonance  frequencies  with  temperature 

♦  RUS:  Lithium  Niobate/Alumina  metal  bonded 

transducers 

T  <  500C 

Measure  and  fit  40-50  lines  ->  Cn  and  C44 
Sensitive  to  symmetry  and  homogeneity 

♦  Rod:  Centre  node  mount,  capacitive/magnetic 

transducers 

T<  1000c 

Measure  fundamental  longitudinal  and 
shear  modes  E  and  G  (Young’s  and 
shear  modulus) 
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MAGNET  AND 
PICKUP  COIL  Fo: 
SHEAR  MODE 
(DETECTOR) 
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Alumina  ceramic 

Stainless  steel 
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Table  4.  RR  measurements  of  Youni 
Temperature  (C)  I  Length/L^ 


1.0000 


1.0009 


1.0014 


1.0020 


1.0035 


1.0053 


1.0053 


1.0054 


1.0072 


1.0085 


1.0091 


1.0000 


1.0001 


1.0013 


1.0032 


1.0081 


1.0117 


1.0117 


1.0110 


1.0102 


1.0099 


1.0000 


1.0109 


1.0129 


1.0129 


1.0124 


1.0117 


1.0110 


1.0104 


1.0104 


1.0102 


1.0101 


1.0101 


1.0098 


1.0095 


1.0095 


1.0095 


1.0095 


1.0094 


1.0094 


1.0093 


1.0094 


1.0094 


1.0094 


1.0113 


1.0107 


1.0105 


;*s  modulus  -  5180  Pearlite  i 
I  Longitudind 


laEUHiSitaiQj 


25.581 


25.626 


25.362 


25.277 


24.935 


24.429 


24.423 


24.389 


23.819 


23.226 


23.113 


25.629 


25.633 


25.445 


24.992 


23.561 


20.378 


20.358 


22.108 


22.53 


22.743 


23.166 


25.656 


22.295 


19.398 


19.438 


19.638 


19.913 


20.103 


20.273 


20.283 


20.383 


20.412 


20.383 


20.453 


20.643 


20.803 


20.933 


21.013 


21.193 


21.413 


21.734 


21.777 


21.813 


21.82 


20.68 


21.659 


22.07 


normalized) 
Young's  Modulus 
GPa 


211.3 


211.8 


207.4 


205.9 


200.0 


191.7 


191.6 


212.5 


158.7 


119.9 


120.4 


123.0 


126.5 


129.0 


131.3 


131.4 


132.8 


133.2 


132.8 


133.7 


136.3 


138.4 


140.1 


141.2 


143.7 


146.6 


151.1 


151.7 


152.2 


152.3 


136.5 


149.9 


155.6 
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Sample  msOl  5180  steel  ms  per  TWD  250C 

mass=  0.8919  gm  rfao=  7.643  gm/cc  fraction  of 

Cll  C44 


n 

f-expt 

f-model 

%err 

weight 

k 

i 

df7d(moduli) 

1 

0.273840 

0.272363 

-0.54 

1.0 

4 

1 

0.00 

1.00 

2 

0.303360 

0.302188 

-0.39 

1.0 

4 

2 

0.00 

1.00 

3 

0.367920 

0.369609 

0.46 

1.0 

6 

1 

0.05 

0.95 

4 

0.380640 

0.380712 

0.02 

1.0 

7 

1 

0.06 

0.94 

5 

0.383040 

0.384115 

0.28 

1.0 

3 

2 

0.01 

0.99 

6 

0.393230 

0.393999 

0.20 

1.0 

2 

1 

0.01 

0.99 

7 

0.000000 

0.399060 

0.00 

0.0 

1 

1 

0.04 

0.96 

8 

0.400830 

0.401842 

0.25 

1.0 

8 

1 

O.OI 

0.99 

9 

0.000000 

0.434440 

0.00 

0.0 

5 

1 

0.00 

1.00 

10 

0.444150 

0.443232 

-0.21 

1.0 

6 

2 

0.32 

0.68 

11 

0.000000 

0.446697 

0.00 

0.0 

5 

2 

0.00 

1.00 

12 

0.000000 

0.457719 

0.00 

0.0 

5 

3 

0.00 

1.00 

13 

0.469230 

0.468069 

-0.25 

1.0 

8 

2 

0.15 

0.85 

14 

0.000000 

0.468456 

0.00 

0.0 

7 

2 

0.29 

0.71 

15 

0.479870 

0.479327 

•0.11 

1.0 

2 

2 

0.15 

0.85 

16 

0.000000 

0.479395 

0.00 

0.0 

1 

2 

0.30 

0.70 

17 

0.498000 

0.497998 

0.00 

1.0 

3 

3 

0.15 

0.85 

18 

0.000000 

0.535291 

0.00 

0.0 

8 

3 

0.01 

0.99 

19 

0.544800 

0.542180 

•0.48 

1.0 

5 

4 

0.35 

0.65 

20 

0.000000 

0.552534 

0.00 

0.0 

5 

5 

0.34 

0.66 

21 

0.553120 

0.554548 

0.26 

1.0 

4 

3 

0.05 

0.95 

22 

0.556240 

0.555254 

•0.18 

1.0 

2 

3 

0.01 

0.99 

23 

o.oodooo 

0.566521 

0.00 

0.0 

6 

3 

0.10 

0.90 

24 

0.000000 

0.569224 

0.00 

0.0 

7 

3 

0.09 

0.91 

25 

0.571580 

0.569744 

-0.32 

1.0 

1 

3 

0.08 

0.92 

26 

0.576260 

0.574978 

-0.22 

1.0 

3 

4 

0.01 

0.99 

27 

0.000000 

0.618932 

0.00 

0.0 

7 

4 

0.08 

0.92 

28 

0.000000 

0.619188 

0.00 

0.0 

1 

4 

0.10 

0.90 

29 

0.000000 

0.626533 

0.00 

0.0 

6 

4 

0.05 

0.95 

30 

0.648000 

0.650363 

0.36 

1.0 

1 

5 

0.09 

0.91 

31 

0.658500 

0.661049 

0.39 

1.0 

7 

5 

0.09 

0.91 

32 

0.661200 

0.663693 

0.38 

1.0 

6 

5 

0.10 

0.90 

33 

0.703800 

0.705123 

0.19 

1.0 

5 

6 

0.89 

0.11 

34 

0.706200 

0.706173 

0.00 

1.0 

6 

6 

0.06 

0.94 

35 

0.000000 

0.714476 

0.00 

0.0 

3 

5 

0.11 

0.89 

36 

0.000000 

0.721312 

0.00 

0.0 

2 

4 

0.09 

0.91 

37 

0.730800 

0.730772 

0.00 

1.0 

7 

6 

0.06 

0.94 

38 

0.000000 

0.733186 

0.00 

0.0 

8 

4 

0.07 

0.93 

39 

0.000000 

0.736673 

0.00 

0.0 

4 

4 

0.00 

1.00 

40 

0.738900 

0.739970 

0.14 

1.0 

3 

6 

0.09 

0.91 

41 

0.000000 

0.742701 

0.00 

0.0 

2 

5 

0.10 

0.90 

42 

0.000000 

0.743813 

0.00 

0.0 

5 

7 

0.10 

0.90 

43 

0.747900 

0.747665 

-0.03 

1.0 

8 

5 

0.10 

0.90 

44 

0.751200 

0.749611 

-0.21 

1.0 

1 

6 

0.06 

0.94 

cll  cl2  c44 

2.5530  1.0972  0.7279 

dl  d2  d3 
0.51150  0.49000  0.46560 

lo^  5  rms  error®  0.2799  % 

chistpiare  increased  2%  by  the  following  %  changes  in  independent  parameters 

Cll  :  0.40% 

C44  :  0.09% 
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Fr«qu«ney  variation  of  lino  33  (94%  C  f,)  eemparad  to 
■  *  C,,  fittad 


O 


o 

-q 

0 


Comparison  of  fuU  RUS  data  fit  and  single  pure  mode  fit  for  Cn 


Froqudticy  variation  of  lino  1  (100%  C^)  comparad  to 


Comparison  of  full  RUS  data  fit  and  single  pure  mode  fit  for  Cm 
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(Bdo)  9 


113 


Poisson's  ratio 
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Variation  of  Poisson's  ratio  (v)  from  RUS  and  RR  data 
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Mode  frequency  kHz 
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5180  steel  sample  A9  -  resonance  frequencies 
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51 80  steel  sample  A9  T-30C,  First  pass 

TWO  4.14-9$ 
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Amplitude 
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5180  steel  sample  A9  T=30C,  Second  pass 


400000  410000  420000 

Frequency  Hz 
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+  Fulfilled  the  “engineering”  aspect 

E(T) ,  G(T)  for  a  number  of  phases 


+  Happy  agreement  between  RUS  and  rod  method 


?  Did  not  look  at  attenuation  -  this  may  be  the  most 
important  piece  of  ultrasound  information  when 
watching  transitions  in  microstructured  materials 


?  Connection  between  observed  splitting  and 
distortion  of  the  RPR 
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RESONANT  ULTRASOUND  STUDY  OF  TEXTURE 
IN  POLYCRYSTALLINE  MATERIALS* 

Keir  Foster  and  Robert  G.  Leisure 
Department  of  Physics 
Colorado  State  University 
Fort  Collins,  CO  80523 

and 

George  Alers 
NIST 

Boulder,  CO  80303 


*  Work  supported  by  the  Colorado  Advanced  Materials  Institute 
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OUTLINE 


I.  Introduction 

II.  Effect  of  texture  on  elastic  constants  of  polycrystalline 
materials. 

III.  RUS  measurements  on  textured  copper:  elastic  constants  and 
orientation  distribution  function  coefficients. 

IV.  Ultrasonic  pole  plots 

V.  Conclusions 
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INTRODUCTION 

•  Texture 

•  Nonrandom  orientation  of  crystallites  (grains)  in 
polyciystalline  aggregates. 

•  Arises  in  fabrication  processes  (rolling,  drawing),  welding, 
etc. 

•  Has  strong  effects  on  formability,  and  on  nondestructive 
testing. 

•  Methods  of  Investigation 

•  X-ray  diffraction 

•  Neutron  diffraction 

•  Ultrasound 

•  Objectives  of  Present  Study 

Use  Resonant  Ultrasound  Spectroscopy  to  investigate  the 
weak  anisotropic  elasticity  due  to  texture. 
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TEXTURE  AND  ELASTIC  CONSTANTS 


For  plates,  =  “rolling  direction,  X2  =  “transverse  direction”,  and 
“normal  direction. 

•  Euler  Angles 


Sample  Coordinates 


rystallite  Coordinates 
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•  Transformation  between  the  crystallite  and  sample  coordinate 
systems 

Vector 


Elastic  constant  tensor 


3  3  3  3 


;7=l^=lr=l^l 


^  jCl 

pi  qj  rk  si 


The  transformation  matrix  is 


a- 


cos  y/  cos  ^cos^zJ-sinj^/’sin^^ 
-  cos  cos  ^  sin  -  sin  cos ^ 
cos  sin  ^ 


sin  y/  cos  9  cos  ^ + cos  sin  ^ 
-  sin  cos  ^sin ^ + cos  cos ^ 
sin^/sin^ 


-sin^cos^"' 
sin^sin^ 
cos^  i 
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•  Elastic  constants  of  crystallite  in  sample  coordinate  system 

Transform  the  crystallite  elastic  constants  to  the  sample  coordinate 
system. 


Results  for  a  crystallite  with  cubic  symmetry. 


C  ' 

C  '  = 

^22 

C  '  = 


JJ 


C  ' 

^44 

C  '■ 

^55 

C  ' 

^66 

C  '  = 

2j 

C  ' 

^13 

C  ' 

^12 


C^i-2Cr^ 

■■q,-2Cr^ 

Cji-2Cr3 

Q4+C/-3 

C44+Crg 
Cj,+Cr^ 
Cj2  +  Cr3 


'•1=  ^1^21+^1^31+44 

^2  ~4‘^22 +44 +44 
^3  ~ 44 +44 +44 

r  ^(p-  cP  ’{‘Cp  Cp  -^-Cp  Cp 
^4  2^1 3  ~  z2^23  ~  j  2^3  3 

^5  ^23^21  '^^33^31 

^6  ^11^12  "^^21^22  "^^31^32 


'V  “  ^2^^93  ^32^' 

^8  ~  ^n^f3  “^^21^23  '^‘^31^33 
+*^2 1'^22 


2 

:>o 


=  Cj2+C/-g 


C=q.-Q-2C 


✓44 


Note:  Other  C,y  are  ignored  because  orthorhombic  symmetry  for  the 
sample  will  be  assumed  later. 
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•  Orientation  distribution  function  ir(5,i|;,(|))  (5  =  €Os(8)) 

=  probability  of  finding  crystallites  with  orientations 
between  (ij;,0,(j))  and  (i|;  +  ,0  +  <i0,  (J)  +  d^). 


It  will  be  convenient*  to  expand 


l=0m=-l  n--l 


-imij/^-in(p  ^ 


where  the  are  generalized,  associated  Legendre  functions.  The 
are  the  Orientation  Distribution  Coefficients  (ODC*s)  and  are 
determined  experimentally. 


*R.-J.  Roe,  J.  Appl.  Phys.  36, 2069  (1966) 
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I 

•  The  <  C.>  of  the  sample,  in  the  Voigt  approximation,  will  be  the 

average  of  the  C^.  over  the  orientation  distribution  function.  These 
will  be  of  the  form 

<  C[  j  >=<  Cj  j  >  -2  <  C  ><  Tj  > 


etc. 


To  calculate,  expand  the  r, 


Calculate  the  averages 


Ik 


2k 


<r.>=  JJJ 
0-10 


00  rn=+ln=+l 

=4;z-2y  y  yR‘  w. 

Imn  Imn 

/=0w=— /  «=-/ 
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•  Finally! 

<C^j  >=A+2/z+  +^'W 

<q,-  >=A+2/H|M^^j(^+2|0)^+:^i^ 

<q.'>=Ai+2^4jM^^J^ 

<q^'  >=^-|26i2<^j^_ 

<q,' 

<q.  >=^-Ji6^<^j(F^+ 

<q.' 

<q,' 

The  theoretical  values  of  the  Lame  constants  X  and  \i,  and  the  anisotropy 
constant  <C>,  depend  on  the  averaging  methods  used  (Voigt,  Reuss, 
Hill,  etc.)  for  the  texture-free  material.  The  sample  is  assumed  to  have 
orthorhombic  symmetry. 
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Resdnant  Ultrasound  Spectroscopy 


I.  Sample-Transducer 
Arrangement 


1  -2  mm 


II.  Spectrometer  _ 

Transmitter 


-  Receiver 


III.  Typical  Modes 


Signal  (voll^) 
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Olin  Type  Cl  10  polyciystalline  copper 
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EXPERIMENTAL  RESULTS  FOR  OLIN  TYPE  Cl  10  COPPER 

(Copper  +  0.02 1  wt  %  Oxygen) 


Measurements  on  three,  independently  prepared,  copper  samples  at 
room  temperature. 


Elastic 

Constant 

Sample  1 

Sample  2 

Sample  3 

Mean 

%  Error 

Cn(GPa) 

199.7 

198.4 

201.9 

200.0 

0.51 

C22(GPa) 

199.9 

199.1 

201.0 

200.0 

0.28 

C'33(GPa) 

195.7 

197.6 

197.4 

196.9 

0.31 

C'44(GPa) 

44.06 

45.24 

45.47 

44.96 

1.02 

C-„(GPa) 

44.04 

45.04 

45.66 

44.91 

1.05 

C'«(GPa) 

44.59 

45.79 

45.41 

45.26 

iO.78 

C'23(GPa) 

109.8 

107.5 

112.0 

109.8 

1.18 

C,3(GPa) 

109.2 

106.1 

109.0 

108.1 

0.93 

111.7 

109.7 

110.0 

0.82 

A  least  squares  fitting  of  the  sample  1  results  to  the  theoretical  expressons, 
using  <C>  =  -100.25  GPa,  gives. 


W400 

-2.36  X  10-* 

W420 

+4.56  X  10*^ 

W440 

-5.11x10-"' 

X  (GPa) 

109.1 

\x  (GPa) 

44.37 

Computed  values  using  single  crystal 
C„  for  copper. _ 


Voigt 

Average 

Reuss 

Average 

X  (GPa) 

101.5 

111.0 

p  (GPa) 

54.6 

40.0 
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Using  the  values  of  the  JT  X  and  we  compute  the  and  compote  to  the 
results  for  sample  1. 


Experimental 

Computed 

C'„(GPa) 

199.7 

C%2(GPa) 

199.8 

C'33(GPa) 

195.7 

C'44(GPa) 

44.06 

44.26 

C'ssCGPa) 

44.04 

44.17 

C'66(GPa) 

44.59 

44.34 

C'23(GPa) 

109.8 

109.1 

C',3(GPa) 

109.2 

108.9 

C',2(GPa) 

108.6 

109.1 

RMS  Error  =  0.7% 
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•  Reciprocal  lattice  vector  orientation  distr&iitioB  fiinctioii 

=  cos(Xi)) 


Consider  a  reciprocal  lattice  vector  tj 


Crystallite  Coordinates 


Sample  Coordinates 


q.{g,rj)dgdri~  Probability  to  find  vector  t|  between  (x,t|)  and  (%+  dx  ,j] 
+  (ivi) 


It  can  be  shown* 

?,(^'.7)=^+5(Hp®,)[i^(O»'400  +  2P4'(Ocos(2;7)ff;20+2iy(Ocos(477)fF^ 


Ej.  =cos(0.) 


where  F”i(z)  are  the  normalized,  associated  Legendre  fimctions. 

•  Ultrasonic  pole  plots 

Contour  plots  of  ^i(C»ri)  for  a  particular  reciprocal  lattice  vectors  t* 


*M.  Hirao,  K.  Aoki,  and  H.  Fukuoka,  J.  Acoust  Soc.  Am.  81, 1434  (1987). 
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[100]  pole  figure  for 


Olin  type  Cl  10  copper 
(Sample  1) 


Unanneaied 


1=  Q 


Rolling 

Direction 
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[100]  pole  figures  for 
Olin  type  Cl  10  copper 
(Sample  1) 


unannealed 


p 


Annealed  at  250  C 
for  1  hr. 


p 


Annealed  at  300  C 
for  1  hr. 


p 


Roiling 

Direction 
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[110]  pole  figures  for 
Olin  type  Cl  10  copper 
(Sample  1) 

Unannealed 


Annealed  at  250  C 
for  1  hr. 


p 


Annealed  at  300  C 
for  1  hr. 


Roiling 

Direction 
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CONCLUSIONS 
PRESENT  RESULTS 

RUS  elastic  constant  measurements  on  Olin  type  Cl  10  copper 
described  by  orthorhombic  elastic  constant  matrix,  but  indicate 
very  small  anisotropy  -  in  disagreement  with  sound  velocity 
results. 

Qj'  appears  to  deviate  from  the  theoretical  expressions. 


OUTLOOK 

Measurements  on  materials  with  stronger  anisotropy  needed,  with 
a  careful  comparison  to  ultrasonic  velocity  measurements. 

RUS  has  potential  to  characterize  texture. 


INTRO  TRANSPARENCY 


COMPARISON  OF  RUS  MEASUREMENTS 

with 

Conventional  Methods 
C.M.  Fortunko  &  H.  Ledbetter 


IMIST 

Materials  Science  and  Engineering  Laboratory 
Division  853.05 
325  Broadway 
Boulder,  Co  80303 

ph.  (303)  497-3062  e-mail  fortunko@boulder.nlst.gov 

fax  (303)  497-5030 
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COLLABORATORS 


•  M.C.  RENKEN 

•  K.W.  HOLLMAN 

•  S.  KIM 

•  N.  SIZOVA 
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TECHNICAL  APPROACH: 

•  ESTABLISH  “TRUE”  LONGITUDINAL-WAVE  VELOCITY  IN 
A  HOMOGENOUS  ISOTROPIC  SPECIMEN  USING  A 
NUMBER  OF  SHORT-PULSE  EXPERIMENTAL 
GEOMETRIES  AND  SIGNAL-ANALYSIS  METHODS. 

•  ESTABLISH  LONGITUDINAL  VELOCITY  USING  RUS  ON 
THE  SAME  SPECIMEN. 

•  CALCULATE  THE  ERROR 
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WE  NOW  USE  RESONANT  ULTRASONIC 
SPECTROSCOPY  FOR  ROUTINE  DETERMINATION  OF 
ELASTIC-STIFFNESS  COEFFICIENTS.  IT  IS  THEREFORE 
OF  INTEREST  TO  ESTABLISH  THE  ACCURACY  OF  THE 
METHOD. 
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COMPONENTS  OF  ACCURACY 


ACCURACY 


BIAS 


PRECISION 
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Fused-Silica  Specimen 


^ - 25.4  mm - > 


p=  2.202  g/cm^ 
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SHORT-PULSE  EXPERIMENTAL  GEOMETRIES 

(PULSE-ECHO) 


•  WATER-IMMERSION 

•  BUFFER-ROD 

•  DIRECT-CONTACT 
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Some  Physical  Configurations 
(Pulse-Echo) 

(1)  Liquid  immersion 


i 


. ' 


XDucer 


(2)  Buffer-Rod 


XDucer 


; : ! :  j  M  •  • . 

Sample 

1 _ 

i=ii 

Sample 


Direct  Contact 

XDucer 


Sample 
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SHORT-PULSE  SIGNAL  ANALYSIS  METHODS 

•  PULSE-ECHO  SUPERPOSITION/VISUAL  OVERLAP 

•  AMPLITUDE-SPECTRUM/PI-POINT 

•  PHASE-SPECTRUM/DECONVOLUTION 

•  PULSE-ECHO  CROSS-CORRELATION/AUTOMATED 

•  TIME-FREQUENCY/(UNDER  DEVELOPMENT) 
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RESOLUTION  OF  SHORT-PULSE  METHODS 


At 


1)  L  is  known  to  0.008% 

2)  At  is  affected  by  Bandwidth,  sampling  rate  and  signal  to 
noise  (0.04%) 
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"As  Measured"  Velocities  (No  Corrections) 


[Fused  SlUca 

Oveiall  Chart 

(mmAjs) 

No  Corrections 

1(25.634  mm) 

Sional  Processina 

\mmg^ 

1 

Method 

IQSSSiSEH 

Amplitude  Spectnim 

Phase-Spectrum 

Pulse-Echo 

Tim  A/Pront  lonru 

Pl-Point 

Deconvolution 

Immersion  Tank 

— 

Pulse>Echo  Confia. 

5.9309354 

5.927888889 

5.935156707 

5.929244444 

Dry  Contact  PZT 

Puise-Echo  Confia. 

5.827608667 

5.9144 

5.870609784 

Oiv  Contact  Qua/tz 

Pulse-Echo  Confiq. 

5.928316667 

5.638533333 

5.609666667 

5.846210537 

Dry  Buffer  Rod 

1 

Pulse-Echo  Confiq. 

5.93035 

5.9221 

5.933796296 
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Physical  Sources  of  Error 
(Short-Pulse  Methods) 


Transducer  has  finite  thickness 


i - - - 

— 

r 

(1) 


■»l  h 


Temperature  Differences 


Puise-Shape  /  Bandwidth  Effects 
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Velocities  (with  Corrections) 


Fused  Silica  1  Overall  Chart 

Ledbetter's  iWIth  Corrections 

(25.634  mtnt  ISipnat  Pfocfissifia 

iMsmod 

Puise-Echo 

Amdituda  Soectnmi 

Phase-SoBCtruni 

Pulse-Echo 

Tima/FreouBncv 

Suoemosftion 

Pf-Point 

Cross-Correiation 

Analysis 

! _ 

Puise-Echo  Confiq. 

5.9309354 

5.927888889 

5.9247602 

5.935156707 

5.929244444 

Drv  Contact  P2T 

1 

Pulse-Echo  ContiQ. 

S.fl42104123l  5.87992341 

—  1 1 1 1  1 1  Il'TT 

5.665337594 

Orv  Contact  Quaitz 

Pulse-Echo  Confici 

5.93256504 

5.914261234 

5.685182592 

1  Sa328B2349 

5.93035 

5.9221 1 

5.9337962961  1 

•  Without  Diffraction  Correction  (estimated  to  be  less  than  0.1%) 

Coefficient  of  thermal  expansion  is  5.5e-7.  Temperature  Range  =  2.7  C 
Temperature  correction  is  less  than  0.01  % 
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"TRUE"  Value  for  Velocity  =  5931 .0  m/s 
(as  calculated  from  below) 


Average  Velocitv  without  outliers 

_ with  corrections _ 

_ (mttvus)  _ 

_ 5.93101 6935 


Standard  Deviation  without  outliers 


0.00336829 


With  Error  of  +/-  0.11%  using  a  95%  confidence  intervai 
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Short-Pulse  Velocity  Comparison 


Residual 

T  Velocity 

Residual 

% 

% 

Immersion  Tank 

IHHI 

Dry  Contact  PZT 

Pulse-Echo  Confiq. 

iBHUM 

3EEHEQ 

Pulse-Echo  Confiq. 

5.842104 

1.499116 

Pulse-Echo  Suoerpositjon 

Pulse-Echo  Superoosition 

Immersion  Tank 

Dry  Contact  PZT 

Pulse-Echo  Confiq. 

0.05274 

Pulse-Echo  Confiq. 

Ampiituda  Spectrum 

Ampiituda  Spectrum 

Pl-Point 

Immersion  Tank 

Dry  Contact  PZT 

Pulse-Echo  Confiq. 

5.92476 

0.105492 

Pulse-Echo  ConfiQ. 

5.920362 

0.179639 

Phase  Soectrum 

Phase  Soectrum 

Immersion  Tank 

Dry  Contact  PZT 

Pulse-Echo  Confiq. 

5.935157 

0.069799 

Pulse-Echo  Confiq. 

0.770177 

Cross-Correlation 

F 

Cross-Correlation 

Immersion  Tank 

Dry  Contact  Quartz 

Pulse-Echo  ConfiQ. 

Pulse-Echo  Confio; 

EPJcW.-f.i-il 

Time/Frequency 

Pulse-Echo  Superposition 

Al  Buffer  Rod  (Dry.  PZn 

Dry  Contact  Quartz 

Pulse-Echo  Confiq. 

5.93035 

0.011245 

Pulse-Echo  Confiq. 

5.914261 

0.28251 

Pulse-Echo  Superposition 

Amplitude  Spectrum 

1 

Al  Buffer  Hod  (Dry,  PZTI 

Pulse-Echo  Confiq. 

5.9221 

0.150344  1 

Dry  Contact  Quartz 

Phase  Soectrum 

1 

Pulse-Echo  Confiq. 

5.685183 

0.772786 

1 

^hase  Spectrum 

Al  Buffer  Rod  (Dry,  PZTl 

Pulse-Echo  Confiq. 

0.046861  1 

Dry  Contact  Quartz 

Cross-Correlation 

T 

1 

>ulse-Echo  Confiq. 

_ _ L 

i 

Dross-Correlation 
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RUS  Experimental  Geometry 


Pinducers 


Specimen 
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RUS  SIGNAL  ANALYSIS  WE  USE  METHODS 
DEVELOPED  BY  P.  HEYLIGER  (COLORADO  STATE 
UNIVERSITY)  CALCULATES  ELASTIC-STIFFNESS 
MODULI  FROM  OBSERVED  RESONANT-FREQUENCY 
SPECTRUM.  THE  METHOD  CAN  HANDLE  A 

VARIETY  OF  SHAPES  AND  CRYSTAL  SYMMETRIES, 
INCLUDING  TRIGONAL. 
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ERRORS  IN  RUS  METHOD 

VALUE  OF  C  IS  IMPACTED  BY  TRANSDUCER 
CLAMPING  AND  SPECIMEN  GEOMETRY  (0.08%) 

p  IS  KNOWN  TO  0.004% 
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COMPARISON  OF  VELOCITIES 

SHORT-PULSE:  5931.2  m/s  ±  6.7  m/s 
which  is  0.1 1  %  with  a  95%  Confidence  Intervai 

RUS:  5939.8  m/s 

ERROR  is  0.15% 


157 


TRANSPARENCY  22 


ANISOTROPIC  MATERIALS 


COMPARE  RUS  ELASTIC  STIFFNESS 
CONSTANTS  FOR  HIGH-QUALITY  SILICON 
SPHERE  WITH  THOSE  OBTAINED  BY 

MCSKIMIN  (0.1%)  4  DECADES  AGO. 
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Silicon's  Voigt 


present 

literature 

ratio 

C-11 

165.61  ±0.03 

165.78 

0.9990 

^12 

63.84±0.03 

63.91 

0.9989 

C44 

79.55±0.01 

79.61 

0.9992 

c 

50.89 

50.94 

0.9990 

B 

97.76 

97.87 

0.9989 
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NEED  FOR  ACCURATE  VELOCITY  MEASUREMENTS 


Internal  Strain 

0.1% 

0.01% 

Pressure  Effects 

0.1% 

0.01% 

Texture 

1  % 

0.1% 

Phase  Mixtures 

0.1% 

Hardness 

0.1% 

Defects 

0.1% 

0.1% 
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INTRO  TRANSPARENCY 


Sampled  CW  investigation  of 
vortex  viscosity  and  phase 
transitions  in  YBa2Cu307.5  single 

crystals. 


D.  Dasgupta,  J.  R.  Feller,  C  Hucho*, 
Bimal  K.  Sarnia  and  M.  Levy 


Department  of  Physics 
University  of  Wisconsin-Milwaukee 

Work  sponsored  by  Office  of  the  Naval  Research 


Present  address  :  FB  Physik,  Germany 
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OUTLINE 


•  Sampled  CW  technique 

•  Attenuation  vs  Magnetic  Field 

•  At  low  magnetic  fields  ( <  0.5T ) 

•  At  high  fields  upto  1.6T 


Attenuation  vs  Temperature  in  different 
constant  fields 


TRANSPARENCY  2 


What  we  do .. 

Ultrasonic  studies  on  high  temperature  superconducting 
single  crystals. 

Sampled  CW  technique. 


Setup  schematic  : 


RF  Contact 


Fine  threads 
ofGE  Varnish 


a ,  b 


!s'.- 


YBCO  crystal 


Fine  gold  leads 
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How  we  do  it ... 


Excitation  frequency  :  5  MHz 
Longitudinal  velocity  ~  5000  m/s 

Typical  crystal  thickness  :  50  pm 

long  wavelength  limit 
bulk  measurement 


X  ~  1  mm 


Current  densities  : 

Sampled  CW  method  ~10'^A/cm^ 

Transport  measurement  ~  0.5  A/cm^ 


Attenuation  units 
Sensitivity 


1  dB/jj.s  or  2dB/cm 
lO^dB/ps 


L65 


^t^enuAftoo  C 


TRANSPARENCY  5 


General 


^pplie^l  AtoancKe  PieM 

B  C  TeeU  ) 


66 


ATTENUATION(dB/us) 


transparency  6 


Affpinufttipn  vs.  •motgng'h'c  j^tet 

(  H  c>>»b  plane.  ;  T  g  8t_iJS5  K) 


^  lattice 

oC  £:  B*'* 

yoYttfx  -  voyIsx  'in'tv>9Ai'oo 
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log  a(dB/us) 
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Log  V5«  lops  ot 


f 


er£ 


ni  const  X 


log  B(T) 

increases  wtHi  increasing  T 
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B(T) 


1.7 


Temperaturc(K) 
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Summary 


•  Vortex  pinsite  interaction  is  dominant  at  iow  fields 

(B.pp<0.3T) 

•  Intervortex  interaction  is  dominant  at  higher  fields 

(B.pp>0.3T) 

•  Low  frequency  ultrasound  measurement  shows 

transition  in  vortex  system  : 


•  At  iow  fields  flux  lattice  is  soft 

•  At  high  fields  flux  line  lattice  is  rigid 
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Sampled  Continuous  Wave  Investigation 
of  the  Vortex  State  in  YBa2Cu307.5  Single  Crystals: 
Measurement  Considerations 


University  of  Wisconsin-Milwaukee 
J.  R.  Feller,  D.  Dasgupta,  M.  Levy 


Work  supported  by  Office  of  Naval  Research 
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Possible  Reasons  for  High  Sensitivity 

A  O 

•  Low  current  density  (J  <  10  A/cm  ) 

J  ~  1  AJcvr?  for  transport  (dc  resistivity)  measurements 

•  Long  acoustic  wavelength  (A-  ~  1  mm)  compared  with 

sample  thickness  (~  50  pm) 

•  Low  amplitude 

•  Weak  acoustic  coupling 
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Sampled  Continuous  Wave  Measurement  of  YBCO 

Single  Crystals 


Magnetic  Field 

- ► 


RF 


■^a,b 


YBCO 


Ground 


5  MHz  Quartz 
Transducer 


GE  Varnish 


Resolution  ~  5  x  10’^  dB/|is 
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Weak  Coupling  and  Sensitivity 

Define  a  coupling  parameter:  0  <  a  <  1 

Effective  bond  area:  =  aA^^p,,  0  <  ^  A^p,, 

Assume 

a  =a, +aa,  a„:  measured  attenuation 

m  t  s  m 

a  t :  transducer  attenuation  coefficient 
a^:  sample  attenuation  coefficient 

Q.=^7— 7 - ^ 

2(a,  +  aa3) 

Define  the  sensitivity  as 

s  =^  =  -“— ?— 

5a,  2  (a,  +  aa,)^ 

|S„|  is  maximized  by  varying  a: 

^  =  0=>  at=aa3  a^=2cL^  or 
odi  ^ 

For  Qt  =  10,000  and  Qs  =  300,  the  maximum  in  sensitivity  occurs 
when 

a  =  0.03  and  Qm  =  5000. 
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Weakly-coupled  system 


force-free 

boundary 


transducer 


bond 

sample 

force-free 

boundary 


A,  Area  of  transducer 

Aj  Area  of  sample 

A^  =  aA,  Effective  area  of  bond 

D  S 

a  Coupling  parameter 

0<a<  1 


Acoustic  impedance:  Z  =  Apco/k 
Wave  number:  k  =  ca/v  -  jot/2 


The  bond: 

Effective  area: 
Effective  density: 
Effective  sound  speed: 
Effective  impedance 


Ab  =  aA3 

apb 

avb 

(aA3)(apb)(av^)  =  a% 


L,(a)  oc  a2  Cb(a)  oc  a"* 
As  a  0,  0,  C(,  00 
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Resistance 


Reactance 


- r 

- 1 - 1 - 

1 

_ i_ 

r 

_ L _ \ - 

r~ 

_ 1 _ 

3.5  4  4.5  5  5.5  6 

Frequency  (MHz) 
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Amplitude  Amplitude 
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Sampled  Continuous  Wave  Measurements 


Frequency  (MHz) 


Frequency  (MHz) 
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Conclusions 


•  Sampled  continuous  wave  technique  demonstrated  to  be 
highly  sensitive  to  vortex  state  transitions  in  YBCO 
platelets 

•  Factors  contributing  to  high  sensitivity:  Low  current 
density,  long  acoustic  wavelength,  low  amplitude,  weak 
acoustic  coupling  between  transducer  and  sample 

.  Sensitivity  )  found  to  exhibit  a 

maximum  in  the  weakly-coupled  regime 

•  Agreement  with  preliminary  frequency-swept 
measurements 
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Thin  Film  Characterization  Using 
Resonant  Ultrasound  Spectroscopy 

Jason  White  and  J.  D.  Maynard,  Department  of  Physics  and  the  Acoustics 
Program,  The  Pennsylvania  State  University,  University  Park,  PA  16802 
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.  Survey  of  thin  film  characterization  techniques 
.  Theory  of  RUS  on  thin  films 
.  Numerical  calculations 
.  Experimental  methods 
.  Conclusion 
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Thin  Film  Characterization 
Techniques 


I.  Direct  measurement  of  deflection  in  response  to  an 
applied  force 

a.  Bulge  test 

b.  uniaxial  tension  test 

c.  nanoindentatlon 

II.  Ultrasound  methods 

a.  Free-standing  films 

b.  Film  and  substrate 

1 .  SAW  techniques 

2.  Resonance  method  ^ 

ill.  Brillouin  Scattering 
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Basis  function: 


i 


fz) 

\b) 


m 


fiT 

\c) 


Matrix  Eigenvalue  Equation: 

(D^Ea  =  r  a 
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©^Ea  =  r  a 

t  t 


1  1  1 


E  =5..  fffpO  ^dV 

W  ^7  J  J  ^ 


-1-1-1 
1  1  1 


=  5. 


u 


1  1  1 


-1-1-1 


r.= 


y/r  a®,  a®^ 

-a-b-c 


dXj^  dxi 


a  o  c 

=4  JJf 


“  *  ^  ^<D  5<1)  “r^f  f 


-a-b-c 


dxj^  dx, 


-a 


-b(c-ac) 
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Si  specimen  500x400xl00jLim 
A1  film  on  500x400  face  from  10  to  0.01  jim  thick 


no  film 
(Mhz) 

10|im  film 
(%A/) 

1  |im  film 
(%A/) 

0.1  fim  film 
(%A/) 

0.01  fim  film 
(%A/) 

2.5791 

9.571 

1.031 

.1037 

.0104 

2.7515 

7.329 

.8026 

.0851 

.0081 

4.2864 

6.823 

.7433 

.0740 

.0074 

5.3460 

8.584 

.9233 

.0928 

.0093 

5.8465 

8.106 

.8757 

.0881 

.0088 

6.0633 

3.459 

.3371 

.0335 

.0033 

6.9605 

6.489 

.7060 

.0706 

.0071 

7.1318 

3.346 

.3353 

.0335 

.0033 

8.0462 

4.291 

.4204 

.0420 

.0042 

8.3537 

3.388 

.3301 

.0328 

.0033 

8.5105 

3.367 

.3224 

.0320 

.0032 

8.7763 

7.568 

.8196 

.0824 

.0082 

9.7565 

6.496 

.7075 

.0709 

.0071 

11.0223 

6.742 

.7281 

.0731 

.0073 
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Experimental  Setup 


INTRO  TRANSPARENCY 


System  Identification 
Alternatives  to  Swept 
Sine  and  Curve  Fitting 


Ralph  Muehleisen 
ASEE  Fellow 
Department  of  Physics 
Naval  Postgraduate  School 
muehleis@physics.nps.navy.mil 
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Purpose 


•  Provide  some  alternative 
approaches  to  Swept  Sines  and 
Curve  fitting  for  RUS 

•  Why?  ....  Why  not? 

While  the  traditional  approach 
works,  there  might  be  another 
approach  which  is  better  for 
some  circumstances 
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Traditional  Approach; 
Advantages/Disadvantages 

•  Advantages 

-  It  works  and  works  well 

-  Not  difficult  to  implement 

•  Disadvantages 

-  Not  particularly  fast 

-  Requires  some  specialized 
equipment  (e.g.  lock-in  amp) 

-  Not  “real  time”  or  even  close 

-  Not  especially  suitable  for 
production  line  testing 
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Step  1 :  Replace  Swept 
Sine  with  Sampled  Noise 

•  Wideband  noise  excites  sample 

•  Instead  of  H(a)),  ARM  A  System  ID 
methods  are  used  to  find  impulse 
response  h[n]  and  hence  the  discrete 
time  poles  and  zeros  of  the  system 

•  Resonance  freq.  and  damping  are 
extracted  directly  from  h[n]  and  the 
standard  inversion  process  is  used  to 
find  Cjji^ 
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•  Advantages: 


-  A  standard  high  speed  data  acq. 
card  can  be  used 

-  No  lock-in  amplifier  is  required 

-  Should  be  faster  than  Swept  Sine 
but  avoids  FFT  problems 


-  Technique  has  worked  well  for 
structural  modal  analysis  problems, 
and  and  in  other  systems  with  well 
separated  poles. 

-  Literature  is  replete  with  papers  on 
ARMA  system  and  spectral 
estimation. 
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Step  2:  Replace  Curve  Fit 
with  Adaptive  Filter  system 


•  Instead  of  finding  resonances  and 
doing  a  non-linear  curve  fit,  use  an 
adaptive  filter  technique  for  online 
identification  of  the  elastic  constants 
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•  Operation 

-  The  filter  is  adjusted  until  the  error 
goal  is  reached. 

-  There  are  many  algorithms  for 
adjusting  (LMS,  RLS,  etc.) 
most  of  which  are  forms  of  a  non¬ 
linear  least  squares  curve  fit. 

•  Advantages: 

-  True  online  measures  of  are 
obtainable  (great  for  production 
testing) 

-  Adaptive  algorithms  are  usually 
faster  than  their  standard  least 
squares  counterparts  because  they 
use  efficient  routines  for  updating 
the  error  gradient  estimate 
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The  Step  Beyond: 
Directly  Determine  C^/^/ 

•  Reformulate  the  problem  entirely  -  do 
not  solve  the  variational  based 
eigenvalue  problem 

•  Instead  convert  the  bounded  elastic 
wave  equation  to  a  discrete  time 
equivalent  (maybe  a  state  space 
system?) 

-  Why?  This  might  yield  an  exact 
ARMA  filter  with  coefficients 
proportional  to  the  Cp 

-  Determination  of  the  ARMA  model 
would  then  determine  the  Cjj^ .  This 
might  be  done  in  real  time. 
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Summary 

•  Use  ARMA  system  ID  methods  to 
replace  the  Swept  Sine 

-  Faster  determination  of  resonances 

-  Requires  only  a  fast  data  acq.  board 

•  Use  Adaptive  Filters  to  replace 
the  Curve  Fitting 

-  Faster  determination  of  Cp 

-  True  online  operation 

•  Use  alternate  form  of  problem 

-  Cjjid  might  be  found  directly  and  in 
real  time. 

•  Techniques  useful  for  real  time 
production  line  testing  of  samples 
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Resonant  Ultrasound  Spectroscopy 
Measurements  for  Spheres  Under 
Pressure  With  Different  Gases 


D.G.  Isaak,  J.D.  Games,  O.L.  Anderson 

UCLA 

H.  Oda 

Okayama  University,  Japan 
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Goals 

•  Use  RUS  to  measure  3rd  order  elastic  moduli. 

•  Pressure  dependence  of  elastic  moduli  from  RUS. 

•  Effects  of  gas  loading  on  frequency  shifts? 

•  Test  dCij/dP  from  pulse-echo  experiments. 


I 
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“The  effect  of  the  bond  has  been  one  of  the 
most  abominable  problems  plaguing  high 
precision  ultrasonic  measurements.  There 
has  been  a  long  history  in  studying  bond 
effects  . . .  Conclusive  results  on  the  nature 
and  the  effect  of  the  bond  are  still  lacking.” 


Spetzler,  Chen,  Whitehead,  and  Getting,  Pure  and  Applied 
Geophysics,  141,  341-377,  1993 
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Report  on  Fused  Silica 

•  Two  spherical  specimens  (r=2.4154, 1.945  mm) 

•  Three  gases  (helium,  nitrogen,  argon) 

•  Two  specimen-holding  configurations 

•  Two  torsional  modes,  T02  and  T(,3 
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Dependence  of  (dG  /dP  )j  on  Molecular  Mass  of  Gas 
-  Helium,  Nitrogen,  Argon  - 
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^(de/oe)  *0  apnimBeifli 


201 


Molecular  Mass  (g) 


Frequency  (Hz) 
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Amplitude  (arb.  units) 
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Large  Sphere 
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(s:^|un  'qje)  epn^nduiv 


203 


Frequency  (Hz) 


Frequency  (Hz) 
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Amplitude  (arb.  units) 
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Small  Sphere 


Small  Sphere 
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(s)|un  -qje)  epn^nduiv 
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Frequency  (Hz) 


Frequency  (Hz) 


TRANSPARENCY  9 


Amplitude  (arb.  units) 
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(s\[un  -qje)  apn;i|duiv 
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Frequency  (Hz) 
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Pressure  Dependence  of  Frequencies 
-  Helium  Gas  - 


Pressure  Dependence  of  Frequencies 
-  Nitrogen  Gas  - 
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Pressure  (bar) 
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f(P)-f(0)  (Hz) 
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770217  1.23038  49.2 

1190124  1.90116  76.1 


Molecular  Mass  (g) 
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Magnitude  of  (dG/dP).^ 


CO  CO  CO  CO  ^  ^ 


o 
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Dependence  of  (dG  /dP  on  Molecular  Mass  of  Gas 
-  Helium,  Nitrogen,  Argon  - 
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Conclusions  on  Fused  Silica 

•  Small,  measurable  effect  of  pressurizing 
gas  on  dG/dP. 

{dG’l6m  =  6x10'^  g"*,  OT~moIe.  mass  of  gas) 

•  Magnitude  of  dG/dP  is  3.70  (for  m=Q), 
8%  higher  than  for  pulse-echo. 

•  Evidence  of  bond  effects  in  pulse-echo? 
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Gas  Pressure  Dependence  of  Radiation  Impedance 

of  a  Resonating  Sphere 


University  of  Wisconsin-Milwaukee 
J.  R.  Feller,  R.  Sorbello,  D.  Beck,  M.  Levy 

University  of  California,  Los  Angeles 
Institute  of  Geophysics  and  Planetary  Physics 

D.  Isaak 


Work  at  UW-M  supported  by  Office  of  Naval  Research 


214 


TRANSPARENCY  1 


Outline 


Part  1. 

•  Discussion  of  Isaak’s  RUS  data  on  the  silica 
sphere. 

•  Pressure  dependence  of  Q  and  resonant  frequency 
of  the  T02,  SI  1,  SOO  modes. 


Part  II. 

•  Theoretical  examination  of  the  radiation 
impedance. 

•  Dependence  of  resonance  Q  on  pressure  and 
molecular  weight. 

•  Expected  frequency  shift  due  to  mass  loading. 

•  Comparison  with  experimental  results. 
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Fused  Silica  Sphere  (Radius  2.414  mm) 


T 

S 

S 

S 

T 

S 


N 

0 

0 

1 

0 

0 

0 


L 

2 

2 

1 

0 

3 

3 


Freq.  (kHz) 

620.57691 

652.04321 

804.52548 

944.70449 

958.90296 

961.98376 


Measurements  of  D.  Isaak,  UCLA: 


Frequency  (kHz) 
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Fused  Silica  Sphere  (Radius  2.414  mm) 


955  957  959  961  963  965  967  969 

Frequency  (kHz) 
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Amplitude 
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Fused  Silica  Sphere  (Radius  2.414  mm)  T02  Mode 


Degenerate  modes:  n  =  0,  1  =  2,  m  =  0,  ±1,±2 


m  =  0: 

Ue  =0 

=  A  sin 20  coscot 


m  =  ±l: 

Uq  =  A  COS0  cos(a)t  ±  (j)) 

=  Acos20cos(cot±(|)) 

m  =  ±2: 

Uq  =  A  sin0  cos(cot  ±  2(|)) 

=  yAsin20cos(cot±2(|)) 
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Frequency  (kHz) 


TRANSPARENCY  5 


Fused  Silica  Sphere  (Radius  2.414  mm)  T02  Mode 


0  20  40  60  80  100  120 


Pressure  (bar) 
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a 


measured 


-  tto  +aR 


1 


Q  measured 


1 

Qo 


Residual  loss: 


=  constant « 10 


-5 


Radiation  loss: 


_j_  JVp  T02 
Qr'^Ip  SOO,  Sll 


f 

measured 


=  fo  +  Af 


f-  ^/P  T02 

I  -P  SOO,  Sll 
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Amplitude 
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Fused  Silica  Sphere  (Radius  2.414  mm)  SOO  Mode 


n  =  0,  1  =  0,  m  =  0 
=  A  coscot 

Ue  =  =  0 

(Pure  radial  breathing  mode) 
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Fused  Silica  Sphere  (Radius  2.414  mm)  SOO  Mode 


Amplitude 
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Fused  Silica  Sphere  (Radius  2.414  mm)  SI  1  Mode 


Frequency  (kHz) 


Degenerate  modes:  n  =  1,  1=1,  m  =  0,  ±1 
m  =  0: 

=  A  COS0  coscot 
Uq  =  Bsin0cos0t 
=  0 

m  =  ±l: 

Ur  =  Asin0cos(cot±(|)) 

Uq  =  Bcos0  cos(cot  ±  (|)) 

=  Bcos(©t  ±  ((|)  +  f )) 
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Fused  Silica  Sphere  (Radius  2.414  mm)  SI  1  Mode 
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WAVE  EQ: 


THEORY 

92ui  3Tij 
"  dxj 


Ui 


=  f-th  component  of  displacement  u 


=  X  V  u  5ij  +  |J,  ( 


auj  du 


3xj  axj 


)  =  stress  tensor 


SOLUTIONS  (See,  e.g..  Lapwood  &  Usami.  1981 1: 

I.  IRROTATIONAL;  Uj  =  V  $  ei®t 


V^O  +  =  0  h  =  co/cl  =  co/'v/  (A,+2|i)/p 

O=j;(hr)Y;„(ft0) 

II.  EQUIVOLUMINAL  with  0^=0:  Un=Vxr\|t  ei®t 

V2\)/  +  =  0  k  =  co/cj  =  (d/a/  |j,/p 

t|/=j;(kr)Y;^(a0) 


III.  EQUIVOLUMINAL:  Um  =  (1/k)  V  x  V  x  r  x  e>“t 
V^X  +  k2x  =  0  X  =  j/(kr)  {6.  (p ) 


SPHEROIDAL  MODES :  u  =  A  Uj  +  B  Um 
TOROIDAL  MODES;  u=  u,. 
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EIGENFREOUENCIES  &  RADTATTON  IMPEDANCE 

FREE  SURFACE  B  .C. :  T^.  =  T0j.  =  =  0 

Determines  resonance  frequency  cOq  for  unloaded  sphere. 

LOADED  SURFACE  B.C.: 

Gas  exerts  force  F  on  surface  of  sphere  (r  =  a), 


and 

^Gr  =  =  0  (neglecting  viscosity  of  gas). 

RADIATION  IMPEDANCE: 


Z  =  — ^  (evaluated  at  r=a) 

Vgas 


=  ^ -  (evaluated  at  r=a) 

1  CO  Uj. 

Outgoing  radiation: 

p  =  Po  h/2)(kgr) 


3v 

where  kg  =  co/cg,  together  with  eq.  of  motion  -Vp  =  ^ 

°  di 

gives: 

Z  =  -ip^  Cg  h/2)(kg  a)  /  h/2)'  (kg  a) 
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COMPLEX  FREQUENCY  SHIFT: 

BC  at  r=a:  Tj^-  =  -i  co  Uj.  Z 

or  F(co)  =  -i  CO  Z,  where  F(co)  =  Tj^  /Uj.. 


Solution  gives  complex  eigenfrequency  co  =  cOq  +  Aco. 
Small-loading  expansion: 

F(o)o)  +  f'  ((Hq)  Ao)  =  -i  COq  Z(cOo) 

/.  Aco  =  -i  cOq  Z(cOo)  /  f'  (cOq) 

All  information  about  surrounding  gas  is  contained  in  Z(cOq). 
Further,  since  F'  (cOq)  is  real. 


and 


ReAco/coQ=  ImZ(c0o)/F'(c0o) 
Im  Aco/cOq  =  Re  Z(cOo)  /  f'  (cOq) 


The  Quality  factor  Q  =  (Oq/(2  Im  co)  can  then  be  obtained  from 


Q 


=  2  Re  Z(cOq)  /  f'  (cOq) 
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SOO-MODE  (I  =  0): 

Uj.  =  A  jQ((or/cL)  cos  cot  (Uq  =  =  0) 

Tjj  =  X  (duj./dr  +2  Uj./r)  +  2  duj/dr  ;  F(co)  =  T^-  /  Uj- 
Transcendental  equation  as  before:  F(co)  =  -i  co  Z 

where  Z  =  -  i  h/^^(kg  a)  /  hp^  (kg  a) 

-nr  r  1  +  i  (kg  a)->  1 
s  L  1  +  (kg  a)-2  J 


For  reported  expts  ,  kg  a  =  16  for  He  and  49  for  Ar.  So 
dependence  on  molecular  mass  M  of  gas  and  static  pressure  p  of 
gas  can  be  obtained  by  taking  kg  a»l  limit  for  Z: 

Re  Z  =  p^  c^  oc  p  M  ^4-1/2=  p  m1/2 
Af  oc  Im  Z  =  Pg  !  (kg  a)  =  p^  cp  /(co  a)  «=  p  MM“^  =  p 


Theoretical  results  reported  here  for  1/Q  and  Af  are  obtained  using 
full  expressions  for  Z  and  F  (cOq). 
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2  for  SOQ-y^ode^ 


?/f,S 
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1000/Qp  Af(kHz) 
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Fused  Silica  Sphere  (Radius  2. 


Nitrogen 


Argon 


Helium  (Experiment) 
Nitrogen  (Experiment) 
Argon  (Experiment) 
Helium  (Theory) 
Nitrogen  (Theory) 
Argon  (Theory) 


Pressure  (bar) 
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Sll-MQDE  (I  =  l.m  =  0±n: 

u  =  A  Ui  +  B  Uiji 

At  r=a:  Tj^  =  (a A  +  pB)  (0,(|)) 

T0,  =  (yA  +  5B)  (e,(^  )/a0 

T(j)r  =  (yA  +  6B)  im  Y/^  (0,(|)  )/sin  0 
a  =  -(X  +  2|i)h2ji(ha)  -  (4jih/a)j/(ha)  +  (4^/a2)jj(ha) 
p  =  (4^/k)(kj/(ka)/a-ji(ka)/a2) 

Y  =  2^(  h  j  1  *(ha)/a  -  j  j  (ha)/a2) 

5  =  |Lik  j  1  "(ka)  (h  =  co/cl  ;  k  =  co/cj) 

UNLOADED  SPHERE: 

B.C.  at  r=a:  Tj^  =  Tqj-  =  =  0 

/.  (aA  +  pB)  =  0  and  (yA  +  5B)  =  0. 

Consistency  requires  determinant  =  0.  a6  -  Py  =  0  gives  cOq. 
LOADED  SPHERE: 


B.C.  at  r=a:  1)  /  Uj-  =  -i  co  Z 

'^0r  ~  '^(|)r  ”  ^ 

Using  Uj- ,  Tit  and  T0J.  in  these  equations,  get  usual  form: 

F(co)  =  -i  CO  Z 

where 

a5  -  py 


F(co)  = 


h  ji’(ha)/a  -  (2y/ka)  ji(ka)/a2 


Numerical  results  for  Af  and  1/Q  shown  in  figure. 
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Fused  Silica  Sphere  (Radius  2.414  mm)  SI  1  Mode 


V^«: 


Mass  Loading  Af 


..  '"A-..,  Nitrogen 

.  Argon  ^ 


•  Helium  (Experiment) 
A  Nitrogen  (Experiment) 
■  Argon  (Experiment) 

-  Helium  (Theory) 

.  Nitrogen  (Theory) 

.  Argon  (Theory) 


Helium 


A  A  ^  ^  A _ ^ 


Pressure  (bar) 
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TOR  OTP  AL  MODES 

Uj.  =  0  no  coupling  to  a  non- viscous  gas.  1/Q  =  0  and  Af  =  0.) 

Need  to  include  viscous  coupling. 

Onnsider  toroidal  modes  with  m=0: 

Uf  =  ug  =  0  =>  Tn-  =  Ter  =  0 

u,j,  =  j;(kr)  3P;  (e)/90 

T(j,j.  =  p.(k  j/>(kr)  -  j;(kr)/r)  aP;  (e)/ae 

UNLOADED  SPHERE  BC  at  r=a;  T^r  =  0  gives  (Oo. 
LOADED  SPHERE  BC  at  r=a: 

S3.S  (b 

-fig — ^  where  =  viscosity  coeff.  of  gas. 
Assuming  viscous  length  5  «  a: 

Vgas  (|)  °^  e  with  k’  =  and  5  =  V  Pg®/2}J.g 

/.  =  vgas(l, 

Defining  Z  =  'piscous  /  y  ^  the  BC  at  r=a  has  usual  form 

F(co)  =  -ia)Z 

where  F(co)  =  T(|)j.(U(j))/U(|)  =  fi(k  j/^^ka)  -  j;(ka)/a)/j;(ka),  and 

Z  =  *^(l+i) 

Note  Zoc  pin  since  jig  is  independent  of  p. 
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Fused  Silica  Sphere  (Radius  2.414  mm)  T02  Mode 


23A 
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CONCLUSION 


•  Our  radiation  impedance  analysis  gives  results  in  rather  good 

agreement  with  measured  Q  values.  Our  frequency  shifts  Af 
are  much  smaller  than  the  measured  values,  which  implies  that 
the  measured  Af  are  not  due  to  mass  loading  by  the  surrounding 
gas.  This  is  consistent  with  the  data  on  the  T02  mode,  where 
there  is  very  little  dependence  of  Af  on  gas  species.  However, 
the  measured  Af  for  the  SI  1  and  SOO  modes  show  large 

contributions  to  Af  that  depend  on  gas  species,  and  this  is 
perplexing. 

•  All  effects  of  gas  loading  are  contained  in  Z,  from  which  we 
can  obtain 


Af/f=  Im  Z(coo)  /  F' (cOq) 

1/Q  =  2  Re  Z(coo)  /  F'  (cOq) 

where  all  relevant  details  of  the  elastic  medium  for  the  unloaded 
case  are  contained  F  (cOq)* 


•  The  last  two  relations  give  the  very  general  result 

Af/f  Im  Z(cOo) 
l/Q  "  2  Re  Z(cOo) 

which  can  be  applied  to  samples  of  more  general  shape, 
provided  that  the  appropriate  Z  is  calculated  for  such  samples. 


235 


TRANSPARENCY  22 


VJ 


VJ  l'iiv  ^^cc_  ^e.  lo  c^'  “O e  % 


l'\’g.\  Vvt  V\e\Vz.  O  ^  •^KOWT  OL  "1^  Ko*.  He. /g.^  g  d 

'•Ay  S  i  <j  e>  ^  o.v\  j  C  orV^r^ 

<x  v  S  -Vru  c,  -Lu  ^re,  V  i  t  v-a  -Lt  vt  ^ 

lo  cJ  “Oe  S  ^  /Vc-v> 

^ ^  ^ (x.'V  z  «  "TU<.  y3V'ess'*rt.  oo’(“s»A.«.  -I-U  «_ 

f  cP;  =  -z  kyoc  ff<jc;  y^-^Ck) 


ZttR 


w  Lys.  1?  =  +(y-j,f  +(z-zj 

CLyid  “{"lye.  'f  e.  ^  vo. -4>io v\  is>  CYcy 


^ur-Co-ce^  o 


t  fU 


fLi 


f 


V'a.Hel 


j0(]f9( 


.a. 


^.C.,  o~\ 


^onstJcV'  4^^  +a.c.e_  ex'\' 

"Wn  S  ■^<3L  C  <_  C«Ov/\  V  r  '.  \>  vj  Ve  ^  (X  •Vev 

”  ^  ^  c-  ^ 


+< 


A 


•  Vf, 


It  \°^ 

I  Veve,  cuYc..  A\S6 


iiovi  ■1'''' l^U^V  >  0v>  S  "Vo 


'^YO'^  “I'Lc,  d 'HieA'  ■^O-C.-e,  S  d'f'  Va  llcje. 


j-'  e. 

O  — 


{*^j 


ZrR 


y\-U 


W  ll  trc.  I  vy-l'e  ^  kci-”£^  ovy  /‘s  ower  cx// 

•fie-  ?■  ^(Le-S  p  "1"  ovn^  ^4"  Z- 
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Oo  vv^  yo  u  "V  \ovA  cx  \  Sc.U 


e.w\^  • 


C-V\  J 


\J  SG^ 


T  kg 


Hi 


2irR 


JT 

n-  ^ 


(7= 


JX 


e+c . 


■Vo  '\"VcVolV^  "("O  6e-l  f  —  Co  V\  s  i  ‘t  4'<-v^*V  So\u’0«^^S 

o  “V^e.'S  C.  e.^  C^<SL  -^‘owr^. 

I  -C-ic-vi  \  V  ieS  J  IXve.  i  \A-Vecjk  \r<xls  4L«- 


4-. 


cjcc-es  C.e/w\-f  4^  -fKe-  "Z  =  t'z 

kokV«-  loi^OLvi  WA.'C-  Si\AQ  «j  loLVi o.-f 

ticlaJ  :  Z=5^<^X=^_,Z^%'4J-^,«^‘. 


=  94  ^ 


cx 


>3-  > 

C-oVk  si  ole-^ 

-  ^  fK  cj 


>iig 


S-rR 


-  L 


k 


% 


k 


Jj(,  f  t/zv 


kg 


grg 


'U^  Cjtj  z«) 


x-f 


-  /, 


y. 


?'4  ^ 


'-t 


+  y--) 

7.=  -%  ‘t-7--|- 

5o\u'Viov\  «  Covvi^u'Ve.  7"  *  ■§■ 


W  U<2,Y<2,  Vv  '  ^  <>-  ^oSl4\V<> 

S/-.4  l\  =  O  -(•ov'  4Le.  c-a\  dul<3L'(^'ov\  o  ^• 


41,. 
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INTRO  TRANSPARENCY 


Determining  the  Radiation  Impedance 
for  Arbitrarily  Shaped  Surfaces 


J.  D.  Maynard  and  Jason  White 
The  Pennsylvania  State  University 


Technology  Issue 

Calculate  the  sound  radiation  (acoustic  impedance)  for  a  vibrating  surface. 

Applications 

•  Noise  control 

•  Improving  sound  radiation  from  hi-fi  speakers,  musical  instruments,  etc. 

•  Fluid  loading  (Resonant  Ultrasound  Spectroscopy  for  samples  under  pres¬ 
sure) 


Objective 


Develop  a  method  for  calculating  the  acoustic  impedance  for  a  vibrating 
surface  with  an  arbitrary  shape. 
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Acoustic  impedance  for  a  vibrating  surface 


Physical  surface  with  surface  normal  velocity  v(fs)e~^‘^  =  n  ■  v  (fs) 

Wave  Eqn.  V^P-ljc^d^PIdt^  =  0  -)•  Helmholz  Eqn.  V^P+k^P  =  0,  fc  =  uj/c. 

Boundary  value  problem;  v(fs)  determines  unique  solution  P(f) 

Evaluate  P  (f)  at  r  =  fs.  Acoustic  impedance  is  Z  (fs)  =  P  (fs)  /v  (rs) . 

In  RUS,  Im{P(fs)}  =  Im  {Z  (fs)}  v  (fs)  gives  stress  boundary  condition, 
shifts  frequencies. 

l/2Re  {P  (fs)}  V  (fs)  gives  dEldt,  the  acoustic  energy  radiated  to  infinity.  In 
RUS,  quality  factor  Q  =  uE/  (dE/dt)  . 
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The  Green's  Function  for  the  Boundary  value  Problem 
Green’s  function  G(r)  satisfies: 

+  k^G  =  -AtzS  if) 

AND  G{rs)  =  0  or  n-VG{fs)  =  0 
on  the  boundary  Ts. 

Solution:  P(rO  =  J  J  v(fs)G(f-fs)drs.  NOTE:  G(f-fs)^^\ 


The  Green  s  function  may  be  found  if  the  boundary  surface  corresponds  to 
the  level  surface  of  a  separable  coordinate  system. 
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The  Heimholz  eqn.  has  Eleven  separable  coordinate  systems 


Cartesian  coordinates  x,y,z 

^i{k,x+kyy+k,z) 


Cylindrical  coordinates  z 

H„  (k,p) 


Spherical  coordinates  r,6,4) 
hm(kr')e^^‘^P[^  (cose) 
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Cartesian  coordinates  x,y,z 


General  soln. 
Surface  velocity 


P  (2:,  y,  ?)  =  j  J  A  (kx,  ky) 


xdky 


,  ,  1  dp  / 

V  (x,  y)  —  ^  (z  =  0) 


ipck  dz 
1 

47r2 


Fourier  Transform 


Solution  P(x,y)  = 


v(kx,ky)  =  j  j  v(x,y) 


e^(k.x+k,y)dk  dk 


■XU»A/y 
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Arbitarily  Shaped  Sources  not  level  surface  of  a  sep.  coord,  sys. 


Must  find  the  Green’s  function  numerically 


Use  discrete  points 
OR 


Use  isoparametric  boundary  elements 


Green’s  Thm.  P  (r)  =  J  J  (fs‘)  n  ■  VG  (r- fs‘)  -  iv  (fs‘)  G  (f- fs‘)jdrs' 
where  G  (R)  =  with  R  =|  r  —  rs  |  . 

Evaluate  at  r  =  fs^  Surface  Helmhofz  Intergral  Equation  (SHIE) 


P  (fs‘)  n  ■  VG  (fs  —  fs')  —  iv  {fs‘)  G  (fs  —  rs') 


dr  s' 


Integral  eqn  for  PCr^)-  Discretization  -4  matrix  eqn.  Invert  to  find  P(rs)- 
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Problems  with  the  Surface  Helmholz  Integral  Equation 

•  is  singular  when  fs‘  =  fs 

•  At  certain  frequencies,  there  are  interior  resonances  corresponding 
to  eigenfunctions  which  are  zero  on  the  surface. 


There  is  no  unique  solution  to  the  SHIE;  the  matrices  are  singular  and  cannot 
be  inverted! 

There  are  methods  for  avoiding  the  interior  resonances;  e.g.  the  Combined 
Helmholz  Integral  Eqn.  Formulation  (CHIEF),  where  the  field  at  selected 
points  in  the  interior  are  required  to  be  zero.  All  such  methods  are  messy, 
because  they  require  user  intervention  to  select  interior  points,  etc.,  and  the 
computer  codes  are  long  and  complicated. 
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Solution  to  the  Arbitrarily  Shaped  Surface  Problem 

With  accessible  high-speed  computers,  there  is  an  Easy  Solution 
Index  surface  points  with  ‘f.  Let  Pj  =  P  (fsj) . 

Any  solution  to  the  Helmholz  eqn.  may  be  expanded  as  series  of  eigenfunc¬ 
tions,  e.g.  spherical  wavefunctions: 

=  e*’^*Pr  (COS 

I  m 

Let  the  index  u  represent  (/,m) .  Let  Mju  =  hi  (krj)  (cos0j)  .  Then 

Pj 


Invert  matrix  to  find  Au  =  [Mju]~^  Pj.  FINISHED! 

Added  Benefit!  By  using  eigenfunctions  which  only  radiate  outward,  the 
problem  of  interior  resonances  is  entirely  avoided! 


Why  not  before?  if  not  level  surface,  then  many  evaluations  of  hi,  etc., 
and  no  FFT’s.  Not  viable  for  on-site  noise-controi  engineers. 
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Determining  the  change  of  Q  in  RUS  due  to  fluid  loading 


i- 

Q  uE  dt 

E  =  —ahc  ^  +  r o-no-u 

n  V 

dE  1  f  f 

1 

=  2 

i 

MfjVj]  ViAi 
i 

~  2^^^  MijXpi'ipjAi 

»■  j 
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Determining  the  change  of  Q  in  RUS  due  to  fluid  loading 


Normalize  with 


pscsioA 

2E 
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Solving  Boundary  Value  Problems  with  Spherical  Wavefunctions 
J.  D.  Maynard.  The  Pennsylvania  State  University 


Assume  a  time  dependence  exp  {—iut)  and  a  complex  sound  pressure  field  P  =  Re  {P  (f)} 
+tlm  {P  (f)}  satisfying  the  Helmholz  equation  with  wavevector  k  =  u/c..The  complex 
particle  velocity  field  is 


Using  separation  of  variables,  we  find  a  complete  set  of  complex  linearly  independent 
functions  which  radiate  out  to  infinity  (an  important  boundary  condition).  With  three 
space  variables,  there  must  be  two  independent  constants  of  separation  and  two  mode 
labels.  Thus  we  can  write  for  these  basis  functions: 


V“$im  (r)  +  k^^im  {f)  =  0  (2) 

Any  solution  to  a  radiation  problem  may  be  written  as  a  linear  combination  of  these 
functions; 


=  (3) 

I  m 

Let  Aim  =  (iim  +  ibim  and  ^im  (0  =  Pim  i'f^  +  iSim  (0  5  where  the  new  constants  and  fields 
are  purely  real.  Now 


=  [(a/mP«m  (f)  -  bimSim  (0)  +  *  {aimSlm  (0  +  bimR{m  (0)]  (4) 

/  m 

The  known  boundary  condition  is  the  normal  component  of  the  particle  velocity  at  a 
surface  given  by  points  fg .  Note  that  the  known  surface  velocity  is  purely  real.  We  let  the 
known  data  by  given  by  the  real  function  defined  by 


f  (fs)  =  pc  h  (fg)  ■  V  (fs)  (5) 

where  h  (fg)  is  the  unit  normal  at  the  surface  point  fg.  Combining  Eqs.  1,  4,  and  5,  we 
have 


/k)  =  EE 


aim  ( rn  (f.)  ■  VS|„  (f,))  +  6im  ( In  (f,)  ■  (r,) 


248 


(6) 


-EE 


1 , 


^Im  1  ^  (^s)  *  ^ Rim  (^s)  )  ^/m  (  ;  ^  (^s)  *  ^ Rim  (^s) 


/  m  ^ 


=  EE  [(aimSim  (j's)  +  ^ImRlm  i'^s))  *  {^ImRlm  (^s)  (^s))] 


I  m 


where  R'l^  {fs)  =  (1/fc)  n  (fs)  •  VRim  {ts)  and  S[^  {fs)  =  {l/k)  n  (fg)  •  VSim  (r®) 


Since  /  (fs)  is  purely  real,  the  imaginary  part  of  the  left  hand  side  of  Eq.  6  must  vanish, 
and  we  have 


E!  E!  ^IrnRlm  s)  —  E!  E^  i^s)  (7) 

I  m  I  m 

Now  assume  that  the  normal  component  of  the  particle  velocity  field  is  known  at  a  finite 
set  of  discrete  points  fs  with  s  =  1, 2, 3,  ...N.  We  simplify  notation  by  writing  /  (fs)  =  fs, 
and  similarly  for  other  functions  of  fg.  We  use  the  symbol  u  as  an  index  similar  to  s.  We 
also  simplify  notation  by  letting  the  one  integer  index  the  two  subscripts  I  and  m  :  /j,  = 
[/  (n) ,  m  (/i)] .  We  truncate  to  a  finite  number  of  basis  functions,  so  that  //  =  1,  2, 3,  ...L. 
Let  u  be  an  index  similar  to  /x.  Finally,  we  use  the  convention  that  an  index  repeated  in 
a  term  indicates  a  sum  over  that  index.  Now  we  can  rewrite  Eqs.  6  and  7  as 


fs  —  (®) 

Rs/i^t^  ~  (^) 

Now  fs  is  a  vector  of  length  N,  and  6^  are  vectors  of  length  L,  and  and  5"'^  are  N 

by  L  matrices.  We  shall  have  N  >  L.  so  that  the  matrices  are  not  square.  The  method  of 
Singular  Value  Decomposition  (SVD)  allows  one  to  find  the  inverse  of  such  matrices.  Let 
T,yu  be  the  inverse  of  5'^.  Multiplying  Eq.  9  by  Tl/u  solves  for  the  vector  b  in  terms  of  the 
vector  a  : 


-  Ti,y,R^^(l^ 

Substituting  Eq.  10  into  Eq.  8  yields 


(10) 


fs  —  i^sfi  Rsu'R^uRu^i)  (^1) 

Use  SVD  a  second  time  to  find  the  inverse  of  the  matrix  in  parentheses  in  Eq.  11.  Indicating 
this  inverse  with  Qfi,s,  we  now  have 


dfj,  —  Qfisfs  (12) 

This  is  equivalent^  to  adjusting  the  coefficients  to  least-squares  fit  the  boundary  data 
fs-  With  Eq.  10  we  also  have  b^.  From  Eq.  4  we  have 
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Re  {P  (fs)}  =  Rs^Uf,  -  Ssfibf^  (13) 

Im  {P  (fs)}  =  Ss(j,a^  +  Rs^bfj,  (14) 

With  P  (fs)  and  f  (fg)  =  pch{fs)  ■  v{fs) ,  all  other  acoustic  quantities  may  be  calculated 
using  closed-form  expressions. 

We  shall  separate  variables  in  spherical  coordinates  (r,  d,  4>)  but  fg  and  n  (r*)  will  be  given 
in  cartesian  coordinates  [x,  y,  z) .  We  shall  make  all  lengths  non-dimensional  by  multiplying 
by  the  wavenumber  k  =  u/c.  We  simplify  notation  by  letting  kx  -y  x,  ky  y,  kz  -f 
z,  kr  r,  and  (1/^)  V  ->  V.  We  also  define  (  =  cos0.  We  have 


r  =  y/x"^  ^y"^  +  z^  <!>  =  tan  ^  [y/x]  C  =  V’' 


dr  x 

dr  y 

dr 

z 

dx  r 

dy  r 

dz 

T 

dd 

-y 

d(f>  _ 

X 

d(j) 

dx 

x"^  + 

dy  x^ 

+ 

to 

dz 

—xz 

_  -yz 

2  1 

X“  +7/' 

dx 

y.3 

dy 

^3 

dz 

r3 

With  h  {r)  =  (rii,  Uy,  n,) ,  we  make  the  following  definitions: 


(15) 


(16) 


dr  dr  dr 

(17) 

^  +  -f-n,— 

d(f)  d(j)  d4> 

(18) 

ac  ac  ac 

(19) 

We  now  have 


d 


d 


n.V  =  K^^+h,-  +  h,- 

Separating  variables  in  spherical  coordinates  yields  the  spherical  wave  functions: 


(20) 


^Im  {r.  $,  (j))  =  {ji  (r)  -t-  zyi  (r))  (cos  m<f)  -I-  i  sin  mcj))  Yim  (cos  6)  (21) 
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where 


(C)  = 


'2/  4- 1  (/  —  m)l 
Alt  {l  +  m)\ 


pm 


{<) 


(22) 


The  ji  and  yi  are  the  spherical  Bessel  functions,  and  the  P[^  are  the  associated  Legendre 
polynomials.  The  Yim  are  similar  to  the  spherical  harmonics,  but  with  the  0  dependence 
removed;  the  (j)  dependence  is  included  explicitly  in  Eq.  21. 


Note  that  —l<m<l.  Eq.  22  is  for  m  >  0  only;  for  m  <  0,  we  use  the  Condon  and 
Shortley  phase:  Yi-rn  (C)  =  (-1)'”  Yim  (C)  • 

We  truncate  the  number  of  basis  functions  by  letting  Lmax  be  the  highest  order  for  the 
spherical  Bessel  functions.  The  the  total  number  of  basis  functions  is 


•^max 

L  =  ^  {21  +  1)  =  (Lmax  +  1)^  (23) 

/=0 

The  real  and  imaginary  parts  of  are 


Rim  {r.  0,  C)  =  {ji  {t)  cos  m</»  -  yi  (r)  sin  mcj))  Yim  (C)  (24) 

Sim  (r,  </>,  C)  =  Ui  (r)  sin  m(f)  +  yi  (r)  cos  mif))  Yim  (C)  (25) 


We  now  have 


dRim 

dr 

dRim  ^ 

dcp 

dRim 

dC 


=  (j'l  (r)  cos  m<l)  -  y'l  (r)  sin  m.<b)  Yim  (C) 
-m  {ji  (r)  sin  m(t)  +  yi  (r)  cos  m4>)  Yim  (C) 
=  {ji  (r)  cos  m(f>  -  iji  (r)  sin  m(f))  Yl^  (C) 


r>i  f  j.  u  ^Rlm  I  7  ^Rlm  ,  dRim 


dr 


dtp 


dC 


and  similar  expressions  for  Sim  and  S',. 


/m* 


(26) 

(27) 

(28) 
(29) 


It  should  be  noted  that  is  singular  when  x  =  y  =  0.  Also,  (^)  is  singular  when 
^  =  ±1  (i.e.  when  x  =  y  =  0  ).  When  these  singularities  occur,  Yim  is  zero  and  is 
zero.  The  zeros  dominate  except  when  m  =  ±1.  After  some  algebra,  we  find  that  when 
2;  =  y  =  0(C  =  ±l  and  (p  is  undefined)  and  m  =  ±1,  then 
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(30) 


R'lm  (^>  0  =  -  J  {ji  (r)  no,  -  yi  (r)  muy)  (O) 

=  Ui  (0  -  yi  (r)  muy)  {I  +  1)  C'l 

s'lm  (-^’^  (’■)  +  2//  (r)  n^)  (^/l-c2y;;„  (o) 

=  (ii  (^)  ”^”3/  +  yi  (0  ’^x)  (y^  (/  +  1)  c']  (31) 

The  spherical  Bessel  functions  and  the  associated  Legendre  functions  may  be  calculated  by 
computer  using  the  upward  and  downward  recursion  relations  as  in  Ref.  1.  The  derivatives 
of  the  functions  may  also  be  found  by  computer  using  the  following  relations: 


fi  (r)  =  ji-i  (r)  -  (I  +  l)ji  (r)  jr 
and  when  (  ^  ±1  and  0  <  m  <  I, 


YL  (0  =  W  +  m)  Ki-i,™  (C)  -  (0) 


Note  that  if  |  m  |>  /,  then  Vim  (0  =  0.  In  a  computer  program  which  calculates  the 
associated  Legendre  functions  care  must  be  taken  when  C  =  ±1-  In  particular,  the  limiting 
value  indicated  in  Eqs.  30  and  31  should  be  used  rather  than  the  singularity. 

1.  W.  H.  Press,  S.  A.  Teukolsky,  W.  T.  Vettering,  and  B.  P.  Flannery,  Numerical  Recipes 
(Cambridge  University  Press,  Cambridge,  1992). 
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ABSTRACT: 

The  resonant  sphere  technique,  RST,  was  applied  to  measure 
elastic  moduli  of  spinel  single  crystal  sphere  with  3.6675  mm  in 
diameter.  Resonance  frequency  of  23  modes  were  measured  between 
293  and  1200  K.  As  similar  to  previous  measurement  for  a 
rectangular  prism  specimen  (Suzuki  and  Anderson,  1982),  resonance 
frequencies  of  all  modes  decrease  with  temperature,  and  some  of 
them  show  discontinuous  change  of  slopes  near  904  K,  whereas  some 
other  modes  do  not  show  the  chemge  clearly. 

Inversion  of  the  frequency  data  provide  us  a  set  of 
independent  elastic  moduli  of  and  C44  (or  another  set  of 

Cii,  C12  and  C44)  as  a  function  of  temperature.  The  standard 

error  of  observed  frequency  after  iterative  calculation  is  a 
=0. 16kHz  (ca.  0.012%).  Therefore,  estimated  probable  errors  of 
these  moduli  are  very  small  as  compared  with  those  in  previous 
measurements  with  RST. 

Elastic  moduli  of  the  spinel  specimen  with  density  ^=3.5846 
g/cm^  at  293K  (20  ®C)  and  their  probable  errors  are: 

Cii  =  281.31  ±0*01 

C12  =  155.43  ±0.01 

C44  =  154.588  ±0.007 

and  C„  =  62.936  ±0.003  in  GPa.  The  anisotropy  factor  is 

A=C44/Cg=2.4561,  which  is  much  larger  them  that  of  MgO  (A*1.55)  and 
CaO  (A=0.98).  Isotropic  properties  cure 
Ks=197.39  ±0.01  , 

and  ;r=107.81  (Hill  average). 
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Combining  with  thermal  expansion  data, 
anharmonic  parameters .  For  example : 

we  can  evaluate 

Gruneisen  parameter 

Anderson- 

-Gruneisen  parameter 

r  = 

= 

8^  * 

at  300K 

1.16 

2.98 

4.72 

1,200K 

1.17 

4.46 

6.37 

One  of  sheez*  moduli,  do  nof  have  clear'  discontinuity  in 

the  slope  against  temperature,  whereas  another  shear  modulus 
C3(  =  (Cii-Ci2)/2)  has  distinctive  bend  at  904  K,  and  bulk  modulus 

Kg  (  =  (Cii+2Ci2)/3)  has  intermediate  between  those,  because  of 
opposite  change  of  slopes  in  and  0^2  above  904  K. 

Deviation  at  904  K  from  assumed  linear  tendency  is  about  2  %  for 
Cg  and  less  for  the  other  moduli.  The  aJarupt  change  of  the 

Gruneisen  parameter  at  904K  is  Ay—  0.1,  mostly  due  to  leunbda 
transition  in  thermal  expansion  coefficient. 

Preliminary  MD  calculation  shows  that  the  second  order 
transition  may  be  caused  by  the  site  preference  of  magnesium  ions 
to  octahedral  sites  at  high  temperature. 


KEY  WORDS: 

resonant  sphere  technique,  RST,  anisotropy,  elasticity, 
spinel,  thermal  expansion 
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Introduction 

Elasticity  and  its  pressure  and  temperature  dependency  of 
minerals  is  greatly  important  in  study  for  constituents  and  their 
state  in  the  Earth  interior.  For  the  measurement  of  elasticity, 
we  take  notice  of  a  resonance  method  adopting  a  sphere  shape  for 
specimens,  no  matter  which  they  are  isotropic  or  anisotropic. 

This  method  may  be  called  the  resonant  sphere  technique,  or  RST 
and  is  advanced  in  these  years  for  anisotropic  materials .  The 
sphere  shape  of  specimens  in  the  resonance  method  brings  many 
advantages  in  measurement  of  elasticity  and  cuielasticity  (Suzuki 
et  al.,  1991,  1992). 

In  the  Earth's  mantle,  one  of  major  constituents  of 
orthosilicate,  olivine,  changes  its  structure  having  high  density, 
as  demonstrated  by  many  high  pressure  experiments.  The  phase  is 
called  as  the  gamma-phase  and  has  been  greatly  taken  notice 
of(e.g.,  Suzuki  et  al.,  1979-a).  Its  crystal  structure  is  the 
same  as  spinel  MgAl20^.  The  gamma-phase  is  not  stable  at  ambient 

pressure,  and  much  attention  has  been  focused  on  spinel,  MgAljO^ 

as  its  analogue 

Spinel  MgAljO^  is  also  one  of  important  minerals  for  practical 

uses,  and  many  physical  properties  have  been  measured.  A  clear 
inflection  point  is  found  at  933  K  in  the  linear  expansion  curve 
of  spinel  (Suzuki  and  Kvunazawa,  1982),  and  such  inflection  was 
also  found  in  electrical  conductivity  near  the  temperature  (Weeks 
and  Sender,  1980).  Crystal  structure  was  studied  extensively 
(e.g.,  Yamanaka  and  Takeuchi,  1982)  and  they  suggested  an  order- 
disorder  phase  change  due  to  preferred  site  occupation  of  cations . 

Temperature  change  of  elasticity  was  once  measured  for  the 
rectangular  prism  specimen,  and  the  similar  change  was  confirmed 
in  elasticity  change  (Suzuki  and  Anderson,  1982),  though  some 
systematic  deviation  of  frequency  had  brought  in  some  eonbiguity. 

In  the  present  measurements ,  we  used  a  synthetic  specimen  of 
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spinel  shaped  into  an  accurate  sphere.  Frecpiency  changes  of'  it 
were  traced  up  to  1,167  K  well  above  the  transition  temperature. 

Data  Reduction 

In  order  to  deduce  elastic  moduli  ^ij  from  resonant 

frequencies,  it  is  necessary  to  compute  accurately  enough  the 
eigenfreguencies  of  the  specimen,  since  the  deduction  is  made  by 
comparing  the  measured  resonant  frequencies  with  reference 
eigenfreguencies  given  by  the  elastic  moduli  assumed  on  the 
specimen.  Theoretical  evaluation  of  eigenfreguencies  for 
anisotropic  materials  is  described  elsewhere  (Michizuki,  1988; 
Isoda  et  al.,  1990) 

In  order  to  obtain  elastic  moduli  from  n's  observed  resonance 
frequencies,  the  residual  equations  are  defined  as: 


(k  “  1,  ^  3 . . .  n) 


where  4-  4°  is  difference  of  observed  frequency  4  and 

computed  reference  frequency  4°  k  -  th  mode,  = 

-  is  difference  of  specimen's  elastic  moduli  ^ij  and  assumed 

values  of  reference  elastic  moduli  (e.g.,  Chang  cind  Barsch, 

1973),  is  residual  of  the  k  -  th  mode.  Summation  is  made  for 

possible  numbers  of  i  euid  j  satisfying  a  relation  i  .  The  4° 
and  34  are  computed  by  using  the  reference  elastic  moduli 

,  dimension  and  mass  of  the  specimen.  The  are  the 

unknowns  to  be  determined  by  minimizing  the  squared  sum  of 
residuals,  Z  The  elastic  moduli  of  the  specimen  are  finally 

obtained  by  ''' 
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EKperimental 

Apparatus  and  specimen 

The  block  diagram  of  the  RST  measurement  system  is  shown  in 
Figure  1.  This  system  is  basically  not  different  from  those  shown 
in  Fig. 4  in  Ohno  (1976)  or  Figs.  1  and  2  in  Sumino  et  al.(1976), 
when  continuous  wave  (CW)  method  is  applied  for  data  acquisition. 
Recently  further  developments  were  made  to  realize  the  FT-method 
(Inouye,  1993?  Fujio,  1995)  and  high  pressure  measurement  (Ohno  et 
al.,  1993). 

A  synthetic  boule  of  spinel  was  cut  and  shaped  into  a  sphere 
with  the  so-called  two-pipe  method  (Fraser  and  LeCraw,  1964). 

For  accurate  determination  of  elasticity,  we  try  to  make  the 
relative  error  of  diameter  less  than  a  few  micro-meters. 

Measurements 

Resonance  frequency  under  no  external  forces  was  obtained  by 
measurements  changing  applied  force  to  zero  (Table  1).  At  higher 
temperature,  frequency  was  measured  under  constant  supporting 
force,  of  which  effect  was  corrected  with  the  frequency  of  zero- 
force  at  room  temperature. 

Temperature  was  measured  with  chromel-alumel  thermocouples  set 
very  near  to  the  specimen.  Measurements  of  frequency  were  made 
for  23  modes  between  room  temperature  293  K  amd  1173  K  (Figure  2). 
Temperature  intervals  were  about  30  K  around  the  transition 
temperature  and  some  tens  Kelvin  otherwise.  Reproducibility  was 
confirmed  by  measurements  at  several  temperatures  of  the  cooling 
run  made  at  the  close  temperatures  of  the  heating  run.  They  do  not 
show  much  difference  from  the  heating  run.  Hysteresis  is  quite 
small,  if  possible. 

Figure  3-a  and  -b  compare  the  temperature  dependency  of 
frequency  near  the  temperature  of  904  K.  Some  of  modes  change 
clearly  their  slopes  in  both  sides  of  the  transition  temperature. 
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Results 

Elasticity 

Elasticity  of  the  spinel  specimen  is  obtained  by  the  least 
squares  calculation  mentioned  above  and  listed  in  Table  2. 
Aggregate  isotropic  properties  are  obtained  by  the  Voigt-Reuss- 
Hill  averaging  scheme  (Kumazawa,  1964)  and  also  listed.  The 
iterative  calculation  was  made,  and  this  could  reduced  the 
standard  error  very  much  as  shown  in  Table  2. 

Temperature  dependency  of  elastic  moduli  are  shown  in  Figures 
4  and  5.  In  the  cubic  symmetry,  three  of  elastic  moduli  are 
independent.  Least  squares  calculations  using  the  modes  with 
clear  change  of  slopes  and  without  those  show  that  results  are 
very  close  each  other  within  errors,  and  this  confirms  the 
consistency  of  observed  frequency  data. 

Anisotropy  factor 

Anisotropy  factor  of  cubic  symmetry  is  defined  by  the  ratio  of 
two  shear  moduli  as 
A=  C44/C3 

where  Cs  is  a  shear  modulus  defined  by  Cg=(Cj^2  -^2)  /2.  The  A  of 

spinel  is  about  2.5.  For  isotropic  materials,  A  is  unity,  and  MgO 
and  CaO  have  1.11  and  1.22  respectively.  Elastic  properties  of 
spinel  is  said  to  be  quite  anisotropic.  Also,  the  A  shows  clear 
inflection  at  904  K  as  shown  in  Figure  6. 

Thermal  expansion 

Thermal  expansion  data  of  spinel  is  reproduced  in  Figure  7. 
Smoothed  values  of  expansion  coefficients  are  obtained  applying 
the  Gruneisen’s  theory  of  thermal  expansion  (Suzuki  et  al.,  1979- 
b) .  Expansion  coefficients  just  below  the  transition  are 
expressed  by  a  quadratic  equation  as  shown  in  Figure  7.  Inflec¬ 
tion  temperatures  from  thermal  expansion  measurement  and  frequency 
measurement  are  a  little  different,  probably  due  to  differences  in 
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measurement  system  and  specimen.  In  the  calculation  using 
thermal  expansion  data,  this  small  discrepancy  is  adjusted  by 
assignment  of  o- (rj.)=28.6X10"^/K  and  a  (To)=25-6X IQ-^/K,  where 

is  the  transition  temperature,  Tq  =  800  K,  arid  a  function  a  (r- 
Tf^)^+b  (r-r(j)  +  aQ  is  assumed  between  Tq  and  T^..  The  Gnineisen  para¬ 
meter  is  evaluated  combining  bulk  modulus,  thermal  expansion  and 
heat  content  data  (Robie  et  al.,  1978),  and  shown  in  Figure  8. 

Discussion 

In  the  present  measurements,  we  can  obtain  accurate  values  of 
elastic  moduli  of  spinel,  and  also  we  can  trace  temperature 
dependency  of  them  precisely.  Much  of  these  are  due  to  the 
sphere  shape  of  the  specimen.  The  xyz  algorithm  (Visscher  at  al., 
1991)  is  useful  for  iterative  calculation,  when  combined  with  the 
variational  method  for  evaluation  of  df  (Oda  et  al.  1993). 

This  reduced  the  standard  error  very  much  as  shown  in  Table  2. 

This  magnitude  of  error  is  much  smaller  than  those  in  previous 
reports  by  the  resonance  methods. 

Two  shear  moduli  of  spinel,  ^44  cind  Cg  are  quite  different  in 

magnitude  as  shown  in  Figure  5  (b)  and  (d)  and  in  Figure  6,  eind 
this  causes  large  elastic  anisotropy  as  compared  with  other  cubic 
crystal,  such  as  MgO  and  CaO.  The  Cg  is  less  than  a  half  of 

which  means  the  structure  is  very  compliant  to  shear  stress 
parallel  to  [110]  direction  cuid  this  may  be  also  the  reason  of 
large  temperature  effect  on  Cg  and  the  phase  transition. 

Preliminary  calculation  with  the  MD  method  shows  that 
magnesium  ion  in  tetrahedral  sites  is  easy  to  move  to  octahedral 
sites  at  high  temperature  (personal  communication,  Matsui  1997). 
This  may  be  the  reason  of  the  transition  at  about  900  K. 

Both  of  observed  thermal  expansion  coefficient  a  and 
isothermal  bulk  modulus  have  discontinuities  at  the  transition 
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temperature  under  ambient  pressure.  The  boundary  of  the  low  and 
high  temperature  phases  can  be  estimated:  difference  of 
isothermal  bulk  modulus  at  the  transition  temperature  is  1  GPa  or 
smaller,  the  boundary  is  about  100  GPa  or  higher  at  room 
temperature . 

Acknowledgement :  The  authors  are  grateful  to  Dr.  H.  Oda  and  Mr- 
K.  Yamada  for  the  computation,  and  Messrs.  S.  Sawane  and  K.  Ohno 
for  their  help  in  high  temperature  measurements .  Discussion  with 
Dr.  A.  Yoneda  helps  improve  the  contents. 
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Table  1.  Observed  frequency  of  the  spinel  specimen 
with  diameter  3.6675mm. 


no. 

mode 

f®/kH2 

1 

Eg  -1 

977.25 

2 

T2u-1 

1003.16 

3 

Eu  -1 

1390.13 

4 

Tlu-1 

1414.93 

5 

A2g-1 

1432.84 

6 

T2g-1 

1465.34 

7 

T2g-2 

1748.48 

8 

Tlg-1 

1834.20 

9 

Tlu-2 

1919.90 

10 

Eg  -2 

1933.64 

11 

T2u-2 

1987.20  * 

12 

T2u-3 

2096.48 

13 

Alg-1 

2107.10 

14 

A2u-1 

2202.08 

15 

Tlg-2 

2238.40 

16 

Eu  -2 

2280.10  = 

17 

Alg-2 

2280.10  = 

18 

Tlu-3 

2400.88 

19 

Eg  -3 

2435.85 

20 

T2g-3 

2507.80 

21 

T2u-4 

2546.30 

22 

Tlg-3 

2608.55 

23 

T2g-4 

2652.35 

24 

Eg  -4 

2691.70 

25 

Alu-l 

2708.00 

*:  not  observed  at  high  temperatures 
=:  cannot  be  distinguished 
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Table  2.  Elasticity  of  single  crystal  spinel  MgAljO^ 
with  density  p  =3584. 6  Mg/m^, 
and  its  aggregate  properties,  at  293K. 

Single  crystal  elastic  properties  in  6Pa 


Ks 

197.395 

± 

0.013 

Cs 

62.936 

± 

0.003 

^44 

154.587 

± 

0.007 

•^11 

281.310 

± 

0.014 

^12 

155.437 

± 

0.013 

standard  deviation  in  the  least  squares 
calculation  a=0.16/kHz 


Isotropic  properties  of  the  Hill  average 
bulk  modulus  Kg=197.39  GPa 

rigidity  modulus  jU =107. 806  GPa 
Young's  modulus  E  =273.61  GPa 
Poisson's  ratio  =0.26898 
elastic  wave  velocity  Vp  =9.755  km/s 

Vg  =5.484  km/s 
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Figure  1.  Computer-controlled  system  for  data  acquisition  in 
the  Resonant  Sphere  Technique.  Thick  arrows  show  computer 
control  of  the  devices  and  data  acquisition  through  the 

Thin  arrows  show  transmission  of  analogue  signal. 
A/D:  Analogue-digital  converter,  P/C:  Personal  computer, 

P/T:  Pressure  and/or  temperature  sensor,  O/P:  Oven  and/or 
pressure  vessel.  S:  Specimen,  T:  Transducer,  x-y/R:  x-y 
recorder,  M/F:  Main  Frame  for  the  least  squares  calculation 


Figure  2.  Observed  frequencies  of  23  modes  between  293  and  1167  K 


Frequency  f/kHz  Frequency  f/kHz 
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Figure  3  .  Magnified  illustration  of  Figure  2  auround  1,400  and  2,500 
kHz.  (a)  From  lower  to  higher  frequency;  Eu-1,  Tlu-1,  A2g-1  and 
T2g— 1  modes.  Tow  modes  Tlu— 1  and  A2g— 1  show  clear  inflection 
at  904  K,  and  the  others  are  almost  linear  with  convex  curva¬ 
ture.  (b)  From  lower  to  higher  frequency;  Tlu-3,  Eg-3,  T2g-3  and 
T2u-4  modes.  The  modes  Tlu-3  and  T2u-4  show  clear  inflection  at 
904  K. 


Elastic  moduli,  M/GPa 
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Figure  4,  Elastic  moduli  of  spinel  in  the  temperature  range 
between  293  and  1167  K.  Three  of  those  are  independent. 
Adiabatic  and  isothermal  bulk  moduli  are 

Kp  =K  /(1+ff  7  T) ,  respectively.  Shear  modulus  is  C  =(C  - 

^  B  8  11 

0^2) Rigidity  modulus  are  obtained  by  the  Voigt-Reuss- 
Hill  averaging  schemes - 
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Figure  5.  Magnified  illustration  of  Figure  4. 

Thermal  expansion  data  and  the  Gruneisen  parameter 
are  used  to  evaluate  isothermal  bulk  modulus  in 
the  panel  (c).  <a)  (b)  ^44 

(c)  bulk  moduli,  and  (d)  sheeu:  modulus  C^, 

and  (e)  rigidity  modulus  fi  . 
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Shear  modulus,  Cs/GPa  Bulk  moduli,  K/GPa 


POSTER  18 


270 


Anisotropy  factor  A  Rigidity  modulus  u/GPa 
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Figure  6.  Anisotropy  factor  A  changes  its  slope  clearly 
at  904  K.  The  A  is  defined  as  A=C../C 

44  S 
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Figure  7*  Thermal  expansion  coefficient  of  spinel  in  10^-6/K* 

square:  obtained  by  numerical  differentiation  of  observed  expansion 
data;  smooth  curve:  by  least  squares  fit  of  the  Gruneisen  equation; 
dot:  by  a  quadratic  equation  just  before  the  transition. 
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Resonant  frequency  of  an  elastic  sphere  at  high  pressures 
Hitoshi  Oda,  Kouichi  Yamada,  and  Isao  Suzuki 
Department  of  Earth  Sciences,  Okayama  University,  Okayama  700,  Japan 
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1.  Introduction 

In  order  to  extend  pressure  range  for  application  of  the  resonant 
sphere  technique  (RST),  the  cavity  resonance  method  has  been 
developed  by  Ohno  et  al.  (1993).  In  this  method,  a  spherical  specimen 
is  placed  in  the  center  of  a  spherical  container  which  is  filled  by  high 
pressure  gas,  and  resonant  frequencies  of  the  spherical  specimen  are 
measured  under  the  confining  pressure.  Since  the  elastic  constants  are 
determined  by  minimizing,  in  a  least  squares  sense,  the  differences 
between  measured  and  computed  resonant  frequencies  of  the  specimen, 
the  free  vibration  frequencies  have  to  be  computed  for  three  layered 
sphere  which  is  composed  of  the  inner  spherical  specimen,  high 
pressiu'e  gas  layer  and  outer  solid  layer.  However  the  free  vibration  of 
a  radially  heterogeneous  sphere  having  a  gas  layer  has  not  been  known 
in  detail  so  far.  In  the  present  study,  we  investigate  how  the  resonance 
spectrum  of  the  inner  spherical  specimen  is  affected  by  the  gas  layer, 
solving  a  problem  of  free  vibration  of  the  three  layered  sphere. 

The  spherical  shell  structure  composed  of  three  layers  is  shown  in 
Figure  1.  The  outer  layer  is  a  solid  container  of  tungsten  carbide,  WC, 
and  helium  gas,  He,  fills  the  middle  layer  as  pressure  medium  to 
compress  an  inner  spherical  sample.  This  model  is  called  “cavity 
model”.  The  problem  of  free  vibration  was  solved  for  two  cavity 
models,  where  the  inner  spherical  specimen  is  elastically  isotropic 
(model  A)  and  anisotropic  (model  B).  The  elastic  constants,  density  and 
radius  of  both  models  are  given  in  Table  1.  The  computer  program 
coded  by  Saito  ( 1988)  was  used  for  computing  the  frequencies  and 
eigenfunctions  of  normal  modes  of  the  three  layered  sphere  with 
isotropic  elasticity. 
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2.  Free  Vibration  for  cavity  model  with  an  isotropic  specimen 

Resonance  modes  for  model  A  are  described  by  torsional  nT]  and 
spheroidal  nS/  modes,  where  n  and  1  are  the  radial  and  angular  order 
numbers,  respectively.  These  normal  modes  are  characterized  on  the 
basis  of  the  distribution  of  strain  energy  excited  in  the  three  layered 
sphere.  The  strain  energy  of  torsional  modes  is  confined  only  in  the 
inner  sample  or  the  outer  solid  layer  because  of  no  excitation  of  the 
torsional  displacements  in  the  gas  layer.  However  the  energy  of 
spheroidal  modes  is  generally  excited  over  the  entire  space  of  the  three 
layered  sphere.  Figure  2  depicts  the  sti*ain  energy  distribution  of 
spheroidal  modes.  The  strain  energy  is  localized  in  the  inner  sample, 
gas  or  outer  layer  because  of  a  large  difference  in  acoustic  impedance 
between  the  gas  and  solid  layers.  We  call  such  a  mode  “sample”,  “gas” 
or  “container”  mode,  whose  nomenclature  derives  from  the  layer  where 
the  strain  energy  is  localized.  In  these  modes,  the  frequencies  of  sample 
modes  provide  information  regarding  the  elastic  constants  of  the  inner 
spherical  specimen. 

The  resonant  frequencies  of  the  thre(j  layered  sphere  were  computed 
as  a  function  of  pressure  in  gas  layer.  The  changes  in  elastic  constants 
and  density  due  to  pressure  and  the  volume  change  of  the  sample  and 
container  were  taken  into  account  for  the  computation  of  resonant 
frequency.  The  pressure  derivatives  of  elastic  constants  are  listed  in 
Table  1.  In  Figure  3,  the  pressure  dependence  of  the  computed 
frequencies  of  the  ^52  modes  for  mod<jl  A  is  shown  for  the  sample,  gas 

and  container  modes.  The  number  of  gas  modes  is  largest  in  the  three 
kinds  of  resonance  modes.  The  gas  mode  shows  the  largest  slop  in  the 
pressure  change  of  resonant  frequency,  and  its  pressure  curve  often 
intersects  those  of  sample  and  container  modes.  For  a  careful 
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examination  of  the  intersections  of  the  pressure  curves,  a  small 
rectangular  part  in  Figure  3,  where  the  pressure  curves  of  sample  and 
gas  modes  appear  to  cross  at  the  presstu'e  of  22.3  MPa,  is  enlarged  in 
Figure  4.  The  ciuve  of  10S2  mode  approaches  that  of  11S2  mode  as  the 

confining  pressure  increases  up  to  the  intersection  pressime,  and  both 
modes  gradually  separate  with  further  increase  in  pressme.  Figure  5 
depicts  the  change  of  strain  energy  distribution  of  the  10^2  and  nS2 
modes.  As  the  pressure  increases,  the  10S2  changes  from  gas  mode  to 
sample  mode  and  the  11S2  changes  from  sample  mode  to  gas  mode. 

When  om  attention  is  focused  on  the  sample  or  gas  mode,  the 
intersection  is  regarded  as  an  interchange  from  11S2  to  10S2  for  the 
sample  mode  and  from  10S2  to  11S2  for  the  gas  mode.  Such  an 
interchange  occurs  only  between  the  spheroidal  modes  of  nS;  and 
satisfying  the  conditions  of  i=i'and  n-n'*±l. 
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Vibration  for  cavity  model  with  an  anisotropic  specimen 
In  the  case  of  model  B  in  which  the  inner  spherical  specimen  is 
elastically  anisotropic,  the  free  osciBation  of  three  layered  sphere  can  be 
computed  by  the  Rayleigh-Ritz  (variational)  method,  where  the  normal 
mode  displacements  of  spherical  shell  structure  with  an  isotropic 
specimen  are  used  as  the  basis  functions.  Figure  6  shows  the  frequency 
spectrum  computed  for  the  cavity  model  with  a  spherical  periclase 
specimen.  The  resonance  modes,  as  well  as  those  for  model  A,  are 
classified  into  sample,  gas  and  container  modes  on  the  basis  of  the  strain 
energy  distribution.  It  is  foimd  that  most  of  resonance  modes  are  gas 
modes.  For  a  comparison,  the  frequency  spectrum  computed  for  the 
cavity  model  with  an  isotropic  specimen  is  also  shown  in  Figure  6.  The 
elastic  constants  of  the  isotropic  specimen  were  obtained  by  the  Voigt 
averaging  scheme  of  the  elastic  constants  of  the  periclase.  A  larger 
number  of  sample  modes  are  identified  in  a  given  frequency  interval,  in 
comparison  with  free  vibration  for  the  cavity  model  with  the  isotropic 
specimen.  This  increase  in  sample  modes  is  attributed  to  split  of 
degenerate  modes  due  to  the  elastic  anisotropy  of  the  inner  specimen. 
The  pressure  dependence  of  resonant  frequencies  for  model  B  is  shown 
for  mode  groups  peculiar  to  the  cubic  crystal  symmetry  to  which  the 
priclase  specimen  belongs  (see  Figure  7).  The  interchange  between  the 
gas  mode  and  sample  or  container  modes  occurs  only  in  the  resonance 
modes  belonging  to  the  same  mode  group.  The  resonant  frequencies  of 
sample  and  container  modes  increase  only  by  about  IkHz  as  the 
confining  pressure  increases  up  to  60MPa,  while  the  increase  in 
frequency  of  gas  mode  is  about  200kHz  for  the  same  pressure  increase. 
The  large  frequency  increase  of  gas  mode  may  be  attributed  to  larger 
change  rate  in  bulk  modulus  of  gas  layer. 
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Isotropic 


Sample  mode 

Gas  mode 
Containiermode 


Figtire  6.  Resonance  spectrum  (bottom  figure)  for  model  where  a  peiidase  specimen  is  used  as  die  inner 
spherical  specimen  of  three  layered  modeL  For  a  comparison,  frequency  spectrum  computed  for  the  cavity 
model  with  an  isotropic  specimen  is  shown  in  the  top  figure.  The  elastic  constants  of  the  isotropic  specimen 
were  obtained  by  the  Voigt  averaging  scheme  of  the  elastic  constants  of  the  periclase.  The  long,  middle  and 
short  bars  represent  the  resonant  frequencies  of  sample,  gas  and  container  modes,  respectively.  A  larger 
number  of  sample  modes  are  found  in  a  given  frequency  interval,  in  comparison  with  the  fiee  oscillation  for 
the  cavi^  model  with  the  isotropic  specimen.  This  increase  in  the  number  of  sample  modes  is  attributed  to 
split  of  degenerate  modes  due  to  the  elastic  anisotropy  of  the  irmer  sample. 
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4.  Discussion 

As  seen  in  Figures  3, 4  and  7,  the  intersection  between  the  pressure 
curves  makes  the  mode  identification  difficult  because  there  are  two  or 
more  resonance  modes  in  the  vicinity  of  the  intersection.  In  order  to 
avoid  such  a  trouble  in  the  mode  identification,  it  is  necessary  to  make  a 
vibration  system  that  the  number  of  gas  modes  is  as  small  as  possible  in 
a  given  frequency  interval  because  the  intersection  occurs  between  gas 
mode  and  sample  or  container  mode.  Figure  8  shows  the  effect  of 
thickness  of  gas  layer  on  the  number  of  gas  modes.  The  frequency 
difference  between  adjacent  gas  modes  is  larger  as  the  thickness  of  gas 
layer  is  thinner.  This  result  implies  that  the  possibility  of  intersections 
between  gas  and  sample  modes  may  be  lower  for  the  model  of  three 
layered  sphere  with  thinner  gas  layer. 

The  resonant  frequencies  were  measured  for  a  spherical  periclase 
specimen  by  the  cavity  resonance  method.  The  measured  frequencies 
are  compared  with  the  computed  frequencies  in  Figure  9.  Agreement 
between  measured  and  computed  frequencies  are  good.  This  result 
indicates  the  justification  of  our  computation  for  free  vibration  of  an 
elastically  anisotropic  sphere  under  high  pressure.  Therefore  the  RST  is 
applicable  to  elasticity  measurements  of  anisotropic  materials  under 
high  pressure. 
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5.  Condiision 

In  order  to  make  it  possible  to  employ  the  RST  for  the  elasticily 
measurements  under  high  pressure,  we  computed  free  vibration  for 
cavity  models  of  spherical  shell  structure  composed  of  three  layers  of 
solid  sample,  gas  and  solid  container.  The  results  obtained  are 
summarized  as  follows. 

l)It  is  possible  to  compute  the  free  vibration  of  the  three  layered 
sphere.  A  good  agreement  is  obtained  between  computed  and 
measured  resonant  frequencies. 

?.)The  resonance  modes  for  the  cavity  model  are  divided  into  sample, 
gas  and  container  modes,  v^ose  strain  energy  distributions  are  localized 
into  sample,  gas  and  container  layers,  respectively. 

3) Interchange  takes  place  between  adjacent  resonance  modes  belonging 
to  the  same  mode  group,  and  it  makes  the  mode  identifrcadon  difficult. 

4) To  avoid  such  an  interchange  between  resonance  modes,  it  may  be  a 
better  way  to  make  the  thickness  of  gas  layer  as  thin  as  possible. 
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Figure  1.  Model  of  spherical  shell  structure  that  is  composed  of  sample, 
gas  and  container  layers.  Radii  of  sample,  gas  and  container  layers  are 
denoted  by  r^,  r2  and  respectively.  The  elastic  constants  and  density  of 

each  layer  are  given  in  Table  1.  We  consider  two  cases  where  the  inner 
spherical  sample  is  elastically  isotropic  (model  A)  and  anisotropic  (model  B. 
In  model  B,  the  inner  specimen  is  a  single-crystal  periclase. 
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Tabel  1.  Patamater  values  of  model  B 


Density  (g/cm^)  Radius  (mm) 


Elastic  (Constants  (GPa) 


K 

Cji 

Cl2 

C44 

Sample 

Gas 

Container 

MgO 

He 

WC 

3.589 

2.5E-04 

15.1 

2.057 

3.057 

7.057 

162.83 

2.4E-04 

391.00 

133.69 

219.00 

296.15 

95.35 

155.89 

dK/dP 

dfi/dP 

aCii/dP 

dCii/dP 

dCii/3P 

Sample 

Gas 

Container 

MgO 

He 

WC 

4.22  2.08  8.98  2.84  2.09 

(a=0.99239f5.4467P.24.270p2+80.028  P^) 

6.0  Z5  -  ■ 

a:  Sound  velocity  in  km/s. 

The  elastic  constants  of  MgO,  K  and  were  obtained  by  the  Voigt  averaging  scheme  of  Qj. 


Table  1.  Paramater  values  of  model  A 


Density  (g/cnP) 

Radius  (mm) 

piagric  Constants  (GPa) 

K 

H  Cn  C12  C44 

Sample 

Gas 

Container 

Steel 

He 

WC 

7.789 

2.5E-04 

15.1 

2384 

5.0 

8.5 

165.0 

2.4E-04 

391.0 

79.5  -  -  ■ 

219.0 

dK/dP 

dfl/dP  dCii/dP  aCii/aP  dCii/dP 

Sample 

Gas 

Container 

Steel 

He 

WC 

99239+5.4467P-24.270P^+80.028  P®) 

6.0  2.5  -  - 

a:  Sound  velocity  in  km/s. 
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Figure  2.  Distribution  of  strain 
energy  excited  in  the  three 
layered  sphere  with  an 
isotropic  specimen  (model  A). 
The  strain  energy  distributions 
for  rigidity  and  bulk  moduli  are 
represented  by  the  solid  and 
dashed  lines,  respectively.  The 
strain  energy  is  localized  in  the 
inner  sample,  gas  or  outer  layer 
because  of  a  large  difference  in 
acoustic  impedance  between 
the  gas  and  solid  layers.  Such  a 
mode  is  called  “sample”,  “gas” 
or  “container”  mode,  whose 
nomenclature  derives  from  the 
layer  where  the  strain  energy 
is  localized.  In  these  modes, 
the  frequencies  of  sample 
modes  provide  information 
regarding  the  elastic  constants 
of  the  inner  sample. 
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Figure  3.  Pressure  dependence  of  resonant  frequencies  of  sample  (solid 
line),  gas  (dashed  line)  and  container  (dotted  and  dashed  line)  modes,  which 
are  equivalent  to  the  ^82  modes  of  model  A.  The  number  of  gas  modes  is 
largest  in  the  three  kinds  of  resonance  modes.  The  gas  mode  shows  the 
largest  change  in  the  pressure  dependence  of  the  resonant  frequencies  and 
its  pressure  curve  often  intersects  those  of  sample  and  container  modes.  A 
small  rectangular  area  is  enlarged  in  Figure  4. 
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Figure  4.  Pressure  dependence  of  the  frequencies  of  10S2  (solid  circles)  and 
11S2  (open  circles)  modes  in  the  pressure  interval  from  20  to  25  MPa.  The 
pressure  curve  of  ioS2  mode  approaches  that  of  11S2  mode  as  the  confining 
pressure  increases  up  to  about  22.3  MPa,  and  then  both  curves  gradually 
separate  with  further  increase  in  pressure.  When  our  attention  is  focused 
on  the  sample  or  gas  mode,  the  intersection  seen  in  the  small  rectangular 
area  of  Figure  3  is  regarded  as  an  interchange  from  11S2  to  10S2  for  the 
sample  mode  and  from  20S2  to  11S2  for  the  gas  mode.  Such  an  interchange 
occurs  only  between  the  spheroidal  modes  of  nS;  and  n'S/' satisfying  the 

conditions  of  1-1'  and  n-n'=±l,  where  n  and  n’  are  the  radial  order  numbers 
of  the  spheroidal  oscillation,  and  1  and  V  are  the  angular  order  numbers. 
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Figure  5.  Change  of  the  strain  energy  distribution  of  10S2  and  11S2  modes 
against  confining  pressure.  The  resonant  frequencies  of  both  mndes  are 
nearly  equal  at  about  22.3  MPa.  As  the  pressure  increases,  the  jqSz  changes 

from  gas  mode  to  sample  mode  and  the  11S2  changes  from  sample  mode  to 
gas  mode. 


POSTER  15 


E,-! 

40  60  so 

//MP» 


Rgure  7.  Pressure  dgjendence  of  resonant  firequendes  for 
model  B.  The  top  panel  shows  the  pressure  dq>endence  of 
the  frequencies  of  gas  (dashed  lines),  container  (thin  lines) 
and  1^-lf  T2u-1»  Eq«1  and  T2^-l  modes  (thick  lines)  which 
are  namenriamres  peculiar  to  free  oscillation  of  the  periclase 
specimen  that  belongs  to  cubic  crystal  symmetry.  The 
pressure  dependence  is  separately  represented  for  the  mode 
groups  of  T2u  and  T2^  in  the  lower  three  ptaneis.  The 
interchange  between  the  ga.^  mode  sample  or  container 
modes  occurs  only  in  the  resonance  modes  bdonging  to  the 
same  mode  groim. 
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Figure  8  The  effect  of  the  gas  lay^s  thickness  on  the  frequency  of  gas 
mode.  Pressure  dependence  of  the  frequency  is  shown  for  different 
thickness  of  the  gas  layer,  dgas-  The  left,  middle  and  right  panels 
correspond  to  the  cases  of  dgas=2.85,  =1.0  and  =0.1mm,  respectively.  As  the 
thickness  of  gas  layer  is  thinner,  the  frequency  difference  between  adjacent 
gas  modes  is  larger  and  die  number  of  gas  modes  decreases  in  a  given 
frequency  interval.  This  result  implies  that  the  possibility  of  intersections 
between  gas  and  sample  modes  may  be  lower  for  a  model  of  three  layered 
sphere  with  thinner  gas  layer. 
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0  20  40  60 

P/MPa 


Figure  9  Comparison  of  calculated  frequencies  with  resonant  frequencies 
measured  for  a  spherical  pehclase  specimen.  The  thick,  thin,  and  dashed 
lines  represent  the  pressure  dependence  of  die  resonant  frequencies  of 
sample,  container  and  gas  modes,  respectively.  The  measured  frequencies 
of  Eu’l  and  Eg-1  modes  are  shown  by  circles.  The  relative  difference 
between  both  frequencies  is  about  0.4%.  Frequency  increase  of  the  sample 
mode  due  to  pressure  changing  from  0  to  60  MPa  is  about  1  kHz,  while  the 
increase  in  frequency  of  gas  mode  is  about  200kHz  for  the  same  pressure 
increase.  If  the  frequency  increases  are  measured  for  a  number  of 
resonance  modes,  die  pressure  derivatives  of  the  specimen's  elastic 
constants  are  estimated. 
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Acoustical  resonator  frequency  shift  due  to  the 
migration  of  suspended  particles 
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Particle  Migration: 


A/2 
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Microsphere  CCD  Photographs: 
Glass: 


^  •'ift  tiOl: 
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Probe  beam 

i  ^ 

Layers  of  microspheres 


Reflected  beam 


Collinear  Four-Wave  Mixing: 

•  At  the  Bragg  condition  layered  particles  scatter 
sound  with  the  same  phase  [modulo  (2n)]. 

•  Equilibrium  volume  fractions  as  small  as  0.0005  can 
cause  appreciable  reflectivities  (or  frequency  shifts) 
once  the  particles  migrate  to  nodes  or  antinodes. 

•  At  the  Bragg  condition  both  the  pump  and  probe 
waves  are  In  the  "stop  band"  of  the  spatially  periodic 
medium  caused  by  the  particle  migration. 
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Number  density: 


|“p(z-‘)  =  D^np(z,t) 


2^1(„,(z,.)  d„(2b)) 


\ 


Diffusion 


Radiation  Pressure 


Steady  state  solution: 


n„(z) = 


k  Hq  exp[-q  cos(2kz)] 


n  Io(q) 


VpPoPor, 


SksT  1,  Poj”  k^T 


Acoastic  EnaCgY 


a 

Particle  Radius 

Solvent  Compressibility 

k 

Pump  Wavenumber 

Particle  Compressibility 

nD(z,t) 

#  Particles  A^olume 

SB 

Solvent  Viscosity 

no 

Mean  #  of  Particles/Area 

kfi 

Boltzmann’s  Constant 

Temperature _ 


Particle  Volume 


Standing  Wave  Energy 
Density  _ 


Pump  Pressure  Amplitude 


Particle  Diffusivity 


Compressibility  Factor 
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Number  Density: 
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Reflection  coefficient  for  N  layers: 

Bom  approximation  ->  Incident  wave  is  not  depleted 


_  sin[Nkd] 
sin[kd] 


~80  layers 
glass  microspheres 


pump  ^ 


*9 
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Received  Probe  Amplitude  Received  Probe  Amplitude 
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Normalization  of  the  Received  Probe  Wave: 


Source 


•~r“ 

- 

_ 

-i_ 

— 

— 

f  1 

— 

I 


— 

_ A_ 

__ 

_ 

—j 

1 

First  Interval  =  2  Paths 


Receiver 


Second  Interval  =  4  Paths 


TRANSPARENCY  11 


Frequency  Shift  Experimental  Setup 
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Gated  pump  amplitude  frequency  shift: 


plastic  microspheres 
1.05  gm/cc  sugar  solution 
volume  fraction  =  0.001 
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Ramped  pump  amplitude  frequency  shift: 


Pimp  Amplitude  (kPa) 


plastic  microspheres 
1.05  gm/cc  salt  solution 
volume  fraction  =  0.001 
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Shift  in  Frequency  from  Excluded  Volume 

Maximum  Kinetic  Energy  _  ^  ^ 
Maximum  Potential  Energy 


=  JlV(t)l^dv/|l(|)l^dV  =  ( ~) 


Z 


Rigid  particle 
of  volume  Vp 

Chamber  volume  = 


Scattering  correction:  T.  Wang  et  al.  JASA  (1982). 
Adiabatic  Invariance/radiation  pressure:  Putterman  (1989). 
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Frequency  Shift: 
Adiabatic  invariance: 
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.Eo 


dE  =  ^dCD 

(Do 


Change  in  frequency  for  a  single  particle: 


Peyilu-lu  W) 


((Oo-(»i)  =  fa[U(zo)-U(Zi)] 

■*^0 


Potential  energy  felt  by  particle:  GoT  fCOV  ( l962) 


E  V  r  3  ^ 

U(z)  =  =^--2.  cos(2kz) 

2  V..  3  J 


Av_  iVpr  3, 


2V. 


1-^  [cos(2kzo)-cos(2kzi)] 


Change  in  frequency  for  a  distribution  of  particles: 


Av  _  1  Vp 

Vn  2  L 


3  ^  /  7  N 

J  “p(z)cos  27t-  dz 

PoJ-L/2  l; 


System  energy 


Particle  volume 


Grating  spacing 


Resonator  volume 


Unperturbed  frequency 
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Maximum  frequency  shift: 

np(z)  =>  ^5(z  — Zjj) 
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(|)i  (z)  *  A  exp(ikiz)  +  B  exp(-ikiz)  <j)2(z)  =  C  exp(ik2z)  +  D  exp(-ik2Z) 


"1  I 

z 


Continuity  of  pressure  and  velocity  at  boundary: 


-i(opi(|),(0)  =  -ioap2(|)2(0) 


dz 


2 

dz 


<l>2(z) 


lzs=0 


r  p2  p2  1 
«  • 

fC^ 

f  Pi  Pi  'I 

icon2  icon2 

iconj  icoiii 

^  ^0  ^0  J 

^  ^0  ^0  J 

Definition  of  discontinuity  matrix: 


Til 

ni 


VP2  ^2 


Pl^  V 
Pz 

£l+5l 
p2  Hz; 
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Transfer  matrix  approach  for  frequency  shift: 


Refractive  index 


hi  L2  L3 


t(CO)  —  p3  •  d23  •  p2  •  di2  •  Pi 


Pressure  release  boundaries 


(1  l)t(co)n 


=  ti  1  (CO)  -  ti2  (CO)  +  t2i  (CO)  - 122  (CO) 


Rigid  boundaries 

(1  -l)t(co)n 


=  til  (CO)  +  ti2(C0)  -  t2i(C0)  -  t22(C0) 


309 


Normal  Mode  Frequency 
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Normal  Mode  Frequency 
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Length  of  Cell 


width  of  Layer 
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Conclusions  on  Frequency  Shifts: 

•  Even  for  volume  fractions  as  small  as  0.001 , 
radiation  pressure  Induced  particle  migration 
Introduces  frequency  shifts  Af  as  large  as  200  ppm. 

•  Adiabatic  invariance  model  based  on  Gorkov 
radiation  potential  well  yields  the  general  magnitude 
of  the  final  Af. 

•  The  approach  to  equilibrium  Af  following  pump  wave 
turn-on  has  been  modeled  with  moderate  success. 

•  Frequency  shift  of  a  simulated  resonator-layer 
system  computed  exactly  from  the  transfer  matrix 
agrees  will  with  adiabatic  invariance  model  for  low 
simulated  volume  fractions. 
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Resonances  and  other  mechanisms  for  elastic 
contributions  to  the  scattering  of  sound  by  objects  in  water 


Philip  L  Marston 
Scot  F.  Morse 
Karen  Gipson 

Department  of  Physics 
Washington  State  University 
Pullman,  WA  99164-2814 


Work  supported  by  the  Office  of  Naval  Research 


Assistance  in  measurement  for  large  shell: 

B.  T.  Hefner  (WSU) 

K.  L.  Williams  (Applied  Physics  Laboratory,  University  of 

Washington) 
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Outline 

1 .  Experiments  on  resonances  of  a  large  spherical 
shell  in  Lake  Union. 

2.  Normal  mode  description  of  resonances: 
Resonance  Scattering  Theory  and  fluid  loading 
effects  for  perfectly  elastic  spheres. 

3.  Ray  descriptions  of  resonances 

•  Leaky  or  supersonic  waves 

•  Subsonic  waves 

•  Anomalies  for  shells 

4.  Scattering  amplitude  for  high  Q  resonances  of 
spheres 

shapes 

5.  NoncanonlcaHHilB:  High  frequency  meridional 
ray  backscattering  enhancements  for  high 
frequency  sonar  (not  a  resonance). 
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RA^  HENDERSON 


PVDF  sheet 


source 


Sphere:  stainless  steel  316 

outside  radius  a  =  29.85  cm 

thickness  h  =  1 .52  cm 

h/a  =  5.1% 

dry  weight  130  kg 

density  r  =  7.98  g/cm^ 

longitudinal  sound  speed  Cl  =  5770  m/s 
shear  sound  speed  Cs  =  31 1 0  m/s 
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pC)ru)<i.ro 

measure 
OYi  YHrP 
Tr>od€. 
YieSoToahct- 
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“flexural”  resonances 


Frequency 


n 

calculated  (kHz) 

measured 

2 

1.46 

1.50 

3 

1.76 

1.79 

4 

2.05 

2.03 

5 

2.31 

2.31 

6 

2.60 

2.60 

7 

2.96 

2.96 

8 

3.40 

3.39 

9 

3.93 

3.93 

(kHz) 
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BACKSCATTERING  BY  ISOTROPIC  SPHERES 


plane  wavefronts 


Total  far-field  pressure  Y'  » 

p  =  e  '^  +  (a  e  '^r/  2r)  f 
f  =  form  function 

Backscattering 

fpws  =  (2/i(ka))  X  (-1)n  (2n  +  1)  Bn{ka)/Dn(ka) 

Bn  and  Dp  are  5  x  5  determinants 

(empty  shells) 

[or  3  X  3  determinants  (solid  spheres)] 

containing  spherical  Bessel,  Neumann,  and  Hankel  functions. 
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RST  analysis  of  monostatic  and  bistatic  acoustic  echoes  from  an 
elastic  sphere 

G.  C.  Gaunaurd  and  H.  Oberall*’ 

Naval  Surface  Weapons  Center.  R43  White  Oak,  Silver  Spring,  Maryland  20910 

J.  Acoust.  Soc.  Am.  73  (1 ).  January  1983 


FIG.  2.  Form-function  versus  k^a  for  the  tungsten  carbide  sphere.  (A)  Mea¬ 
surements  (from  Ref.  13),  (B)  classical  theory,  (C)  comparison. 
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Resonance  Scattering  Theory:  Flax  et  al..  JASA  63, 723(1 978). 

Flax,  Gaunaurd,  and  Oberall,  Physical  Acoustics,  Vol.  15(1981). 

Product  Expansion  of  the  S  function: 

Marston,  KargI,  and  Sun,  in  Acoustic  Resonance  Scattering, 
edited  by  H.  Uberall  (Gordon  and  Breach,  1992)  pp.  305-333. 


f(x)  =  Z  f„(x) .  =  r'  ,  X  =  ka  j  ) = 0 

Coyr)pte)i  ^ 

(-)  =  Infinite  Cylinder 


n=0 

(+)  =  Sphere, 


...-1 


fn(x)  =  (2i)  Sjj(x)  “  1]  gjj »  =  S  function  with  IS^I  =  1 


322 


TRANSPARENCY  32 


3,  SeptfatioaofihcfimitvmpgtUI  waves  (left)  Into  montnrrf|oemcr)  and  rigid  background!  (right|.  Resonances  are  labeled  by  an  index /with 
>><Kle  n.  Arrows  show  resonance  fcatnrca  (or  dips)  in  the 


A  Acoost  Soc.  Anu  VoL  73,  No.  1,  January  1fi83 


G.  C.  Gaunaurd  arKj  H.  OberalL*  Echoes  from  elastic  sphere 
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Backscattering  amplitudes  from  a  tungsten  carbide  sphere 


specular  form  function  Ifsl  =  (Pe^l  "  Pw^VCPe^l  +  Pw^) 


Gr  is  a  coupling  coefficient  which  accounts  for: 

(1 )  interaction  for  the  acoustic  field  with  the  Rayleigh 
wave  and  (2)  axial  caustic  for  sphere. _ 
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Aka. 


^j,yc 


ut 


-^'k- 


Fabr^j -Perot  t^r-  +Ke.  Xtts.  Sar^c^ 

—  e,)0,  iTji  oA A vA£.# 

f  —  ^  r 

•//  — TT”; - — — - - -  )  x^ka 

[14-  exp  (  —  iTTpi  4-  Httxc/ci  )  ] 


1.  U/or*I<s  ho'ti^  o«a”  a.nA  OSS’  ^‘fiSO^^o^.v^t^  , 

;z.  The  dar>xpCv»^  j^eed  not  be  ri^A.IC. 

s.  T^esononcfis need  not  be  -iSclA.'Ve<4. 

GTD  for  backscattering  from  elastic  spheres  and  cylinders 
in  water  and  the  coupling  of  surface  elastic  waves  with  the 
acousUcfleld  Philip L Marston  Tfl58  {Tartuary.nae 
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Williams  and  Marston  JASA(1986) 

^lid  Sphere  I  I  I  I  , 

RaySynthesis  ,  I  I 


t  I  \  / 1 


KargI  and  Marston  1989,  JASA(1991) 
Thick  Shell 
Ray  Synthesis 


L  £%act 

r  pcoi 


(3,  ,  C  /C  from 

^  comcLe>-  ^Y'€al 
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G/from  comparison  with  RST 

Ray  theory  near  a  resonance  reduces  to: 

f  -i' 

f,w  .  X ;  , 

where  x  =  Rn/  denotes  the  resonance  conditions 

Require  that  complex  frequency  plane  residue  match 
the  RST  residue  in 

(+)  (-) 
sphere:  n  =^xc/c^^ -(1/2);  cylinder:  n  =xc/c^  , 

gives  =  1 

g|^\x)  ~  SttP^  c/c^  ,  g|  \x)  *  87cP^(7Cx)^^^  . 

Marston  JASA  (1987, 1988).  For  rigid  background: 

r  (-)i  r  (+)i 

As  ka  oo,  arg  G^  ^  k/4,  arg  G^  -»  0 
General  background:  Marston,  Wave  Motion  (1995) 
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\f\ 


form  function  for  thin  SS304 
spherical  shell  with  h/a=0.025 


so  leaky  wave  interference 
with  specular 


s„  wave 
0 
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Marston  and  Sun  ,  "Resonance  and  interference 
scattering  near  the  coincidence  frequency  of  a  thin 
spherical  shell:  An  approximate  ray  synthesis,"  JASA 
92,  3315-3319(1992) 
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High  frequency  back  scattering  enhancement  for  general  shapes 

Experiments 

1 .  G.  Kaduchak  et  al.,  “Elastic  wave  contributions  in  high- 
resolution  acoustic  images  of  fluid-filled,  finite  cylindrical  shells 
in  water,"  JASA  100,  64-71  (1996). 

2.  S.  F.  Morse  et  al.,  "High  frequency  backscattering 
enhancements  by  thick  finite  cylindrical  shells  in  water  at 
oblique  incidence:  experiments,  interpretation  and 
calculations,”  (submitted  JASA). 

3.  K.  Gipson,  Experiments  in  progress. 

Theory 

1 .  P.  L.  Marston,  “Leaky  waves  on  weakly  curved  scatterers.  II. 
Convolution  formulation  for  two-dimensional  high-frequency 
scattering,"  JASA  97,  34-41  (1995). 

2.  Marston,  "Approximate  meridional  leaky  ray  amplitudes  for 
tilted  cylinders:  end-backscattering  enhancements  and 
comparisons  with  exact  theory  for  infinite  solid  cylinders” 
(accepted  in  JASA). 


332 


TRANSPARENCY  42 


Applied  Research  Laboratory 
The  University  of  Texas  at  Austin 


G.  Kaduchak  and  C.  M.  Loeffler,  "Eiastic  wave  contributions  in 
high  resolution  acoustic  images  of  fluid-filled,  finite  cylindrical  shells 
in  water,"  to  appear  in  JASA  June  159G.  CvH6}. 


High  resolution  image  of  finite  cylindrical  shell:  b/a=0.9449 


irrcidenceangle:  34.6'  deg. 


X  X 

X  X  X 

X  X  X  X 
XXX 
X  X 

X 


0  20  40  60  80  100 

distaece  -  centimeters 
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^T 


Spatial  surface  convolution  approximation  of  leaky  wave 
contributions  to  three-dimensional  scattering 


•  gives  spatial  dependence  of  surface  wave  amplitude  directly 
(for  image  interpretation) 

•  not  limited  to  thin-shell  mechanics 

•  avoids  singularities  of  "pure  ray  theory" 

•  meridional  ray  example  is  for  an  especially  large  contribution 
and  gives  expected  flat  surface  limit  as  radius  a  diverges. 

'P-LHantRM )  Ca^Asiyophc  ophcs 

ii^  f  PA/sJoU  A  equities 

M^Z) 
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Meridional  leaky  ray  amplitudes.  ^  i 

Example:  Ravieigh  wave  on  solid  stainless  stea  cylinder. 
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S$loH  Cy/inJn'cai  sktU  in  lOA-ht' 


b/a-  =  A/fl  =<5.or6  ka.s32.7 

C^.  Morst  WSU) 

No  a«ljiAst4W  pom<*4C'kfS  />t  »w*clai 
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Exact  and  ray  synthesized  total  scattering  cross 
section  of  a  thick  spherical  shell  in  water.  KargI  and 
Marston  JASA  (1990). 


Note:  (1 )  There  are  no  acoustically  active  high  Q  low 
frequency  modes;  (2)  The  torsional  modes  of  the  shell 
are  not  acoustically  active  (e.g.  T2  is  at  60  kHz  or  ka  of 
4.9  and  T3  is  at  94  kHz  or  ka  of  7.6). 
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Magnetic  excitation  and  acoustic  detection  of  spherical 
shell  modes  in  air  and  in  water 


Brian  T.  Hefner 
Philip  L  Marston 

Department  of  Physics 
Washington  State  University 
Pullman,  WA  99164-2814 


Work  supported  by  the  Office  of  Naval  Research 
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Ultrasonic  spectroscopy  of  metallic  spheres  using  electromagnetic- 

acoustic  traj^duction 

WardL  Johnson,  Stephen  J.  Norton,  Felix  Sendee,*^  andRobert  Pless 

Metallurgy  IXtiaai.  NattoriailriatiitetfStinulardxttrid  Tadmelegy,  Oaitkeaburt  MMylaiid20S99 


0) 


FIG.  1.  (a)  Sample/trxuisducer  configuration.  The  spherical  metallic  sam¬ 
ple  is  supported  between  the  poles  of  two  permanent  magnets  and  surround¬ 
ed  by  a  coil  with  its  axis  parallel  to  the  magnetic  field.  The  pedestal  support¬ 
ing  the  sample  and  insulating  tube  supporting  the  coil  are  not  shown,  (b) 
Schematic  of  transduction  physics. 


<*) 


(b) 


FIG.  2.  Same  as  Fig.  1,  except  the  coil  axis  is  perpendicular  to  the  magnetic 
field. 
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I- — 5  H” — 
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Total  field  8(7 ,t)  =  B,+  BgCt,!)  where  B,  =  B^i 

is  a  uniform  static  field.  What  is  the  distributton  of 
Maxwell  stresses  on  the  shell  in  a  B  field  with  an 
oscillating  component? 

Note:  At  typical  excitation  frequency  of  64  kHz 
the  skin  depth  5  for  aluminum  is  less  than  the 
shell  thickness  h  =  0.5  mm. 

1 .  The  oscillating  B  field  obey  approximately  the 
diffusion  equation  and  only  weakly  penefrates 
the  interior  of  the  sphere 

V^B  =  p^a^  ,  o = conductivity 


2.  As  a  first  approximation  the  shell  may  be 
treated  as  a  perfect  electrical  conductor  (PEC) 
so  that  02(7, t)  is  completely  excluded  from 
the  shell.  \ 
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3.  The  field  pattern  is  roughly  that  of  a 
static  B-fie!d  for  a  perfectly  diamagnetic  sphere 
(pjngjjjg=  0).  For  a  uniform  applied  field  Bq  sincot 
with  B  =  Bpi,  at  the  sphere’s  surface. 

B  (r  =  a,  t)  =  -^B^Csincot)  (3/2)  sine 


Total  B  =  zB  -  ^  B  sin0 

I  o 

B  =  B  (3/2)  sinwt 
3  0 
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Local  Maxwell  Stress  on  the  surface 


A 


0'^ 


— )■■)  4  — >  2  ^ 

I 


T  = 

Normal  (or  radial)  stress  component 


F- 


Mr 


f 


2  x2  ,  2 

where  B^  =  (B^  +  B^j  sin  0 


2X2 
Bj.  =  B^cos^e 
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-  ^  B^CsInef  (sIncDtf  (Te<^^) 

Term  1  describes  Lorentz  force  of  static  field  on 
the  eddy  current.  It  oscillates  at  a>. 


Term  2  describes  the  radial  Maxwell  stress  of  the 
oscillating  term  alone.  It  has  an  oscillating 

component  at  2to  since  (sincot)  =  ^  (1  -  cos2fflt). 

Both  Terms  1  and  2  couple  to  the  n  =  2  shell  mode: 
since  P2(cos6)  =  (3cos%  - 1)/2 , 
sin2e  =  1  -  cos2e  =  -  (^3)  [Pg  (cose)  - 1] 

SfZiieoS 


Mlc.  Signal  (10''-5) 
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Steel  Shell  (offset)  %d.th  Bias  Field 


30000 


Uoooo  SOOOO  60pd0^  70009  IBOCOO  90000  1  {l^OOO 
A  Fri^ency(Hz)  ^  ^ 


Wave  velocities  for  stainless  steel  shell 

Tabulated  440c  (m/s)  Fitted  (m/s) 
Cl  5854  5960 

Co  3150  3321 


Mode 


Calculated  (kHz)  Measured  (kHz) 

fitted  Cl,  Cs  I  tabulated  440c  c^,  c7” 
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Ruid  Loading  Experiment 
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Fluid  Loaded  Shell  (icqpcxise  at  2f) 


Wave  velocities  for  stainless  steel  shell 

Tabulated  440c  (m/s)  Fitted  (m/s) 
Cl  5854  5960 

C3  3150  3321 


Hydrophone  Signal  fl0*>8) 
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Non-magnetic  Huid  Loaded  Shell  | 


Non-magnetic  shell: 

•  No  response  at  f. 

•  No  resonances  excited  other  than  A2.  Most  notable 
is  the  absence  of  the  sharp  resonances  present  with 
the  magnetized  shell. 


TRANSPARENCY  64 


Fluid  Loaded  Shell  w/  m  perpendicular  to  the  axis  of  the  coil 


30000  40000  50000  60000  70000  80000  90000  100000 

Frequency  (Hz ) 


Fluid  Loaded  Shell  w/  m  parallel  to  the  axis  of  the  coil 
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Possible  torsional  excitation  mechanisms: 


1  Eddy  currents  and  “self  field”. 


«  «  Frequency  (MHz) 

FIG«  6.  Vibrstiosal  spectmm  of  the  3. 145^mn  2024  aiumxnum  sphere  with 
the  coO  configured  as  in  Fig.  2(a). 


“Self  field”  of  the  shell  is 
aligned  perpendicular  to  the  field 
of  the  coil. 


•  Drives'!^  mode 
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2.  Torque  on  distributed  magnetic  moment  by 
oscillating  fields. 

(a)  does  not  rely  on  eddy  currents 

(b)  for  uniform  magnetization  M  may  be  modeled 
by  force  on  surface  magnetization  charge 
density  (source  of  static  H) 


•  Uniform  M  drives  T3  mode 

Asymmetries  in  either  the  magnetization  or  the 
position  of  the  shell  in  the  coil  could  drive  the  T2 
mode. 
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3.  M  aligned  with  coil  axis:  Microscopic  domain  wall 
motion  or  domain  rotation  and  the  conservation  of 
angular  momentum.  (Einstein-de  Haas  effect  and/or 
torsional  magnetostriction) 


Suspension 


CD 


Suspended 

ferromagnetic 

specimen 


Magnetizing  coil 


Fig.  15.4.  Einstein-de 
Haas  gyromagnetic  experi¬ 
ment.  When  the  current 
in  the  magnetizing  coil  is 
reversed,  the  magnetization 
in  the  specimen  is  reversed 
and  the  specimen  rotates. 
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How  can  the  torsional  modes  be  radiating  sound? 


Two  possible  mechanisms: 

1.  The  sphere  is  attached  to 
fishing  line  by  a  small  bit  of 
epoxy.  This  could  be  moving 
with  the  shell  surface  and 
radiating  sound. 


For  example,  the  T2J2.  mode  will 
cause  the  epoxy  and  string  to 
move  parallel  to  the  axis  of  the 
coil  and  hence  radiate  sound. 


2.  The  shell  is  composed  of  two 
separate  halves  which  are 
bonded  together.  These  halves 
may  not  be  sufficiently  identical 
for  the  torsional  modes  to  be 
purely  rotational. 
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Conclusions: 

•  Demonstrated  that  the  modes  of  a  spherical 
metallic  shell  could  be  driven  magnetically  and  the 
acoustic  signal  could  be  used  to  measure  the 
response  in  both  air  and  water. 

•  A  static  bias  field  is  not  needed  if  the  acoustic 
detection  is  carried  out  at  twice  the  driving 
frequency. 

•  For  magnetized  spherical  shells  in  water  the  low 
lying  torsional  modes  could  be  excited  and 
detected  and  these  modes  have  very  large  Q. 
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MATERIAL  PROPERTY  MEASUREMENTS 
VIA  DIGITAL  GHz  INTERFEROMETRY 


Hartmut  Spetzler 


Geological  Sciences/CIRES 
University  of  Colorado 
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FIGURE  1 


Procedure  for  CT-EOS 

at  P=0  lo,  po,  Cp(J),  a(T) 
at  P  and  T  tj(T,P) 
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FIGURE  2 


A 


detector 


diamond 
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FIGURE  3 


FIGURE  4 
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Travel  Time  (ns) 


FIGURE  6 


Travel  Time  from  Raw  Interference  Amplitude  Data 
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FIGURE  7 
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FIGURE  8 


Travel  Time  Data 


Frequency  (GHz) 


FIGURE  9 


3.  OiMle^  A  A/noi  (82).  13%  Ca; 

/*  Sumtno  A  Niiftu«w«  (78).  2  2  %  Ca 
5  Sonczar  t  ai  (77).  2^.  Ca 
C;  Sumno  A  Nishuawra  PB).  0  6%  Ca 
7  Su^uto  A  Anoarsoft  (83).  0  6%  Ca 


FIGURE  10 
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Amplitude  Ratio  (Reflection  Coef.-brb/brn) 


FIGURE  12 
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Amplitude  Ratio  (Reflection  Coef.-brb/brn) 


FIGURE  13 
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FIGURE  14 


GHz  Interferometry  in  a  Diamond  Anvil  Cell 

(tinder  development,  CU  and  Bayreuth) 

a. )  Complete  Equation  of  State  over  P  (10  GPa)  and  T(I200  K) 

simultaneously. 

b. )  Small  samples  (50  microns),  thus  solid  solution  series  and  high 

pressure  phases  possible. 

c. )  Simultaneous  measurements  with  X-rays  yield  new  absolute 

pressure  scale. 
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FIGURE  16 
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Travel  lime  (nS) 


T 


FIGURE  17 
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FIGURE  18 


•  Experimental  data 

-  Least  squares  fitting  results 

----  calculated  using  Jackson 
and  Niesler  (1982)  results 
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APPLICATIONS  OF  ULTRASONIC 
INTERFEROMETRY 

Dipen  N.  Sinha,  Kendall  Springer,  Wei  Han,  and 

David  Lizon 

Los  Alamos  National  Laboratory 

Sponsors:  Office  of  the  Assistant  to  the  Secretary  of  Defense, 
Counterproliferation  Program 
Defense  Special  Weapons  Agency 
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Signal  Amplitude  Excitation  Amplitude 
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Determination  of  Sound  Speed  and  Attenuation 
from  Ultrasonic  Interferometry 
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impedance  ratio 
/  =  liquid 
w  =  wall 


d  =  path  length;  T  =  transmission  amplitude;  a  =  sound  attenuation;  f-  frequency 
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Ultrasonic  Interferometry  in  Multi-Layered  Systems 

Transducer 


Frequency  (MHz) 

•Detected  signal  is  a  superimposition  of  interference  in  both  wall  and  liquid 
•Peak-width  includes  effect  of  impedance  mismatch  at  the  wall-liquid  interface 
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Peak  Width  (Sf) 
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Non-Invasive  Liquid  Density  Determination 
Using  Ultrasonic  Interferometry 


Liquid  Density  Determination 


Non-invasive  Liquid  Density  Determination 
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Peak  Width  (KHz) 
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Container  Geometry  and  Material  Independence 


0  20  40  60  80  100 


Peak  Frequency  (f^  MHz^)) 
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Comparison  with  Theory 
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Sensitivity  of  Ultrasonic  Interferometry  Technique 

Single  Peak  Tracking 


Detection  of  Contamination  in  Liquids 


t^umber  of  drops 
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Ultrasonic  Interferometry  in  Cylindrical  Containers 


Tangential 

waves 


Dual-element 

transducer 


Compressional 

Wave 


Determination  of  Liquid  Viscosity  Effects 


Frequency  (MHz) 


Liquid  viscosity  can  be  determined  from  the  interaction  between  tangential 
waves  and  the  liquid  inside.  The  asymmetry  in  the  container  wall  thickness 
mode  resonance  peak  is  due  to  the  contribution  of  tangential  modes. 
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V 


AREAS  OF  APPLICATION  OF  ULTRASONIC 
INTERFEROMETRY 


Chemical  Weapons  Verification 
Basic  Research 
Industrial  Process  Monitoring 
Biomedical 
Food  and  Beverage 
Environmental  Sensors 
Materials  Diagnostics 
Customs  and  Drug  Interdiction 
Geological  studies 
Ultrasonic  Tomography 
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Non-Invasive  Identification  of  CW  Agents 


Artillery  shell 


Portab 


eSFAI 


Instrument 


SFAI  measurements  of  CW  agent 
physical  properties 
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Use  of  Ultrasonic 
Interferometry  in 
Non-Invasive 
Identification  of 
Chemical  Warfare 
Compounds  in 
Sealed  Munitions 


Portable  Instrument 


Chemical  Munitions 
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Amplitude  (mV) 
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Cylindrical  Piezoelectric  Resonator 


Frequency  (MHz) 


39 
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Noninvasive  Detection  of  Aerosol  in 
Sealed  Metal  Container 


Frequency  (MHz) 
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Ultrasonic  Interferometric  Characterization  of  a  Single 

Drop  of  Liquid 


Miniature  transducers  Frequency  (mhz) 
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Broad-Spectrum  Toxicity  Sensor  for  Monitoring 

Water  Pollutants 

Effect  on  liposomes  is  monitored  using  ultrasonic  interferometry 


Suspension  Liposome  Lipid  bilayer 
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! 

Ultrasonic  Interferometric  Characterization 
of  Bones  and  Joints 


Ultrasonic  Spectra  of  Finger  Joints 
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Effect  of  PCP  on  Liposome  Suspension 
Sound  Attenuation 


1  2  3  4  5  6  7  8 

Frequency  (MHz) 
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Particle  Size  and  Volume  Fraction  Determination 
in  Liquid  Suspensions 


Ti02- Water  Suspension:  f  =  2.22  MHz 


Volume  Fraction,  (ji 


Both  sound  speed  and 
attenuation  in  a  suspension 
depend  on  the  particle  size. 
Sound  speed  measurements 
are  shown  as  a  function  of 
volume  fraction. 

The  particle  size  can  be 
determined  by  curve-fitting 
the  measured  data  with  theory. 
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Measurement  of  Concentration  of  Chemicals 
using  Ultrasonic  Interferometry 


Sound  Speed  Vs.  Concentration  of  Sugar-Water  Solution 


Excess  Attenuation  Coefficient  vs.  Sugar  Concentration 


Frequency  (MHz) 


Sound  speed  and  sound  attenuation  as  a  function  of  frequency  can  be 
noninvasively  monitored. 
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Detection  of  Material  Deposition  Using 
Ultrasonic  Interferometry 
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Amplitude  (arb.  units) 
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Noninvasive  Characterization  of  Eggs  using 
Ultrasonic  Interferometry 
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TRANSPARENCY  23 
Not  Included  in  Transcript 


Non-Invasive  Monitoring  of  Petroleum  Products 


Gasoline  Ultrasonic  Properties 


Chemical  composition, 
of  petroleum  products  and 
other  chemicals  can  be  non- 
invasively  monitored 
using  ultrasonic  interfero¬ 
metry 
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RUS  use  in  chemical  weapons  convention  and  resonant 
photoacoustic  spectroscopy  of  aerosols 

by 

W.  Patrick  Arnott 
Desert  Research  Institute,  DRI 
Reno  NV  89506 

pat@sage.dri.edu 


RUS  use  in  chemical  weapons... 

1.  Manufacture  the  Los  Alamos  system  for  use  in  the  verification  portion 
of  the  Chemical  Weapons  Convention  (CWC). 

a.  Technology  transfer  supported  by  the  Defense  Special  Weapons 
Agency  and  DRI  cost  share. 

b.  Can  operate  within  DRI,  can  transfer  to  Desert  Research 
Corporation. 

2.  DRI  team  includes... 

a.  Rick  Purcell,  M.S.  Mechanical  Engineering. 

b.  Morien  Roberts,  Ph.D.  Computer  Engineering. 

c.  Dan  Wermers,  M.S.  Electrical  Engineering. 

d.  W.  Patrick  Arnott,  Ph.D.  Physics, 
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The  PTS  -  Preparatory  Commission  -  OPCW 

Home  Page 

The  Chemical  Weapons  Disarmament  Website 

Info  I  Documents  I  A  Guided  Tour  of  the  CWC  I  Contact  Address 
THIS  PAGE  LOOKS  NICER  IF  YOU  LET  IT'S  ORIGINAL  COLOUR  SCHEME  SHOW! 


Media!  NGOs! 

All  information  on  the  First  Conference  of  the 
States  Parties  (May  1997) 


Documentation  on 
a  common  transmission  file 
structure  is  now  available. 


A  list  of  interesting  jobs 
at  OPCW 


Do  not  forget  to  check  the 
Billboard! 

Most  recent  update  4  April  1997. 


The  forms  for  industrial 
declarations  are  now 
available 


At  this  website  the  Provisional  Technical  Secretariat  of  the  Pieparatoiy  Commission  for  the  Organisation  for 
the  Prohibition  of  Chemical  Weapons  (what  a  name!  just  call  us  the  Secretariat)  provides  information  on  the 
Chemical  Weapons  Convention  but  also  on  other  issues  relating  to  toxic  chemicals  and  chemical  protection. 

The  Chemical  Weapons  Convention 
(the  complete  text  of  the  Chemical  Weapons  Convention) 

tThe  PrepCom  and  its  members 

les  lists  of  points  of  contact  in  most  Member  States) 

What's  happenine  at  the  PrepCom 
(be  prepared  through  the  PrepCom  Calendarium) 

National  Authority  Adviser 

(almost  everything  a  new  or  well  established  National  Authority  needs  to  know) 

Fact-finding  files 

(plenty  on  chemical  weapons  and  other  toxic  chemicals,  protection,  relevant  publications,  etc.) 
to  other  Internet  sources 

(the  hottest  links  on  the  Web  for  anyone  interested  in  CW  matters) 

k  'th/>  .‘/'vtvtV.  OpC  \jJ .  ^J/ 
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Chemical  Weapons  Stored  in  the  US 
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Table  1 


Toxicity  of  Chemical 
and  Biological  Agents 

Estimated  Lethal  Dose 
in  mg/person*^ _ 


1000  rag 

Chemical 

100  mg 

Agents 

10  mg 

—  1  mg  — 1 

.1  mg 

.01  mg  -  Toxin  Agents 

- .001  mg 

.0001  mg  „ 

Biological  .00001  mg 

Agents  “  .000001  mg 

.0000001  mg 
_ .00000001  mg 

*  One  paperclip  weighs  approximately  500  mg. 

Source:  Office  of  Technology  Assessment^  Technolopes  Underlying 
Weapons  of  Mass  Destruction  (Washington,  D.C.:  U.S.  Govemment 
Printing  Office,  December  1993),  p.  77. 
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ARS  Papers 


f  Paper:  124K). 


('Presentation:  1.5M~)  R.  Roberts,  P.  Lewis,  and  O.  Vela.  A  pattern 


recognition  algorithm  for  the  blind  discrimination  of  liquid  and  solid  filled 
munitions.  (LAUR-95-3789)  To  appear  in  Proc.  29th  Asilomar  Conf.  on  Signals,  Systems,  and 

Computers,  IEEE,  Pacific  Grove,  CA,  Nov.  1995.  (Abstract) 


(Paner:  30K)  O.  Vela.  P.  Lewis.  R.  Roberts,  J.  Chen,  and  D.  Sinha.  Acoustic  classiHcation 

^  _  ...  y-m-  t.  -r  -r-w^  ■*  TT  T-k  m  W1 W  1  ^  I  T  >  ^  /»  . 


of  chemical  weapon  munitions.  (LAUR-95-1073)  In  Proc.  ONSITE  3rd  Int.  Conf.  On-Site 
Analysis,  Houston,  TX,  Jan.  1995.  (Abstract) 


(Paper:  154K). 


(Presentation:  242K)  R.  Roberts,  P.  Lewis,  J.  Chen,  and  O.  Vela. 


Techniques  for  classifying  acoustic  resonant  spectra.  (LAUR-94-3977)  In  Proc.  28th 
Asilomar  Conf.  on  Signals,  Systems,  and  Computers,  pp  1195-1199,  IEEE,  Pacific  Grove,  CA, 
Nov.  1994.  (Abstract) 


•  ii<.  l(Paperr  1 69K)  R.  Roberts,  J.  Chen,  O.  Vela,  and  P.  Lewis.  Munitions  classification  using 
an  acoustic  resonance  spectroscopic  technique.  (LAlJR-93-4108)  In  Proc.  27th  Asilomar 
Conf  on  Signals,  Systems,  and  Computers,  pp  991-995,  IEEE,  Pacific  Grove,  CA,  Nov.  1993. 
(Abstract) 
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(Presentation:  591K)  P.  Lewis.  Acoustic  Resonance  Spectroscopy  (ARS)  for 
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ARS  HOLOGRAPHY 

Video  holograms  showing  some  of  the  vibrational  resonance  modes  of  a  105mm  munition.  Each  "contour"  corresponds  to  1/2  micron. 


This  Page:  http:/Amw.esa.lanl9oy/ars/ars-home.html  -  Last  Update:  18  Oct  95 
Randy  Roberts  -  Los  Alamos  National  Laboratory  -  Group  ESA-MT 
Email:  roberts@lanlQov 
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SIMPLE  BLOCK  DIAGRAM 


1.  Digital  Signal  Analyzer  (DSA). 

2.  Acoustic  Resonance  Spectroscopy  Software  (ARS). 

3.  Piezoelectric  source  and  receiver  attached  to  shell  by  magnets. 
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ARS  IS  PORTABLE 

(PC  BASED  PEN  COMPUTER,  WINDOW  3.1) 


SENSOR  ATIACHMENT  TO  MUNITION 
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1 .  Set  up  the  computer 
for  your  particular  site 
and  expected  munitions. 


Acoustic  Resonance  Spectroscopy 


I  Site  Setup  ] 
I  Collect  Qata  j 
I  greate  Template  ; 

!  yiew&  Analyze  i 
i  Exit .  I 


2.  Make  templates  using 
munitions  with  known 
chemical  fills.  ARS  uses 
these  templates  (one  for 
each  fill  agent)  to 
estimate  the  chemical  fill 
in  other  munitions  of  the 
same  size. 


3.  Use  ARS  to  estimate 
the  type  of  fill  contents 
in  other  munitions. 


'0  !?!  '9i 

1^ 

AAAAAA 


:?  ;?j  I?  ;?!  i? 


:?!  !? 
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•  A/e<x//4 


Table  3.  Emission  factors  for  different  sources  of  ekmentai 
(black)  carbon  (from  Ogren  and  Charlson.  1984) 

Source 

Ekmentai  catfaoB  Ig) 

FndlUl 

Vehicles,  diesel  engine 

2 

Fireplace,  softwood 

1.3 

Jet  engine 

1 

Fireplace,  hardwood 

0.39 

Vehkles.  gasoline  engine 

0.02  , 

Solid  fossil  fuel 

0.001 

0.0003  '* 

Natural  gas 

Rqjrcxluced  by  pennistion  of  Miinskgaard  InternatioaaL 


u 


pu 

Ockvhr 


^Micr0/^AM 


/^c  loA^t.'^^r 


^3Z  /><v\ 
i^Lc4x><tccu^i  ^ 
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PRELIMINARY  EVALUATION  OF  THE  DRI  PHOTOACOUSTIC 

SPECTROMETER 

Comparison  with  aethaloraeter  data  from  the  Brighton  site  during  NFRAQS 

Dec  96  thru  Jan  97 

by 

W.  Patrick  Arnott,  Hans  Moosmiiller,  C.  Fred  Rogers 
Desert  Research  Institute 

Assistance  in  the  field  by 
John  Walker,  Pueblo  CO 
Mitch  Walker,  ARS,  Fort  Collins  CO 

532  nm  laser  provided  by 
Professor  Jin,  Professor  Briich 
Physics  Department,  University  of  Nevada  Reno 

Photoacoustic  spectrometer  development  supported  by 
Environmental  Protection  Agency  Office  of  Exploratory  Research 

Field  operation  supported  by 
Desert  Research  Institute,  IPA  grant 

OUTLINE 

1.  The  North  Front  Range  Air  Quality  Study. 

2.  The  photoacoustic  spectrometer. 

3.  Comparison  of  photoacoustic  spectrometer  and  aethalometer  results.. 
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The  North  Front  Range  Air  Quality  Study. 

The  Northern  Front  Range  Air  Quality  Study  (NFRAQS)  is  a  major  effort 

•  to  determine  the  sources  of  existing  air  pollution  in  the  Denver  urban 
region,  attributing  to  each  source  or  source  category  an  estimate  of  its 
emissions. 

•  to  collect  data  necessary  to  support  informed  decisions  leading  to 
attainment  of  state  visibility  gods  and  federal  air  quality  standards. 

DRI  is  the  main  contractor  for  NFRAQ 


OUR  FIELD  SITE  LOCATION: 
BRIGHTON  COLORADO 
(DEG  96,  JAN-FEB  97) 


X 

BOULDER,  CO 


X  BRIGHTON,  CO 

/ 


X 

DEIWER,  CO 
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_ Schematic  view  of  the  photoacoustic  spectrometer. 

Light  Absorption  in  Dimensions  of  Inverse  Distance  =  B 

B  =  -^res  ^  ^0 
^abs  y-1  Q  ’ 

fo  =  Resonance  Frequency. 

Q  =  Resonator  Quality  Factor. 

Pm  =  Peak  Acoustic  Pressure  at  f©. 
y  =  Ratio  of  Isobaric  and  Isochoric  Specific  Heats  For  Air. 

Pl  =  Peak  Laser  Beam  Power  at  fo. 

Q  =  resonator  quality  factor. 

Ares  =  Resonator  Cross  Sectional  Area. 

TYPICAL  VALUES: 

Pl  =  50  mW  to  105  mW  (depending  on  the  laser  used), 

Ares  =  5.07  cm^, 
fo  =  500  Hz, 

Q  =  80, 

Bandwidth  of  spectrometer  =  fo/Q  =  6.25  Hz, 

Pm  =  6  |iPa  (Pm  =  -10.5  dB  Rel.  20  jiPa), _ 

_ Babs  ~  7.5  Mm'^ _ 

Broadband  noise  background  equivalent  to  Babs  ==  0.5  Mm'^ 
ENB  W  of  Phase  Sensitive  Detection  =  7.5  mHz. 
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Block  Diagram  of  the  Photoacoustic  Spectrometer  and  Detection 

Electronics 


Lasers  used  were  both  very  compact  and  efficient: 

•  532  nm  frequency  doubled  NdYAG  laser  diode  pumped  laser. 


=  60  mW  modulated  power. 

•  685  nm  laser  diode,  «  102  mW  modulated  power. 
532  nm  laser  was  only  available  for  2  days. 


The  685  nm  laser  diode  was  confined  to  pass  cleanly  through  the 
spectrometer  only  with  the  use  of  an  initial  collimating  objective 
very  near  the  diode  followed  by  a  cylindrical  lens  to  reduce 
astigmatism  in  the  beam,  and  finally  followed  by  a  converging 
lens  to  reduce  the  beam  width. 
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Photoacoustic  Spectrometer  with  532  nm  Laser 


a) 


NFRAQS,  Brighton  CO  Site,  18  Dec  96 


Photoacoustic  Spectrometer  with  685  nm  Laser 


00:00:00  06:00:00  12:00:00  18:00:00  00:00:00 
NFRAQS,  Brighton  CO  Site,  7  Jan  97 


b) 


Typical  comparison  of  aetfaalometer  and  photoacoustic  spectrometer  determinations  of 
light  absorption  by  aerosols  during  NFRAQS.  In  a),  the  532  nm  laser  was  used  for  the 
photoacoustic  spectrometer,  though  the  685  nm  laser  was  used  in  b). _ 
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Crude  Theoretical  Investigation. 


Light  absorbing  aerosols  are  not  likely  simple  spherical  particles. 


Rayleigh  Scattering  Limit  of  Mie  Theory  for  Spheres: 

Specific  Absorption  Coefficient  ka  (m^/g): 

J,  ^_67C _ 

^  Pp^  (nr2-nj2+2)2+(2ni.nj)2  ’ 

where: 

Pp  is  the  particle  density,  X  is  the  wavelength,  and  nr  and  nr  are  the  wavelength 
dependent  real  and  imaginary  parts  of  the  particle  refractive  index. 


•  Small  particle  limit,  ka  is  independent  of  particle  diameter,  and  is  dependent  on 
wavelength,  both  explicitely  and  through  the  complex  refractive  index. 


MIE  ABSORPTION 


0.01  0.1  1  10 

Sphere  Diameter  (pm) 

Specific  absorption  coefficient  for  spherical  light  absorbing  particles. 
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VARIOUS  EFFORTS  WITH  ‘CYLINDRICAL’ 

RESONATORS 


DRI 

Plane  wave  mode,  =  500  Hz,  Q  =  80. 

(Laser  power,  integration  time  limited) 

532  nm  and  685  rnn  Semiconductor,  Solid  State  Lasers  50  mW  to  100  mW 
Detection  Limit  Bats  =0.4  Mm'^ 

OTHERS 

ADAMS: 

Best  result  with  azimuthal  mode,  5)  =  3970  Hz,  Q  =  220. 

(Window  absorption  limited) 

514.5  nm  Argon  Ion  laser  1  W  (horrible  eflSciency). 

Detection  Limit  Bats  =4.7  Mm'^ 

EUROPEAN  EFFORT  (Petzold  and  Niessner.  1996) 

Azimuthal  mode,  ^  =  6670  Hz,  Q  =  300. 

(Somewhat  window  absorption  limited,  flow  noise  problematic,  higher  Q 
resonator  made  for  lower  integration  time) 

802  nm  laser  diode,  450  mW 
Detection  Limit  Bats  =1.5  Mm^ 


References. 
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“Real-Time,  In  Situ  Measurement  of  Atmospheric  Optical  Absorption 
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Roessler,  D.  M.  (1984).  “Photoacoustic  Insights  on  Diesel  Exhaust  Particles.” 
Appl.  Opt.  23,  1148-1155. 
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Resonant  Photoacoustic  Spectroscopy  of  Optical  Materials 


J.  D.  Maynard,  Chang  Yu,  and  Wei-ll  Lin 
The  Pennsylvania  State  University 


Technology  Issue 

Develop  and  characterize  optical  materials  which  are  highly  transparent 

Applications 

•  Optical  fibers  in  long  distance  communications  systems 

•  Lenses  and  windows  in  high  power  laser  systems 

•  Electro-optic,  magneto-optic,  and  acousto-optic  components 


Objective 

Develop  a  fast  and  sensitive  method  for  measuring  very  small  optical  absorp¬ 
tion  coefficients 
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Measurement  of  Small  Optical  Absorption 


a 

(cm  "b 

Power  loss: 

Transmission 

1 

O 

Temperature  rise: 

Calorimetric 

10 

Re-radiation: 

Emittance 

10  -5 

Strain,  static: 

Interferometric 

10 

Strain,  acoustic: 

Photoacoustic 

10 

Resonant  Photoacoustic 

10 
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Pulse  photoacoustic 


Laser  energy  per  pulse  10“^  J 

Limited  by:  Noise  in  transducer 

Light  scattered  to  attached 
transducer 

Rep  rate  <  10^  Hz 


Resonant  photoacoustic 


CW  laser,  modulated  @  lO^Hz  10“^  J 

Modulate  at  acoustic  resonance  of  sample 

Gain  =  Quality  factor  (Q)  =  -  10^ 

Non-contact  transducers;  no  scattered 
light,  higher  Q 
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Laser  Beam  as  an  Acoustic  Line  Source 


Energy  absorbed 
Temperature  increase 


E  =  Wral 


AT  — 


Wral 

CppV 


Thermal  expansion 

Geometry 


AV  =  VfipAT 

=  TT  (i?  +  AR)^  I  -  ttR^I  =  2-kRARI 


Peak-to-peak  displacement 


Wra^p 

2-KRpCp 


Wppa 

Particle  velocity  amplitude  A  =  ujARf2  —  2T^RpCp 
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Pressure  radiated  by  pulsating  cylinder 
Green’s  function  for  line  source 
Equivalent  source  strength,  small  kR 

Displacement  source  strength 


P(r)  = 


ipcA 


.K  (kR) 

G  (r)  =  IttHo  (kr) 


Ho (kr) 


Sp  = 


pcA 


[7:H',ikR)\ 


c:i  —kRpcA 


S  = 


iL 

pCUJ 

kR  . 

:  - A 

2u 

kR  W(3pa 
2uj  2t^ RpCp 


kWfSpa 
A'KupC p 
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Green’s  function  for  a  finite  rectangle 


Add  source  strength  S,  evaluate  with  cos  =  1,  ka  —  mTr,  n  —  0 
Displacement  amplitude 
Optical  absorption  coefficient 


/  kWfipa  \  16  Q 
\4xa;pC'p/  7ra6(fc/7r)^ 
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vacuum  box 
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response  (mV) 
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frequency  (Hz) 
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Resonant  Bar  Acoustic  Spectroscopy: 

Characterization  of  viscoelastic  materials  using  the 
resonant  modes  of  a  free-free  bar 

David  A.  Browni  and  Qiushuang  Guo2 

1  Assistant  Professor  Department  of  Electrical  Engineering  and 
Center  for  Marine  Science  and  Technology 
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University  of  Massachusetts 
Dartmouth 


Outline 

Resonant  Bar  Acoustic  Spectroscopy 

0.  Introduction  and  Motivation 

1 .  Excitation  and  Detection  of  Modes  in  a  Free-Free  Bar 
Resonance  in  bars 

and  determination  of  elastic  moduli 

2.  Damping 

3.  Instrumentation 

Transfer  Function  Techniques 
Phase  Locked  Loop  Techniques 

4.  Some  More  Data 

(temperature  and  frequency  dependence  of  elastic  moduli) 

5.  Correction  Factors 

Warping  Functions 
Finite  Bar  length 
Transducer  Mass  and  Stif&iess 
Including  Dissipation  in  the  Wave  Equations 


D.A.  Brown,  May  1997 
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Resonant  acoustic  determination  of  elastic  moduli 

Steven  L  Garrett 

Physics  Department.  Code  61  Gx.  Naval  Postgraduate  School  Monterey.  California  93943 

( Received  1 1  December  1989;  accepted  for  publication  23  February  1990) 

A  long  rod  of  circular  or  elliptical  cross  section  can  be  selectively  excited  in  the  torsional, 
longitudinal,  and  flexural  resonant  modes  using  a  single  pair  of  clectrodynamic  transducers 
consisting  of  coils  glued  to  each  end  of  the  bar  and  placed  in  the  field  of  a  magnet.  The  clastic 
moduli  of  the  material  can  be  simply  determined  to  high  accuracy  by  measuring  the  resonance 
frequencies  of  these  modes  and  the  mass  and  physical  dimensions  of  the  rod.  Smee  the 
longitudinal  and  flexural  modes  both  yield  values  for  the  Young’s  modulus,  the  system  has  a 
built-in  redundancy  that  makes  the  interpretation  of  the  data  more  robust.  The  large  signal 
amplitudes  make  automated  measurement  of  the  temperature  dependence  of  the  moduli  simple 
through  the  use  of  a  phase-locked  loop  to  track  the  change  in  the  resonance  frequency  with 
temperature. 

PACS  numbers:  43.40.Cw,  43.20. Ye.  62-20.Dc.  07.10  -h  i 


INTRODUCTION 

The  accurate  measurement  of  the  clastic  constants  of 
materials  and  their  dependence  on  temperature,  static  pres¬ 
sure,  and  other  ambient  parameters  is  important  in  many 
fields  of  science  and  engineering  research,  as  well  as  in  prod¬ 
uct  design  and  quality  control.  This  is  panicularly  true  for 
hydrophones,  since  most  transduction  schemes  involve  the 
measurement  of  the  deformation  of  some  solid  material  in 
response  to  a  change  in  pressure.  The  elastic  modulus  relates 
these  strains  to  the  applied  stresses,  so  its  value  and  its  tem¬ 
perature  dependence  are  important  design  parameters.  In¬ 
terest  in  the  measurement  of  these  moduli  in  the  acoustical 
transduction  community,  particularly  for  castable  poly¬ 
mers,  has  increased  recently  due  to  the  development  of  inter¬ 
ferometric  fiber-optic  hydrophones,*"^  which  measure  the 
strain  in  the  optical  fiber  induced  by  pressure  changes  in  the 
fluid.  Castable  polymers  are  attractive  for  several  hydro¬ 
phone  designs  that  use  shell  structures^’  or  encapsula- 
tion**"*  *  for  strain  enhancement,  since  the  optical  fibers  can 
be  cast  directly  into  those  materials  as  well  as  for  more  con¬ 
ventional  piezoelectric  hydrophones.'^*^ 

Manufacturer’s  specifications  for  clastic  constants  of 
castable  polymers  are  not  particularly  useful  for  the  hydro¬ 
phone  designer  since  they  arc  usually  determined  by  static 
techniques,  rarely  contain  more  than  one  modulus  (two  arc 
the  minimum  required  to  uniquely  specify  the  elastic  re¬ 
sponse  of  an  isotropic  material*^),  vary  widely  depending 
upon  sample  preparation  (c.g.,  catalyst)  and  cure  tempera¬ 
ture,  and  never  contain  information  about  the  temperature 
dependence  of  the  moduli.  It  is  also  important  in  hydro¬ 
phone  design  to  measure  the  clastic  moduli  at  the  frequen¬ 
cies  of  intended  operation,  since  it  is  well  known  that  the 
static  and  dynamic  moduli  of  plastics  can  differ  substantial¬ 
ly*^  due  to  the  existence  of  relaxation  time  effects  that  can 
make  the  modulus  measured  by  conventional  quasi-static 
stress-strain  curves  significantly  lower  than  the  dynaniic 
modulus.  Since  it  is  the  dynamic  modulus  that  determines 
the  acoustic  sensitivity,  and  the  static  modulus  that  deter¬ 
mines  the  deformation  due  to  increased  operating  depth,  it  is 
important  that  both  are  known  to  the  hydrophone  designer. 

A  recent  scries  of  measurements  on  3 1  elastomeric  sam- 

210  J  Acoust  Soc  Am.  88  (1).  July  1990 


pics  by  Lagakos  et  al.  *^*’  illustrate  how  the  assumption  of  a 
simple  form  for  the  frequency  dependence  of  the  Young’s 
modulus  can  lead  to  unacceptable  disagreement  between  the 
measured  moduli  and  the  predicted  value.  For  the  case  of 
Uralite  3130.  a  polyurethane  used  in  a  prototype  fiber-optic 
planar  flexible  hydrophone  design,’  the  variation  between 
the  model  and  measured' modulus  was  in  excess  of  750%. 
One  must  question  the  validity  of  their  temperature  coeffi¬ 
cient  of  the  modulus  measured  at  1  MHz  when  extrapolated 
to  operating  frequencies  of  less  than  a  kilohertz. 

Because  of  the  number  of  variations  in  the  preparation 
that  arc  available  to  the  experimentalist  (cure,  catalyst,  fill¬ 
ing  fraction  for  composites,  etc.),  it  is  important  to  have  a 
convenient  technique  for  measuring  at  least  two  moduli  for 
reasonably  small  samples.  Since  the  derivation  of  a  third 
modulus  from  the  measurement  of  two  others  can  be  very 
sensitive  to  the  values  of  those  moduli  (particularly  for  sam¬ 
ples  with  Poisson’s  ratio  close  to  0.5 ),  it  is  important  that  the 
measurement  technique  have  very  high  precision.  In  the  fol¬ 
lowing  sections,  a  technique  is  described  that  is  convenient, 
accurate,  precise,  and  economical,  that  depends  upon  the 
measurement  of  the  frequencies  of  the  longitudinal,  flexural, 
and  torsional  resonant  modes  of  a  single  rod-shaped  sample 
of  circular  cross  section  using  the  same  two  transducers  to 
excite  and  detect  all  three  modes. 

The  fact  that  the  technique  is  resonant  insures  high  sig- 
nal-to-noise  ratio  while  the  fundamental  measurement  being 
a  frequency  means  that  one  can  obtain  extremely  high  preci¬ 
sion  with  an  inexpensive  instrument  (i.c.,  a  frequency 
counter).  The  technique  can  be  used  with  both  insulating  or 
conducting  samples  that  are  not  ferromagnetic.  An  addi¬ 
tional  attractive  feature  of  this  technique  is  the  fact  that 
Young’s  modulus  can  be  determined  independently  in  two 
ways  by  the  measuring  the  frequencies  of  the  longitudinal 
and  flexural  modes.  This  provides  the  experimentalist  with 
immediate  feedback  on  the  self-consistency  of  his/her  data 
and  analysis.  (A  man  with  one  watch  thinks  he  knows  the 
time;  a  man  with  two  is  never  sure!) 

The  technique  for  measurement  of  the  torsional  mode  is 
a  refinement  of  one  developed  first  by  Barone  and  Giaco- 
mini'**  to  study  the  modes  of  vibration  of  bars  having  van- 
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Garrett,  S.L.,  “Resonant  acoustic  determination  of  elastic  moduli,” 
J.  Acoust.  Soc.  Am.  88(1),  21 0-221  (July  1990) 
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Resonant  frequencies  of  free-free  bar 

If  the  clamping  in  the  bar  is  not  considered,  the  solution  can  be  expressed  as  the 
linear  combination  of  the  “mode”,  plane  harmonic  wave  solutions  whose  frequencies 
are  the  integer  multiple  of  a  basic  frequency  for  nondispersive  modes  (i.e.  longitudinal 
and  torsional  modes)  and  nearly  integer  multiple  of  a  basic  frequency  for  dispersive 
mode  (flexural). 


Under  the  condition  of  free-free  ends,  for  nondispersive  modes, 

=  ncL  /  2L;  n  =  1,2,3,-  • 

f J  =  ncj  /  2L;  n  =  1,2,3,-  -  -. 

where  Cl  =  Ve7p,  Ct  =  ^/t7p, 

For  fiexural  modes,  the  frequencies  are 

=  tih^ClK  /  8L^ ;  n  =  3.01,5.00,7,9, 1 1 


where  l  is  the  length  of  the  bar, 


Cj,  and  Ci  are  the  phase  speed  of  longitudinal 


K  is  the  radius  of  the  gyration 


and  torsional  waves, 
dxdy 


E 


F 


is  the  Young’s  modulus  —  — 

A 


L 


G  is  the  shear  modulus 
P  is  the  density  of  the  bar. 
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Elastic  Moduli  in  terms  of  Resonant  Frequencies 

So  it  is  a  rather  simple  matter  to  determine  the  elastic  moduli  from 
the  resonance  frequencies  of  a  free-free  bar 


G  =  4pL2(fJ/n)^ 


from  torsional  mode, 


E  =  4pL2(f„^/n)2 


from  longitudinal  mode, 


2  2 
71  K 


from  the  flexural  mode. 
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Complex  Modulus  and  Q 

Equation  of  motion  for  a  damped  driven  oscillator 

M  — +  R  — +Kx  =  F(t) 

If  the  applied  force  is  sinusoidal  of  frequency  ©, 

-co2Mx  +  K|l+j^Jx  =  F* 

where  the  coefficient  of  the  effective  stiffness  term  can  be 
expressed  as 

K*  =  K  ^  ^ 

At  resonance,  o)=cdo  =  Yk/m,  the  loss  tangent  is  inversely  related  to 
the  quality  factor,  Q  by 

1  _  K  _YkM  ^cooM_q 

tan  6  ^  ^  ^ 


445 


TRANSPARENCY  11 


Transfer  Function  Method 


u  =  velocity  V  =  voltage  B  =  Magnetic  field 
F=  force  I  =  current  L  =  length 
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Torsional  Mode 

Temp  =  22.4  C  ,  ASA  Sample  Y 


FREQ  RESP 


17 .7601k 
-10.9356  ±J 
-43.072  ±J 
-93 . 6262  ±J 
-ISO . 067  ±J 
410.367  ±J 
-223.471 
-269.036  ±J 


99 . 1363 
525 .016 
1 . 1065k 
1 . 66311k 
2. 172k 
2.32244k 
2.94455k 


-297 . 62 
-41.933 
9. 1555 Ik 
757 . 679  ±J 
-636 . 194  ±J 
651.556 
351.79  ±J 
-732 . 629  ±J 
524.93  ±J 


S45 . 197 
1 . 19935k 
1 . 67S7k 
2.20706k 


2 . 42646k 


(MPa) 
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Data  Analysis  Torsional  Mode 

Mode  #  Frequency  Freq/Mode  #  Quality  factor 


[n] 

[fn]  Hz 

[fn/n] 

[Q] 

1 

525 

525 

6.12 

2 

1106 

553 

5.93 

3 

1683 

561 

5.59 

4 

2322 

580 

5.22 

5 

2944 

588 

5.12 

G  =  4pL2 1 

ft  =  shape  correction  factor,  1  for  a  circular  rod 
ftn  =  transducer  mass  corection 

G  =  4(l025)(0.2733f  (525f(l.20Xl.23)=  124  M  PA 
525  Hz;  23  C 


-G- 


#  toss  tang  % 


0  1000  2000  3000 

Fioqwoncy  (Hz) 
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Phase  Locked  Loop  Method 
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Ilhatco  1201 

prc'Uiq* 

Buid  P&u 
Filter 


Block  diagram  of  the  instnimentation  used  to  track  the  dynamic  modi 
dependence  on  temperature  and  frequency. 
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Shear  modulus  of  PR1592  as  a  function  of  frequency  and  temperature. 


SiearJkbeMtf/Pj/ 
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Transfer  function 


Thus  the  transfer  admittance  function  has  poles  at  y  and  7*  and  a  single  zero  at  (d — 0. 


Y=a+jb 


Pole  zero  plot  of  coiiq)lex  admittance  function. 


If  we  express  the  location  of  the  poles  as  ,  then  Q  can  be  expressed  as  follows 


Q= 


?  +  a 
-2a 


and  the  resonance  frequency  is  equal  to 
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SHEAR  MODULUS  (Pa) 


wm 


CVI  CO  to 

+  +  +  + 

III  tu  III  IIJ 


EQUENCY  (Hz) 


LOSS  TANGENT 
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Complications  associated  with  soft  materials 
Complications  high  loss  materials 


Correction  factors  due  to  added  mass  and  stiffness  of 
Transducers  becomes  very  Important. 

Materials  can  be  so  soft  that  they  can’t  support  a  constant 
shape 

Losses  can  be  so  high  that  there  Is  no  apparent  resonance 
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Corrections  Factors 

The  simple  (thin  bar)  theory  is  often 
not  sufficiently  accurate 

Corrections: 

due  to  coupled  motion  in  a  bar 
(Saint  Venant  warping  function). 


due  to  thick  or  short  bar 


due  to  added  transducer  mass 


due  to  added  stiffness  of  the  transducers 
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Coupled  motion  in  a  bar  correction  factor 
(Saint  Venant  warping  function). 

When  the  bar  is  excited  to  a  particular  mode,  any  element  in 
the  bar  moves  mainly  due  to  the  particular  mode.  But  still  it  has 
to  move  along  the  perpendicular  direction  because  the  volume  is 
constant. 


Longitudinal  Mode:  an  element  that  is  not  at  axis  of  the  bar 
moves  both  along  the  axis  and  in  the  radial  direction. 


Torsional  Mode:  an  element  move  both  In  the  cross  section 
plane  and  in  the  longitudinal  direction. 


The  coupled  motion  and  hence  correction  factor  can  be 
substantial  as  the  elastic  modulus  or  correctional  area  is 
Increased. 

[Ref  Love,  A  treatise  on  the  mathematical  theory  of  elasticity]. 
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Coupled  motion  case  correction  factors 


When  the  coupled  motion  among  different  dimensions  is 
considered,  the  resonant  angular  frequencies  (o>)  can  be 
recalculated  by  the  energy  methods.  The  angular  frequency 
change  is  A(b.  And  the  angular  frequency  of  the  “bare"  bar  is 
M  is  the  mass  of  unloaded  bar.  Am  is  the  added  mass(include 
coils  and  glue),  l  is  the  length  of  the  bar.  a/  is  the  length  of  the 
transducer  along  the  axis  of  the  bar.  The  added  mass  Is  regarded 

as  a  point  mass  at  the  position  away  from  the  end  by  y.  For 


torsional  mode, 


G  =  4pU 


n 


Qs^l 


where  Qt  =  i  -  ■  ^  - —  i 

'  MKp2  I  2L 

which  is  slightly  smaller  than  the  above  correction. 
(Danielson) 


,  AmK^  2f  nJiAl 

1  + - ^cos 


n2 


J 
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Torsional  Examples  for  different  cross-sections 


For  a  cylindrical  symmetric  bar, 
the  warping  function  0=0, 

For  an  ellipsoidal  cross-section  bar 

with  semi-major  a  (along  ;c  axis)  and  semi-minor  b. 


a^-b 


For  a  rectangular  cross  section  bar  of  sides  20  (along  x  axis)and 
2b,  (Ref.  Spinner  and  left) 

.  ,  (2n+l)7cx 

3  ^(2n+l)7ta  2b 

cosh - 

2b 


0=-xy+4b' 


(-1) 


n 


n=0(2n+l) 


459 


TRANSPARENCY  25 


Correction  for  shear  modulus 


For  torsional  mode,  the  corrected  formula  of  shear  modulus 
becomes: 


G  =  4pU 


n 


m 


2x2 


^  _  (a"  +  b") 

For  the  ellipsoidal  cross  section  bar,  iig  = - ^  9"” 

4a  b 
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Correction  for  Young’s  modulus 

Young’s  modulus  can  obtained  either  from  flexural  mode  or  from 
longitudinal  mode.  The  corrected  Young’s  modulus  from  flexural 
mode  is 


E  = 


n 


where  Qt  is  the  correction  factor. 

For  rectangular  cross  section  bar  of  sides  2a  (along  x  direction) 
and  2^(=/2),  we  have  and  the  correction  factor  Is, 


SI. 


=  1 +6.585(1 +  0.0752V +  0.8 - 0.868^^ j 


8340(1  +  03023V  +  1713v^  j 


1  +  6338(1  +  0.1408  Iv  +  1336v-  / 


As  an  example  if  -=o.i,  r»L07,  which  represents  7%  correction. 
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Wave  equation  (w/  Saint  Venant  warping) 

For  the  torsional  mode,  the  wave  equation  becomes, 
ct^  2 

Rs  az2  ■  at^  ^  az^at^ 


One  can  see  that  there  are  4“’  order  derivative  terms  in  the  wave 
equation  when  coupling  Is  considered. 


^  Is  the  angular  displacement  of  a  “disc”  element  at  location  z . 
2 

and  |LX  are  from  the  coupling  effect 


2_/J^^dxdy 
^  I, 


J  is  the  polar  moment  of  inertia  of  the  bar’s  cross  section, 

Ip  Is  the  polar  moment  of  Inertia  of  the  bar’s  cross  section  when 
the  Saint  Venant  warping  function  0(x,  y)  is  zero. 


J  =  JJ 


r 


V 


2  2  30 

X  +y  +x— — y 

ay 


Ip=JJ(x  +y  )dxdy. 


The  solution  to  the  wave  equation  takes  on  the  wave  number 


2 
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Correction  for  transducer  mass 

No  coupling  case  (just  added  mass) 

If  the  mass  of  transducers  at  the  ends  Is  considered,  but  the 
coupling  motion  among  different  dimensions  Is  neglected,  the 
formulas  are  the  same  except  that  the  length  of  the  bar  Is 
replaced  by  an  effective  length.  For  longitudinal  and  flexural 
mode,  the  effective  length  Is 

Leff  -Leff  =L(l+m/M). 

For  torsional  mode,  the  effective  length  Is 

=L(l  +  2m/M) 

(After  Garrett’s  JASA  article  circa  1990) 

The  correction  can  be  obtained  by  means  of  Rayleigh’s  energy 
method.  The  effects  of  the  transducer  Is  only  the  effective 
Increase  of  the  length  of  the  bar  due  to  the  Increase  of  the  mass 
of  the  bar.  The  resonant  frequencies  Is  less  than  that  of  the 
original  bar.  Because  the  mass  increase  Is  may  or  may  not  be  a 
significant  correction  depending  on  experimental  relative  masses 
of  the  transducer  and  bar. 
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ACOUSTIC  RESONANCE  TECHNIQUES  APPLIED  TO  METROLOGY 
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Why  temperature  scales? 

Why  thermodynamic  temperature  scales? 
What  good  is  Boltzmann's  constant 
What  good  is  universal  gas  constant  /?? 
What  is  acoustic  gas  thermometry? 
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SccftfJ 


/'lyfv*.  Ken  a  V*:  I  eyh  /«, 
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Lor-l'Cn^  • 


oi-  Efrch'^oL 


Q=  Az(e-r^) 


Coir^^'f  P'f^c  ^  ^ 


7" 

•  C»LC 


Hot 
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Examples  of  Primary  Thermometers 


Dilute  Gas 
(equation  of  state) 


PV  =  nRT{\  +  — 

V 


Dilute  Gas 
(speed  of  sound) 


Mc^ 


RT(l  +  AjP  +  .  .  . 


Black  Body  Radiation  p 
(power  /  area)  a 


60h^c 


Johnson  Noise 
(power) 


p  = 


r 


u^rA/ 


Dipole  in  Field 
(energy) 


- 

ctnh 

/iB  1 

_ 

pB 

^  } 
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Acoustic  Gas  Thermometry 

dilute  monatomic  gas 


T  ^  T  ^  c\T) 
273.16  K  '  c\T^) 

Resonance  Method 


^  =  ^aI‘b 


3 

5  T 
3Mc^ 
5  T 


JL  =  [  =  L  ^  ^v(T)Y'^  (mY 

I  yw^fw]  i  }  \fw] 


469 


TRANSPARENCY  6 


RADIAL  RBSONANCES  IN  A  sRh^re 


N/aHE'5T  POSSIBLE  Q  AT  LOiu 

pEf^^fTY  (f^o  Vlsc-OOi 

A/OT  S^Ni1T(i/£  TO  SM/iUL  ERHORS 

iAj  cohJSTRucjfoh/  —  Nofj~oe<se-Aj£fi,Are 
HoDe*> 
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Clean  Gas 
and  Vacuum 


Valve 

Operator 


Figure  5.  Cross-section  of  resonator  and  pressure  vessel.  The 
transducer  assemblies  are  indicated  by  “T,”  and  the  locations  of 
the  capsule  thermometers  are  indicated  by  “PRT.” 
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Amplitude  Deviation  Amplitu 
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L  orchtz/Qi^  Liy)< 
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Simplest  Theory 


Rigid  spherical  shell,  ignore  temperature 


-  nth  spherical  Bessel  function 


dj„{pa) 

da 


2-Kf 

c 


(Rayleigh) 

/  =  ^  of  root;  /  =  0,  1,  2,  .  .  . 

/I  =  #  of  Bessel  faction;  n  =  0,  1,  2,  .  .  . 


476 


Figure  20.  Excess  half-widths  of  (0,/?)  resonances  with  argon  in 
the  resonator  scaled  by  lOVfrequency.  Ag  =  measured  g  minus 
calculated  g. 
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2.  0.t  0.6 

ressure  {m 
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f  2- (t/f 

1/p  =  7. 7?r,  /o  yof-^ .  •  • 

Arj  OH 

r  -  5"  Jr/O"  ^ 


PRESSURE  (kPa) 


sr.^fo-r  Tf.eon^  cee-f 
5:5t+o.G'4  D^4o  h-o.' 
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From  table  9  one 
can  conclude,  for  example,  that  1  ppm  of  water  in 
helium  will  decrease  Cq  by  3.93  ppm  and  that  1  ppm 
of  water  in  argon  will  increase  Cq  by  0.12  ppm. 


Table  9.  Sensitivity  of  co  to  impurities 


Impurity 

M 

(g/mol) 

To 

I  d{ci) 

Co  dx 

in  He  in  Ar 

H2 

2 

1.4 

0.23 

0.68 

He 

4 

5/3 

0.9 

H2O 

18 

1.32 

-3.93 

0.12 

Ne 

20 

5/3 

-4.0 

0.5 

N2 

28 

1.4 

-6.27 

0.03 

O2 

32 

i.4 

-7.3 

-0.07 

Ar 

40 

5/3 

-9.0 

CO2 

44 

1.4 

-10.3 

-0.37 

Kr 

84 

5/3 

-20.0 

-1.1 

Xe 

131 

5/3 

-31.8 

-2.3 

Hg 

201 

5/3 

-49.0 

-4.0 
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Coojc  ^  /f^/ 
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Table  1.  One-sigma  uncertainties  (in  parts  per  million)  from  var¬ 
ious  sources  in  the  redetermination  of  R 


I 

(Volume)^^^ 

density  of  mercury  at  20  “C 

0.28 

storage  and  handling  of  mercury 

0.20 

thermal  expansion  of  mercury  (0-20  “C) 

0.67 

random  error  of  volume  measurements 
corrections  from  weighing  configuration 

0.20 

to  acoustics  configuration 

0.10 

mass  of  counterweights 

0.14 

II 

Temperature 

random  error  of  calibrations 

0.8 

temperature  gradient 

0.4 

III 

M/yo 

Ar-40  standard 

0.7 

comparison  of  working  gas  to  Ar-40 

0.4 

IV 

Zero-pressure  limit  of  (fo„/vo„y 

s.d.  of  Co  from  70  observations  at  14  pressures 
thermal  boundary  layer  correction  (0.3%  of 

0.68 

thermal  conductivity) 

0.30 

possible  error  in  location  of  transducers 

0.55 

Square  root  of  the  sum  of  the  squares 

1.7 
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Measurement  of  the  ratio  of  the  speed  of  sound  to  the  speed  of  light 

James  B.  Mehl 

Physics  Department,  University  of  Delaware,  Newark,  Delaware  19716 
Michael  R.  Moldover 

Thermophysics  Division,  National  Bureau  of  Standards,  Gaithersburg,  Maryland  20899 

(Received  2  June  1986) 

Measurements  of  the  resonance  frequencies  of  the  acoustic  modes  and  of  the  microwave  modes  of 
a  single  cavity  can  determine  u  /c,  the  ratio  of  the  speed  of  sound  of  a  gas  to  the  speed  of  light. 
Such  measurements  with  a  monatomic  gas  would  determine  the  thermodynamic  temperature  T  with 
unprecedented  accuracy.  By  judicious  choices  of  cavity  geometry  and  resonance  modes,  u  /c  can  be 
measured  to  part-per-million  accuracy  using  cavities  whose  geometry  is  known  only  to  parts  per 
thousand.  These  techniques  can  also  be  applied  to  measurements  of  the  universal  gas  constant  R. 
A  measurement  of  R  would  also  require  an  accurate  determination  of  the  average  atomic  mass  of 
the  monatomic  gas. 
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AMPLITUDE  DEVIATION 
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FREQUENCY  (MHz) 


1473.5  1474.0  1474.5 

FREQUENCY  (MHz) 
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Table  iol:  The  difference  T- 


T-TJmK 
(isotherm  fits) 

T-TJmK 
(surface  fit) 

Argon 

302.92 

S.Og  ±  0.62 

4.4g  ±  0.3g 

293.13 

2.79  *  0.44 

3.0g  ±  0.2g 

253.15 

-2.83  *  0.53 

-2.5^  ±  0.22 

234.32 

— 3.6y  i  0,5g 

-Z.lj  ±  0.22 

217.10 

-3.95  ±  0.4g 

-3.7g  ±  0.23 

Xenon* 

302.9166 

3.52  ±  0.9g 

3.53  ±  0.8g 

^Errors  for  xenon  include  all  of  the  items  listed  in  Table  1  together  with  an  allowance  c 
0.43  mK  to  account  for  the  effects  of  the  virtual  leak  and  an  allowance  of  0.6  mK  to 
account  for  uncertainty  in  the  values  of  i43(  7). 
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Sonic  Gas  Analyzers 


Steven  L.  Garrett  and  Matt  Golden 
Penn  State  University 
and 

Robert  Keolian 
Naval  Postgraduate  School 


OUTLINE 


Non-specific  gas  analysis 
Direct  applications 
Secondary  applications  (GC) 

Historical  Perspective 
The  Schlagwetter-pfeife 
The  rise  of  the  Wheatstone  Bridge 

Resonator  design 

Isolate  a  single  mode,  suppress  overtones 


ARL  -  Penn  State 


TRANSPARENCY  2 


Sensor  Applications 


•  Mine  Safety 

Mine  sealed  with  148  still  inside 


BY  NIKOS  KONSTANPARAS 
Associated  Press 

KOZLU,  Turkey  —  Despite  ap¬ 
peals  from  weeping  family 
members  waiting  at  the  surface, 
officials  sealed  off  the  Koziu  coal 
mine  yesterday,  ending  all  hope 
for  148  miners  trapped  one-third 
of  a  mile  underground. 

In  Turkey  s  worst  mine  disaster. 


muiation  of  poisonous  gas  made  it 
perilous  to  enter  the  shaft. 

"No  one  down  there  is  alive. 
Yesterday,  they  brought  some  in* 
jured,  but  not  today,"  said  medic 
Sabahattin  Oztas  as  he  watched  a 
truck  dump  tons  of  dirt  on  the 
steel-covered  pit. 

Miners  took  turns  shoveling  the 
dirt  in  order  to  starve  the  mine  of 
oxygen. 


•  Agricultural  Chemical 


Yesterday,  headlines  in  Turkish 
newspapers  said  the  disaster  could 
have  been  averted.  The  stories 
said  computers  had  showed  in¬ 
creasing  accumulation  of  methane 
gas  at  the  pic,  but  that  manage¬ 
ment  failed  to  evacuate  the  mine. 

State  Minister  Omer  Barutcu 
said  an  eariy-waming  system  did 
not  work  because  methane  gas  in¬ 
creased  to  the  dangerous  level  just 


Methyl  bromide  warnings  requested 


BYBgrSYLOftOW 
Hpraa  Mhiv 

A  lecal-aid  group  for 
ranmmricen  has  sent  a  ietier  lo 
Central  Coast  growers 

asking  them  to  warn  neichbon  21 
days  before  ihev  use  the  powerful 
soti  fumigant  metlryf  bromide. 

Growers  are  not  obligated  le* 
gaify  to  provide  sudi  warnings, 
said  Michael  Meoter,  a  lawyer  in 
(he  Salinas  office  of  CaJifomta 
Rural  Legal  Asusunoe. 


quent  (he  area  is  pan  of  being  a 
good  netchbor.  Meuter  said. 

“I  think  if  grawets  want  to  no¬ 
tify  ihcir  nerghtaors  and  people 
(hat  are  nearby  about  ihetr  appli- 
cations,  they  are  free  to  do  that." 
said  Richard  Nutter,  Monterey 
County  agnculturai  commissioner. 
‘'But  there  is  no  legal  reputremem 
that  would  make  them  respnnsMe 
for  doing  that." 

BccauM  gromeis  would  want 
(heir  neighbors  valontanly.  the 
Agricultural  Commissioner's 


WNhout  Hlnn.  a  group  oppoeed 
to  methyl  bromide  use.  "People 
have  the  right  to  know  that  tosie 
cbemictls  are  heng  used  in  their 
back  yard." 

The  CRLA  sent  the  letters  to 
375  Central  Coast  strawbeiry 
growers  who  have  permiti  to  me 
methyl  bromide.  AJkmf  with  the 
ietteis.  gfiwers  received  a  sug¬ 
gested  warnmg  that  they  can  use 
or  modifjr. 

inchtd^  in  the  warning  are; 

■  A  list  of  symptoms  for  ihoK 


local  and  state  agencies  that  han¬ 
dle  pesticide  emergencies. 

Growers  also  received  a 
wammg-aooe  table  for  methyl  bro¬ 
mide  developed  by  the  Caliiormi 
Oepartmem  of  Pesticide  Regiala- 
(ioa:  the  table  is  not  used  by  rego- 
iaiors  under  cuiTenf  law.  For  es- 
ample.  if  a  grower  creaa  10  acres 
at  a  rate  of  250  pounds  per  acre, 
the  table  recommendi  that  anyone 
who  frepuents  the  sone  wiclun  a  1- 
mde  radun  should  be  warned, 
Meoter  said  that  under  Califor- 


Pete  Wilson's  office.  F 
advocates,  environmer 
othen  who  oppose  r 
mide  argue  the  exe 
growers  a  based  o* 
rather  than  soenufic  r< 
‘That's  probsblv  tn 
said.  "But  (the  dual  ; 
the  reality  of  things.  ‘ 

As  for  the  waniine 
Nutter  saidL  T  don  i  k 
information  that  would 
justify  notification  wtthi 
BcCTuie  It  poses  h 
USB  of  methyl  bromide 


•  Other  obvious  applications 


Electric  vehicles  (H2) 

Chemical  synthesis  {e.g.,  H2  in  NH3  or  Cy 
Blast  furnaces  and  metals  processing  (H2) 


Ethylene  oxide  in  medical  equip,  sterilization 


ARL  -  Penn  State 


503 


TRANSPARENCY  3 


S.  Garrett  Sonic  Gas  Analysis 

Sound  Speed  in  Gas  Mixtures 

•  Equipartition  Theorem 

Each  quadratic  degree  of  freedom  gets  kBT/2 
Communism  is  not  dead! 

Temperature  is  related  to  kinetic  energy 

Pythagorean  sum;  <v^>=-  <Vx^>  +  <w^>  +  <Vz^> 
Symmetry:  <v^>  =  <Vy^>=  <w^>  =  (l/2)kBT/m 


Adiabatic  Sound  Speed 


a^-y 


M 


Sound  Speed  in  Gas  Mixtures 
Mean  molecular  mass 

Mn„x  =  xMi  +(1-x)M2 
Mean  polytropic  coefficient 

Tmix- - - - 7 - 1 

XCyi+{l-x)Cy2  ^ 

Approximate  sound  speed  ratio 

(J?  \+[{M^-M^)l  M^]x 

Sensitivity  analysis  for  resonance  freq.  shift,  5x 


£ 

f\ 


Sa. 


^  mix 


a, 


(r^-Yx) 


Yi 


M, 


Sk 
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(lorreii.  ( lolcJcn  (K-  Kealian  Sonic  Gas  Analysis 

Sensitivity  in  Air 


cir*SMUn 


Approximate  sensitivity  factor 

Assume  y-factor  is  small 


y3= 


/ 


Jt,p 


2M, 


Results  for  various  gases  in  air 


0  100  200 


ARL  -  Penn  State 


Molecular  Weight  (a.m.u.) 
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S.  Garrett  Sonic  Gas  Analysis 

Schlaswetter-pfeife 

•  German  mine  safety  whistle  ( 1 884) 

Pump  mine  exhaust  through  a  whistle 
Wavelength,  X  =  2L  for  an  open-open  pipe 
Frequency,  fx-sJX-  a/2L 

•  Methane  frequency  shift 

Molecular  weights 

^methane  ~  Mnr  —  16.043-28.964  —  _Q  /j  /j  ^ 

28.964 

Polytropic  coefficients 

y methane  Y air  —  t3 1  1.403  __  06^^ 

“[403 

Relative  frequency  shift 

^  =  0.3  8^ = 0.1 9  (& 

J  cur  ^air  ^ 

Example:  1%  Methane  (x  =  0.01) 

Assume  =  440  Hz,  L  =  Ay2  =  39.2  cm,  5f  =  0.84  Hz 

What  is  X  for  an  explosive  mixture? 

•  Temperature  compensation  and  “detection” 

Use  two  whistles  in  good  thermal  contact 
Listen  for  beat  frequency  fbeat=  fmine  -  fair 
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S.  Garrett  Sonic  Gas  Analysis 

(Aiinost)Thermal  Conductivity  Theory 


Thermal  Conductivity 
Simple  kinetic  theory  model 


K 


Clo^ 


Also  related  to  thermal  molecular  speed 
Square-root  absolute  temperature  dependence 
Pressure  independence  at  “normal”  pressures 

Complex  combination  of  constants 

a  =  Characteristic  molecular  diameter  (Angstroms) 
Q  =  Collision  function  of  keT/s  (tough  calculation) 

Ko  =  1.9891  X  10-^  cal  gm’^  seC’ 


•  Requires  off-equilibrium  measurement 

Extremely  dependent  upon  experimental  conditions 
“Probably  the  principal  limitation  upon  the  usefulness  of 
the  method  is  the  fact  that  one  can  not  calculate  the  significance  of 
the  indications  obtained,  but  must,  in  every  case,  depend  upon 
empirical  calibration  under  conditions  closely  approximating  those 
of  actual  service.” 

Nat.  Bur.  Standards,  Tech.  Paper  No.  249  (1924) 

Irreversible  transport  process 

No  resonant  enhancement,  Q  =  1/2 
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5.  Garrett  Sonic  Gas  Analysis 

Gas  Thermal  Conductivity  Sensor 


•  “Hot  wire”  cell  at  right 
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S.  Garrett  Sonic  Gas  Analysis 

Resonant  Sound  Speed  Measurement 

•  Complete  system  for  Helium  Recovery  Air  Alarm 
Garrett,  etaL  Physica  107B,  601  (1981) 


To  Recovery  Line 


Precision/Cost  ratio 

For  frequency,  10  ppm  a  $20,  based  on  quartz  watch 
5-digit  voltmeter  a  $1,000 
Availability  of  cheap  transducers 
Electret  microphone  a  $0.59  (lOOO’s) 

Earphone  sound  source  a  $0.20  (lOO’s) 

Availability  of  cheap  circuit  components 
Voltage  controlled  oscillators  a  $1.00  (lOOO’s) 
Multiplier  chips  {e.g.,  AD633)  a  $5.00  (lOO’s) 
Digital  output  (fi-equency) 

More  compatible  with  analysis  and  display  hardware 
No  A/D  conversion  required 
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A.-TO  ^ 

a 

C>U«S<^^  H  €.«  v»^  >' 


V.OS  S- -4sj  ^ 


<-  »^py^ 


fjN\=;S‘oK^ 


i««o 


*•1— 


S  S.'' ??■  <v«.»s,  «'»W  r, 


2.,  <t\w^S 

^T 


♦.f  VM\\ 


3.  F  '«  ~v«n<A\p^<vA.^«v^  «' 
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S.  Garrett  Sonic  Gas  Analysis 


Novel  Resonator  Geometries 

•  Basic  Dual  Half-Wavelength  Dipole  Design 


Does  not  drive  even  modes 

Mike  located  at  pressure  node  for  S^^ode 

Dipole  drive  localizes  flow  at  edge 

Dipole  gas  motion  is  below  cut-off  for  open  tube 

Dipole  mike  is  insensitive  to  incident  plane  waves 
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51  Garrett  Sonic  Gas  Analysis 


Other  Enhancements 

•  Short  Screened  End 


26  21  22  23  24  25  26  26  21  22  23  24  25  26 

25  Electronic  thermometer  (e.g.,  thermistor) 

26  Screen  to  prohibit  physical  frequency  modification 

•  Quasi-Helmholtz 
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S.  Garrett  Sonic  Gas  Analysis 

Summary 

•  Advantages  of  Sonic  over  TC  gas  analyzers 
Resonant  signal  enhancement 
Quality  factor,  10  >  Q  >  1000 
Calculable  response  independent  of  “drive  level” 
Gas  flow  rate  is  not  a  problem 
Recent  availability  of  new  components/transducers 
High  precision-to-cost  ratio 
Intrinsically  digital  (frequency)  output 
Expanded  processing  and  telemetry  options 
Infinite  lifetime 

No  element  bum-out  or  chemical  regeneration 
Pressure  independent 

Second  virial  coefficient  gives  38  ppm/psi 
Also  advantage  of  thermal  conductivity  analyzers 

•  Disadvantages  of  Sonic  (also  shared  by  TC) 
Non-specific 

A  problem  if  there  are  many  possible  “contaminants** 
Humidity  can  be  a  problem 
Requires  active  temperature  compensation 
Much  easier  for  sonic,  but  still  required 
Humidity  changes  effect  frequency 

±20%  change  in  RH  give  ±800  ppm  in  frequency 
May  require  humidity  compensation 
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Measurements  of  Relaxation  Processes 
in  Gases  and  Liquids 
Using  Resonance  Systems 

by 

Henry  E.  Bass 


Jamie  L.  Whitten  National  Center  for  Physical  Acoustics 
The  University  of  Mississippi 
University,  MS  38677 


The 

lUniversUy  of  Mississippi 
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Absoq)tion  and  Dispersion  Due  to 
Microscopic  Processes 


(uA>oF 


J,  frequency  (^iz) 


related  to  relaxation  or  reaction  time/rate 
Magnitude  of  dispersion/absorption  related  to  concentration, 
heat  of  reaction 

The 

Of  Mississippi 
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Dispersion  and  Absorption 
Can  be  Small 

—  for  N2  -  IxlO'"^ 

V 


aforair@lkHz=2xlO  ^  dB/meter 
for  seawater®  10  kHz  ~  10"Mb/  meter 
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Gases 


Liquids 


Resonators 

Evans/Parker 
Cylindrical  Tube 

Rectangular  Room 

Fisher 

Cylindrical  Tube 
Spherical  Resonator 

The 

of  Mtsstssippi 
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Definition  of  Terms 

a  Absorption/unit  length 

5=ac/7i  With  of  Resonance  Peak 

(Half  power  points) 

f 

Q  =  -^  Quality  of  resonance 

2o 

T  Time  required  for  Energy  Density 

to  decay  oOdB  = 

Reverberation  Time 

The 

University  of  Mississippi 
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Boundary  Effects 

Absorption  of  Sound  in  Tube 
Wall  Loss  Proportional  to  j_ 

End  Loss  Few  %  of  Wall  Loss 
5  =  BVn 

Other  Sources  of  Bulk  Absorption 
Viscosity  and  Thermal  Conduction 
5^/  Typical  for  Tubes  -  Larger  for  Sptoe 


The 

m  University  of  Mississippi 
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Interpretation  of  Results 

6  =  5(o  +  5v 

s,=(RCi/c;c) 

•  Determine  8oj  from  pure  measurements 

•  Determine  from  5^ 

fo=  1240HzforX=1.92% 

Collision  number  2  X  10^ 

of  Mississippi 
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Absorption  coefficient  (a  dB/ktn) 


TRANSPARENCY  10 


Frequency  (Hz) 


University  of  Mississippi 
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Design  Considerations  for  Thermoacoustic  Resonators 

Thomas  J.  Hofler,  Physics  Department,  Naval  Postgraduate  School 


Presentation  Outline 

•  Resonator  Topology  &  Efficient  TA  Refrigerators 

-  4  resonator  topologies 

-  Some  theory 

-  Optimal  dimensions 

•  Suppressing  Harmonic  Distortion  in  TA  Prime  Movers 

-  An  example  prime  mover  apparatus 

-  Suppression  criteria  &  numerical  determination 

•  Amplitude  Stability  Issues  in  a  Novel  TA  Prime  Mover 

-  Thermal  resonator  tuning 

-  Modest  amplitude  oscillations 

-  Sudden  large  amplitude  changes 
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Design  Considerations  for  Thermoacoustic  Resonators 

Resonator  Topology  &  Efficient  TA  Refngerators 


4  Thermoacoustic  Refrigerator  Resonator  Topologies 


— - - yj2  - - 

- c? 

b 

■= 

1 

[—  cold  portion 

TRANSPARENCY  3 


Match  p  at  transition  A, 
P20  _ 

P^o 
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Design  Considerations  for  Thermoacoustic  Resonators 

Resonator  Topology  &  Efficient  TA  Refrigerators 


Sample  Loss  Calculations 


Two  examples  at  the  extremes  of  the  resonator  spectrum 


Example  #1:  XI a  «  0.4  radian;  Prandtl  =  2/3  Example  #2:  Xia  *  0.1  radian;  Prandtl «  0.2 


02/01  02/DI 
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Design  Considerations  for  Thermoacoustic  Resonators 

Resonator  Topology  &  Efficient  TA  Refrigerators 


Optimal  diameter  change  &  resulting  reduction  in  loss 


Diameter  ratio  for  minimum  02  loss  Lots  reduction  at  optimum  02/DI 
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Design  Considerations  for  Thermoacoustic  Resonators 

Suppressing  Harmonic  Distortion  in  TA  Prime  Movers 

«i  Impedance  &  work  flow  at  cold  termination 


High-Amplitude  Prime  Mover  Apparatus 


*  work 
-  !2I 
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Half-wave  TADTAR  Design 

Advantages 

•  High  amplitudes  possible  with  thermoacoustic  driver 

•  Proper  orientation  of  heat  &  work  flow 

-  Stirling  cycle  transport  improves  performance 

•  Good  optimization  freedom 

-  Diameter  &  power  of  TAD  &  TAR  are  independent 

-  Stack  locations  for  TAD  &  TAR  are  independent 

•  Performance  is  very  sensitive  to  location  of 
Pressure  AntiNode  (PAN) 

-  Can  be  thermally  tuned  for  easy  starting 

Disadvantages 

•  Performance  is  very  sensitive  to  location  of  PAN 

-  Will  the  acoustic  amplitude  be  stable? 

•  Unknown  heat  leak  from  hot  heat  exchanger 


Ambient- 


gwiQt 


Ambient— t 

Cold -I 


MWiQf 
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Design  Considerations  for  Thermoacoustic  Resonators 

Amplitude  Stability  in  a  Novel  Prime  Mover 


Preliminary  Measurements 

•  Run  driver  only,  with  cooler  sections  empty 

•  Measure  acoustic  amplitude  &  maintain  400°  C 
hot  temperature  via  human  controller 

•  Resonator  tuning 

-  Gross  tuning  with  tube  inserts 

-  Fine  tuning  with  heater  tape  on  upper 
tube  and  heat  gun  on  lower  tube 

-  3  thermocouples  were  attached  to  each  of 
the  two  small  diam.  tubes.  Each  triplet  is 
averaged,  converted  to  Kelvin,  &  the 
upper-to-lower  ratio  calculated.  This 
“Temp.  Ratio”  indicates  the  level  of  fine 
tuning. 
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Design  Considerations  for  Thermoacoustic  Resonators 

Amplitude  Stability  in  a  Novel  Prime  Mover 


0 


53 


Resonator  Smalt>Tube  Temperature  Ratio 


dynamic  pressure  (V  rms) 
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Design  Considerations  for  Thermoacoustic  Resonators 

Amplitude  Stability  in  a  Novel  Prime  Mover 

Oscillatory  TADTAR  amplitude  (rms) 

Period  =  1  sec.  approx. 


0.20 


0.15 


0.10 


0.05 


200  400  600  800 
sample  # 
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Design  Considerations  for  Thermoacoustic  Resonators 

Amplitude  Stability  in  a  Novel  Prime  Mover 

Sudden  increase  in  TADTAR  amplitude 
sample  interval  =  0.4  ms 
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ACOUSTIC  IMEASUREIMEWTS  of  THERIMODYNAWIIC 
and  TRANSPORT  PROPERTIES  in  GASES 


1 1  United  States  Department  of  Commerce 

I  I  National  Institute  of  Standards  and  Technology 


MICHAEL  R.  MOLDOVER 

NIST  Fellow 
Physical  and  Chemical 
Properties  Division 
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THERMODYMAMIC  PROPERTIES 


1.  Ideal-gas  heat-capacity:  C^[T] 

2.  Virial  equation  of  state; 

PV=Rt{i+  B(T)P  -h  C(T)fP'  +  D{T)ff  +  ... ) 

Note:  Resonance  techniques  are  not  recommended 
for  liquids  because  oscillations  of  container  cannot 
cannot  be  separated  from  oscillations  of  fluid 
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TRANSMISSION 


TRANSPARENCY  7 


compound 


temperature 

range 

K 


maximum 

pressure 

kPa 


number  of 
isotherms 


Candidate  refrigerants  for  vapor  compression  cycles 


CF3-CHF2 

R134a 

233  -  340 

600 

10 

CFC12-CH3 

R141b 

260  -  315 

70 

5 

CHC12-CF3 

R123 

260  -  335 

80 

6 

CHFC1-CF2-C1 

R123a 

265  -  300 

50 

2 

CHF2-0-CHF2 

E134 

255  -  327 

170 

6 

CHF2-0-CHF2 

E134 

255  -  374 

90 

8 

CF2-0-CH2-CF3 

E245 

278  -  384 

50 

5 

CF3-HF-CHF2 

CHFC1-CF3 

R236ea 

267  -  380 

600 

8 

R124 

250  -  400 

900 

17 

CHF2-CF3 

CHF2-CH3 

R125 

240  -  400 

1,000 

9 

R152a 

240  -  400 

1,000 

9 

CF3-CH3 

R143a 

240  -  400 

1,000 

9 

CF3-CH2-CF3 

R236fa 

276  -  400 

1,000 

7 

CHF2-CF2-CH2F 

R245ca 

311  -  400 

900 

5 

CF3-0-CF2H 

E125 

260  -  400 

1,000 

13 

1  composition 

R134a/R32/R125 

260  -  400 

1,000 

12 

CF3-CF2-CF2-CH2F 

R338mccq  300  -  400 

Thermoacoustic  Refrigeration 

400 

6 

5  compositions 

He/Xe 

210  -  400 

1,500 

42 

Semiconductor  Processing 

SFc 

230  -  460 

1,500 

16 

CF. 

300  -  475 

1,500 

9 

C2F6 

210  -  475 

1,500 

14 

CI2 

planned 

HBr 

planned 

BCI3 

planned 

WF, 

planned 
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M  \  RT  RT  RT  RT 

First,  consider  limit  asp  -  0 

C^iT)  ^  f(T)  ^  Mu^ 

R  /(T)  -  1  Mu^  -  RT 

For  polyatomic  gases  /  1. 

(e.g.  SFg  ,  f(29S  K)  =  1.094;  C^/R  =  1/10.6) 
Uncertainty  of  M  from  impurities  amplified  10.6  x 

e 

Ansatz:  B(T)  =  ^^[1  +  (r^  -  l)(g  -  1)] 
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Speed  of  sound  data  for  CHFj-CFj. 
Pentafluoroethane  =  refrigerant  R 125 

Temperature  range  0.71  ^  TIT^  s  1.18 


0  0.2  0.4  0.6  0.8  1.0 

p,  MPa 


Deviations  of  data  from  u(P,  71  surface. 
RMS  deviations  are  2  x  1 0  ®  x  u. 

Impurities  cause  larger  systematic  errors. 
The  surface  uses  3  parameters  for  C°{T), 
(ideal-gas  heat-capacity),  3  parameters 
for  each  of  the  virial  coefficients  5(71  and 
C(7)  and  2  parameters  for  D{T). 
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102aC°/C 
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Ideal  gas  heat  capacity  of  R125, 
deduced  from  the  speed  of  sound  data 
through  the  relation: 

Cp°(7)//?  =  Mu/l{Muo^- RT). 

R  is  the  universal  gas  constant;  M  is  the 
molecular  weight  and  is  the  zero- 
pressure  limit  of  u^[P,  71 


Deviations  of  C°[T)IR  from  a  polynomial, 
o  present  data 
■  independent  spectroscopic 
information  (estimated) 
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=  iiu  d)/d  d/d  V  eOt 
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Compression  factor  of  R125  deduced 
from  ulP,T)  data  using  the  virial  equation. 
The  viriai  coefficients  were  assumed  to 
have  the  temperature  dependencies  of  a 
hard-core  square-well  intermolecular 
potential  and  the  parameters  in  the 
potentials  were  fitted  to  the  u{P,T)  data. 


Fractional  deviations  of  the  density  of 
R125  deduced  from  uiP,T)  data  from 
independent  density  measurements  by 
Boyes  and  Weber. 
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TRANSPORT  PROPERTIES 

1.  Greenspan  acoustic  viscometer  for/7; 


hard  sphere  of  diameter  o 


2.  Prandtl  number  machine 


Pr  = 


^  _  viscous  diffusivity 
A  thermal  diffusivity 


=  —  for  hard  sphere 
3 


3.  If  time  permits,  electromagnetic  equivalent  of 
Greenspan  viscometer:  reentrant  resonator 
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Loss  mechanisms  spatially  separated  in  standing  waves 
High  thermal  losses  High  thermal  losses 


High  viscous  losses 


Greenspan  viscometer:  A  >  length 
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Double  Helmholtz  Resonator 


"2V 
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z  (mm) 
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Acoustic  pressure  contours  near  duct  end 

Calculated  from  solutions  of 


(V2  +  A;Xr)  =  0 


r(mm) 
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TRANSPARENCY  20 


'  o  viscometer  0 
j  •  viscometer  1 
*  viscometer  2 
A  viscometer  5 
A  viscometer  13 
o  viscometer  23 


O.ooO 

o.oot  i-o.oo't 
O,ooz±-O,oo4 


A  Argon,  viscometer  23 
•  Argon,  viscometer  13 
^Helium,  viscometer  13 
o  Propane,  viscometer  13 


m 
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Greenspan  viscometer  as  an  absolute  instrument 
(no  calibration) 

Argon  in  4  viscometers 

Experimental  77  compared  with  reference  values 


0.01 

?=■ 

^  0.00 


-0.01 


— 

-1 — J— 

□ 

-  0 

* 

0 

i 

0 

■ 

□ 

0 

I  a 

cx  . 

. 0. . 

. . Q . 

0 

.  t 

0 

0 

■ 

0 

■ 

■ 

i 

s 

□ 

• 

□ 

1 

1 

■ 

s 

• 

• 

■ 

• 

0 

• 

_i _ 

• 

t 

• 

_ 

100  1000 

p  (kPa) 


Viscometer  23  0.002  ±  0.003 


Viscometer  24  —0.005  ±  0.003 


Viscometer  12  0.001  ±  0.006 


Viscometer  13  —0.002  ±  0.003 
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Co-axial  cable 


■niermometer 
well 


Gold  o-ring 
Inductor 


f-^c 


0 

L 


50 

- 1 _ ^ _ I 

mm 


^  ~  lir  VrV 

C  ^  ffntJ  iiTi 
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Two-mode  reentrant  radio  frequency  resonator 
and  equivalent  circuit 
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Temperature  dependence  of  the  dielectric  constant  of  water 


T(K) 


•  our  data 

n  Fernandez  et  at,  Int.  J.  Thermophysics  929(1995)  with  LCR  meter 

O  Femdndez  et  at,  Int  J.  Thermophysics  1^  929(1995)  with  transformer  bridge 

A  Malmberg  &  Maryott,  J.  Research  Nat  Bur.  Stds.  5^  1(1956) 

V  Oshry,  PhD.  Thesis,  Univ.  of  Pittsburgh  (1949) 

S^,  =  296.3M5  -  1.219±0.005(T/K)  +(2.03±0.01)xl0-^(r/K)^ -  (1.31dO.01)xl0-‘(T/K)\  s  =0.007  (la) 
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